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ABSTRACT 

This  project  involved  high  Reynolds  number  testing  and  measurements  of  the  flow  around  a  two-dimensional 
hydrofoil  having  a  Navy-relevant  cross  section.  Three  multi-week  test  campaigns  in  the  US  Navy’s  William  B. 
Morgan  Large  Cavitation  Channel  were  completed.  These  experiments  extended  the  Reynolds  number  range  of 
prior  hydrofoil  tests  by  more  than  an  order  of  magnitude.  The  measurements  included  laser-Doppler 
velocimetry  profiles  and  spectra,  particle-imaging  velocimetry  flow  fields,  static  surface  pressure  distributions, 
dynamic  surface  pressure  fluctuations,  and  hydrofoil  vibration.  This  experimental  effort  emphasized  the  flow 
physics  leading  to  trailing-edge  vortex  shedding.  The  main  outcomes  of  this  work  are  a  complete  fluid 
mechancial  data  set  for  use  in  development,  testing,  and  validation  of  high  Reynolds  number  flow  models;  and 
insight  and  understanding  of  the  fluid  phenomena  that  lead  to  trailing  edge  vortex  shedding. 

RESULTS 

The  main  findings  of  this  project  are  contained  in  the  two  attached  manuscripts.  The  research  work  completed 
here  represents  the  doctoral  dissertation  of  Dwayne  A.  Bourgoyne,  and  the  master’s  thesis  of  Joshua  A.  Hamel. 
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At  high  Reynolds  number,  the  flow  of  an  incompressible  viscous  fluid  over  a  lifting  surface  is  a  rich  blend  of  fluid 
dynamic  phenomena.  Here,  boundary  layers  formed  at  the  leading  edge  develop  over  both  the  suction  and  pressure 
sides  of  the  lifting  surface,  transition  to  turbulence,  separate  near  the  foil’s  trailing  edge,  combine  in  the  near  wake, 
and  eventually  form  a  turbulent  far-field  wake.  The  individual  elements  of  this  process  have  been  the  subject  of 
much  prior  work.  However,  controlled  experimental  investigations  of  these  flow  phenomena  and  their  interaction  on 
a  lifting  surface  at  Reynolds  numbers  typical  of  heavy-lift  aircraft  wings  or  full-size  ship  propellers  (chord-based 
Reynolds  numbers,  Rec  ~107  —  10s)  are  largely  unavailable.  This  paper  presents  results  from  an  experimental  effort 
to  identify  and  measure  the  dominant  features  of  the  flow  over  a  two-dimensional  hydrofoil  at  nominal  Rec  values 
from  near  one  million  (1M)  to  more  than  50  million  (50M).  The  experiments  were  conducted  in  the  U.  S.  Navy's 
William  B.  Morgan  Large  Cavitation  Channel  with  a  solid-bronze  hydrofoil  (2.1  m  chord,  3.0  m  span,  17  cm 
maximum  thickness)  at  flow  speeds  from  0.25  to  18.3  m/s.  The  foil  section,  a  modified  NACA  16  with  a  pressure 
side  that  is  nearly  flat  and  a  suction  side  that  terminates  in  a  blunt  trailing-edge  bevel,  approximates  the  cross  section 
of  a  generic  naval  propeller  blade.  Time-averaged  flow  field  measurements  drawn  from  laser-Doppler  velocimetry, 
particle-imaging  velocimetry,  and  static  pressure  taps  were  made  for  two  trailing-edge  bevel  angles  (44°  and  56  ). 
These  velocity  and  pressure  measurements  were  concentrated  in  the  trailing  edge  and  near  wake  regions,  but  also 
include  flow  conditions  upstream  and  far  downstream  of  the  foil,  as  well  as  static  pressure  distributions  on  the  foil 
surface  and  test  section  walls.  Observed  Reynolds-number  variations  in  the  time-averaged  flow  over  the  foil  are 
traced  to  changes  in  suction-side  boundary  layer  transition  and  separation.  Observed  Reynolds-number  variations  in 
the  time-averaged  near  wake  suggest  significant  changes  occur  in  the  dynamic  flow  in  the  range  of  Rec  investigated. 
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1.  Introduction 

Relative  motion  between  a  submerged  body  and  a  viscous  incompressible  fluid  induces  hydrodynamic  lift 
and  drag.  The  shape  of  the  body  and  the  characteristics  of  its  flow  field  govern  the  distributions  of  surface  pressure 
and  shear  stress  that  lead  to  these  forces.  Typically,  the  steady  or  unsteady  components  of  these  hydrodynamic 
forces  ultimately  determine  the  utility  of  the  body  for  potential  applications.  Air-  and  water-borne  transportation 
systems  are  often  characterized  by  turbulent,  high-Reynolds  number  flows.  Thus,  experimental  research  in  this  area 
commonly  strives  to  deduce  scaling  laws  and  cause-and-effect  relationships,  while  computational  efforts  focus  on 
developing  and  validating  predictive  flow  models.  Unfortunately,  the  scarcity  of  experimental  studies  at  high 
Reynolds  numbers  hampers  all  of  these  approaches.  The  goal  of  this  paper  is  to  address  this  situation  through  the 
presentation  of  experimental  results  for  the  time-averaged  incompressible  flow  over  a  lifting  surface  at  Reynolds 
numbers  approaching  those  of  heavy-lift  aircraft  wings  and  full-scale  ship  propellers. 

The  surface-bounded,  near-wake,  and  far-wake  flows  formed  by  a  lifting  surface  are  phenomenologically 
interesting  and  have  practical  importance.  For  example,  surface  curvature  influences  boundary  layer  development, 
transition,  and  separation,  which  together  determine  the  initial  conditions  for  the  foil’s  wake.  Given  that  these 
phenomena  may  all  behave  differently  with  increasing  Reynolds  number,  any  overall  scaling  for  lift,  drag, 
hydroacoustic  noise,  or  other  foil-performance  measure  may  result  from  the  interplay  of  multiple  phenomena.  On 
the  practical  side,  high-Reynolds  number  lifting  surfaces  are  essential  for  flight,  maneuvering,  propulsion,  and 
control  of  air-  and  water-bome  vehicles.  In  such  applications,  the  lifting  surface  is  designed  to  meet  performance 
criteria.  Given  the  complexity  of  flow  phenomena  and  the  wide  range  of  length  and  time  scales  that  must  be 
managed  in  high-Reynolds  number  turbulence,  it  is  safe  to  say  that  the  development  of  the  needed  design  tools  is  not 
yet  complete.  For  lifting  surfaces  in  incompressible  flow,  this  situation  exists  in  part  because  of  the  lack  of 
controlled  experimental  data  at  full-scale  chord-based  Reynolds  numbers,  Rec  =  U0Clv  (where  U0  is  the  flow  speed 
far  upstream  of  the  foil,  C  is  the  foil  chord,  and  v  is  the  kinematic  viscosity  of  the  flowing  liquid).  Few  experimental 
results  exceed  Rec  of  several  million  (M)  while  many  applications  lie  in  the  Rec  range  of  30M  to  100M. 

This  paper  presents  the  results  of  an  experimental  study  of  the  flow  over  a  two-dimensional  hydrofoil  for 
Rec  values  from  near  1M  to  more  than  50M.  The  foil  has  a  nearly  flat  pressure  side  and  a  NACA-16  suction  side 
modified  with  a  rounded  trailing  edge  bevel  of  apex  angle  44°  (Geometry  I)  or  56°  (Geometry  II).  The  foil  section 
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is  typical  of  naval  propeller  blades  of  moderate  thickness  and  chamber,  and  was  chosen  for  its  application-relevance 
and  its  potential  for  Re-dependent  flow  features.  In  particular,  three-dimensional  propeller  design  techniques  often 
rely  on  two-dimensional  section  performance  data  (Kerwin  1986).  Also,  the  trailing-edge  bevel  —  typically 
introduced  on  propeller  blades  for  structural  integrity  during  severe  off-design  operations  -  leads  to  a  compact 
region  of  flow  separation  whose  characteristics  depend  on  bevel  geometry  and  Reynolds  number.  The  present 
experimental  effort  is  concentrated  on  this  near- trailing-edge  region,  and  the  results  reported  here  augment  the 
available  measurements  of  flat  plate  boundary  layer  flows  (Eaton  and  DeGraaff  2000,  Osterlund  et  al.  2000)  suitable 
to  test  high-Reynolds-number  exterior-flow  models. 

The  principal  goals  of  this  study  were  to  identify  phenomena  that  affect  foil  performance  at  high  Reynolds 
number,  and  to  provide  a  comprehensive  data  set  for  validation  of  high  Reynolds  number  numerical  models.  Here, 
foil  performance  includes  lift,  drag,  and  the  prevalence  of  near  wake  vortex  shedding,  an  important  hydroacoustic 
noise  source  for  non-cavitating  lifting  surfaces  (see  Blake  1986,  Blake  and  Gershfeld  1989,  or  Howe  1998).  The 
measurements  reported  here  document  only  the  time-averaged  flow  over  the  foil,  but  even  these  results  display 
interesting  Reynolds  number  dependences  resulting  from  the  interplay  of  boundary  layer  transition,  flow  separation, 
and  wake  formation.  Featured  prominently  in  this  interplay  is  the  phenomenon  of  near-wake  vortex  shedding, 
though  these  findings  are  beyond  the  scope  of  this  paper. 

Prior  results  documenting  the  mean  flow  at  low-Mach  number  on  foils,  struts,  and  turbine  blades  has  been 
largely  motivated  by  the  need  to  understand  and  predict  lifting-surface  performance,  and  possibly  eliminate  vortex 
shedding.  Early  work  in  this  area,  particularly  for  compact  bodies,  is  reviewed  in  Bearman  (1965);  a  more  recent 
review  is  provided  in  Blake  (1986).  Mean  velocity  profiles  in  the  boundary  layer  and  in  the  near  wake  of  minimally 
lifting  struts  are  provided  in  Blake  (1975)  for  Rec  of  ~  2M.  The  findings  of  the  present  effort  are  similar  to  these  but 
extend  the  range  of  experimental  Rec  to  more  than  50M  and  include  the  influence  of  lift.  More  recent  research  on 
individual  foils  and  struts  experiencing  steady  inflow  conditions  has  been  computational  (Wang  and  Moin  2000, 
2002),  but  again,  at  lower  Rec  than  reached  in  the  present  studies.  This  situation  is  similar  for  unsteady  foil-flow 
studies  (Luire  et  al.  1998,  Ho  and  Lakshminarayana  1997)  where  Recs  only  as  high  as  4M  have  been  investigated. 

Low-Mach  number  airfoil  studies  cover  many  of  the  same  phenomena  presented  in  this  paper.  Boundary 
layer  transition  on  airfoils  is  reviewed  in  Malick  (1997)  and  a  recent  prediction  technique  is  presented  in  Brodeur 
and  van  Dam  (2000).  Experimental  studies  of  steady  (Bastedo  and  Mueller  1986,  Fitzgerald  and  Mueller  1990)  and 
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unsteady  (Covert  and  Lorber  1932)  boundary  layer  separation  have  been  conducted  at  Rec  values  from  0.14M  to 
0.7M.  Mean  and  turbulent  flow  fields  in  the  vicinity  of  boundary  layer  separation  on  an  asymmetric  trailing  edge  in 
a  two-stream  wind  tunnel  are  reported  in  Thompson  and  Whitelaw  (1985).  Here,  at  Rec  ~  2M,  it  was  determined 
that  the  boundary  layer  flow  that  enters  the  separated  region  near  the  trailing  edge  significantly  influences  the 
development  of  the  wake.  The  present  study  shows  that  this  result  persists  at  higher  Rec*  A  general  review  of 
turbulent  boundary  layer  separation  is  provided  in  Simpson  (1989). 

Research  into  the  performance  of  compressor  and  turbine  blades  and  cascades  has  many  similarities  to  the 
current  study  as  well,  including  the  emphasis  on  lift,  drag,  boundary  layer  development,  trailing  edge  characteristics, 
and  near- wake  vortex  shedding.  A  review  of  this  literature  prior  to  1987  is  provided  in  a  series  of  papers  (Deutsch 
and  Zeirke  1987,  1988a,b)  where  experimental  results  at  Rec  =  0.5M  are  reported.  More  recent  experimental 
(Umbaldi  1996,  Halstead  et  al.  1997a)  and  computational  (Halstead  et  al.  1997b)  efforts  cover  single-  and  multi¬ 
stage  cascades  at  Rec’s  up  to  1.6M  and  0.6M,  respectively.  In  addition,  measurements  of  boundary  layer  transition 
and  separation  conducted  on  flat  plates  at  turbine-blade  conditions  are  reported  in  Volino  and  Hultgren  (2000)  for 
plate-length-based  Reynolds  numbers  up  to  0.3M. 

All  of  the  experimental  work  described  above  on  foils,  struts,  and  blades  involves  test  models  having  Rec 
values  at  or  below  4M.  Furthermore,  prior  steady-flow  results  at  or  above  Rec  =  2M  involve  minimally  lifting  struts 
with  section  lift  coefficients  that  fall  below  that  typically  sought  for  propeller  blades  or  aircraft  wings.  By 
comparison,  the  current  test  model  generated  significant  lift  forces,  and  represents  an  important  step  towards  the 
understanding  of  practical,  three-dimensional  propeller  blades.  Moreover,  the  Rec  values  in  the  current  experiments 
are  high  enough  to  comfortably  span  the  gap  between  the  prior  studies  and  many  full-scale  applications. 

The  remainder  of  this  paper  presents  our  results,  observations,  and  analysis  of  the  time-averaged  low-Mach 
number  flow  over  a  two-dimensional  hydrofoil  with  two  different  trailing  edges.  The  next  section  describes  the  test 
model,  flow  facility,  experimental  apparatus,  and  measurement  techniques.  The  third  section  provides  the  results  for 
the  mean  flow  upstream,  on  the  surfaces,  and  in  the  wake  of  the  hydrofoil,  with  emphasis  on  the  trailing-edge  region 
under  variations  of  Reynolds  number  and  geometry.  These  data  are  used  to  determine  how  the  growth  of  boundary 
layers  on  the  suction  and  pressure  sides  of  the  hydrofoil  influence  the  separated  region  at  the  trailing  edge.  The 
fourth  section  discusses  these  results  and  attempts  to  link  the  observed  phenomena  to  well-established  Reynolds 
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number  trends  in  turbulent  wall-bounded  flows  and  wakes.  The  final  section  summaries  the  effort  and  presents  the 
conclusions  drawn  from  this  work.  ,  .  : 

2.  Experimental  Setup  and  Techniques 

The  high  Reynolds  numbers  of  this  experiment  were  achieved  with  a  large  hydrofoil  tested  in  the  U.S. 

Navy’s  William  B.  Morgan  Large  Cavitation  Channel  (LCC).  The  LCC  is  a  low  turbulence  re-circulating  water 
tunnel  with  a  3.05  m  x  3.05  m  x  13  m  test  section  and  is  capable  of  steady  flows  from  0.25  to  18.3  m/s  (Etter  and 
Wilson  1992).  The  hydrofoil  spanned  the  test  section  and  was  mounted  near  its  geometric  center  (see  Figure  1). 

Gaskets  at  the  hydrofoil-wall  junction  prevented  bypass  flow,  and  tunnel  blockage  based  on  the  hydrofoil  s 
maximum  thickness  (0.171  m)  was  6%.  The  ratio  of  the  foil  chord  to  the  vertical  test  section  dimension  was  0.70. 

The  foil’s  angle  of  attack,  a,  was  measured  between  the  tunnel  axis  and  the  flat  pressure  side  of  the  hydrofoil.  For 
all  results  reported  here,  the  foil  was  mounted  at  a  =  0°,  measured  within  an  estimated  uncertainty  of +/-0.050.  At 
the  maximum  flow  speed  of  18.3  m/s  and  a  =  0°,  the  hydrofoil  generates  approximately  590  kN  (60  metric  tons)  of 
lift.  At  this  speed  and  the  facility’s  limit  on  water  temperature  (104°  F),  a  maximum  Rec  value  of  61M  was 
achieved. 

Testing  was  conducted  in  three  experimental  campaigns  (separate  testing  periods  of  approximately  six- 
week  duration)  with  the  first  campaign  measurements  limited  to  laser  Doppler  velocimetry  (LDV).  Between  the  first 
and  second  testing  campaigns,  debris  was  discovered  in  the  LCC’s  flow  management  section,  and  an  assessment  has 
been  made  of  its  impact  on  the  test  results.  Measurements  from  the  first  campaign  (with  debris)  showed  spatial 
variation  in  the  time-averaged  inflow  ffeestream  velocity  of  approximately  +/-1%,  over  the  measured  flow  range  of 
3  to  18.3  m/s.  Measurements  of  the  root-mean-square  (RMS)  of  freestream  velocity  near  the  leading  edge  of  the 
hydrofoil,  though  insufficient  to  fully  characterize  the  inflow  turbulence,  limit  its  value  to  +/-1.5%  of  the  velocity  far 
upstream.  After  debris  removal,  extensive  freestream  uniformity  and  turbulence  measurements  were  made  in  the 
empty  test  section  at  a  streamwise  location  corresponding  roughly  to  the  hydrofoil’s  leading  edge.  Over  the  full 
range  of  facility  flow  speeds,  spatial  variation  of  the  mean  freestream  velocity  was  within  +/- 0.5%,  and  the  RMS  of 
freestream  velocity  was  within  0.5%  of  the  velocity  far  upstream  (Cutbirth  and  Park  2002).  These  conditions  apply 
to  the  second  and  third  testing  campaigns  in  which  selected  trailing  edge  LDV  measurements  were  repeated  and  all 
other  data  were  acquired.  Though  the  debris  measurably  affected  inflow  uniformity  and  may  have  elevated 
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freestream  turbulence,  a  critical  comparison  of  the  LDV  measurements  repeated  between  campaigns  did  not  reveal 
flow  variation  exceeding  experimental  uncertainty. 

The  hydrofoil  test  model  was  cast  from  Ni-Al  bronze,  machined  to  the  specified  shape  (Jessup  1999)  and 
polished  to  a  RMS  surface  roughness  of  0.25  pm.  Using  flat  plate  boundary  layer  scaling,  this  roughness  represents 
k*  <  1  at  the  highest  test  speed,  and  thus  the  clean  foil  may  be  considered  hydrodynamically  smooth  (White,  1991). 
However,  for  PIV  flow  measurements  the  tunnel  was  flood-seeded  with  silver-coated  glass  spheres  of  16  pm  mean 
diameter.  The  gradual  accumulation  of  these  particles  varied  the  effective  surface  finish  of  the  hydrofoil,  and  the 
large  tunnel  drain  and  fill  times  made  frequent  cleaning  impractical.  The  diameter  of  these  particles  was  equivalent 
to  k+  ~  20  at  the  highest  Rec  tested.  As  will  be  discussed  in  Section  3.2,  these  particles  may  have  affected  boundary 
layer  transition.  Otherwise,  no  boundary  layer  tripping  device  was  employed. 

The  hydrofoil  cross  section  is  depicted  in  Fig.  2  and  provided  in  table  form  in  Appendix  A.  Both  the  figure 
and  coordinates  reflect  a  hypothetical  sharp  trailing  edge.  In  reality,  a  0.4-mm  radius  was  applied  to  the  geometry 
specified,  reducing  the  chord  by  a  negligible  amount.  Fig.  2a  also  defines  the  coordinate  frames  used  in  this  paper. 
In  the  tunnel  coordinate  frame,  the  streamwise  coordinate,  x,  is  defined  by  the  tunnel  axis,  and  vertical  coordinate,  y, 
is  taken  normal  to  x.  (The  spanwise  coordinate,  z,  completes  the  set  of  right-handed  Cartesian  coordinates).  The 
vertex  of  the  trailing  edge  bevel  angle  is  defined  as  the  coordinate  (x/C,  y/C)  =  (1,0).  In  the  surface-aligned 
coordinate  frame,  the  t  coordinate  is  defined  by  the  local  surface  tangent  with  the  h  coordinate  taken  normal  to  t. 
Finally,  some  results  are  presented  in  a  rotated  Cartesian  coordinate  frame  defined  by  h  and  /  at  93  /o  chord.  This 
frame  is  designated  the  trailing  edge  coordinate  frame. 

A  bolt-on  trailing  edge  modification  was  used  to  vary  the  suction  side  bevel  geometry.  The  two  bevel 
designs  tested  are  shown  on  Fig.  2b:  the  more  streamline  trailing  edge,  Geometry  I  (bevel  radius,  RB  =  76.2  mm  and 
apex  angle,  |3  =  44°),  and  the  more  bluff  trailing  edge,  Geometry  II  (RB  =  38.1  mm,  {5  =  56°).  Both  geometries 
generate  suction  side  boundary  layer  flow  separation  in  the  last  2%  of  the  chord.  Laser-based  metrology  on  the 
finished  Geometry-I  hydrofoil  confirmed  a  typical  tolerance  of  +/-0.2  mm  and  a  maximum  tolerance  of  +/-0.5  mm 
between  the  as-designed  and  as-built  geometries.  The  trailing  edge  of  Geometry  II,  which  was  applied  in  the  field, 
could  not  be  constructed  to  as  high  a  tolerance.  Therefore,  the  geometry  given  in  Figure  2b  and  Appendix  A  is  the 
as-built  geometry,  measured  to  within  +/-0.5  mm. 
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The  hydrofoil  experiments  were  conducted  using  well-established  water  tunnel  testing  techniques  for  fluid 
velocity  and  static  pressure  measurements.  Test  section  pressure  was  held  constant  and  sufficiently  high  to  suppress 
significant  cavitation.  Test  section  flow  speed  was  controlled  through  the  rotational  velocity  of  the  tunnel’s  axial- 
flow  impeller  and  monitored  by  a  fixed  LDV  probe.  Tunnel  water  temperature  was  monitored  but  not  controlled;  it 
increased  by  as  much  as  1.3°C/hr  during  tests  at  18.3  m/s  and  varied  between  24  and  40°C  over  the  course  of  the 
tests.  As  a  result,  water  viscosity  and  the  Rec  achieved  at  a  given  test  speed  varied  by  as  much  as  20%.  Table  1  lists 
the  nominal  Rec  values  based  on  an  average  water  temperature  of  32°C.  Measurements  of  vibration  in  a  bandwidth 
from  2  Hz  to  1  kHz  were  made  with  eight  accelerometers  (Wilcoxon  754-1)  distributed  within  the  hydrofoil.  These 
measurements  confirmed  that  the  hydrofoil  was  effectively  rigid  even  when  unsteadiness  developed  in  the  foil’s  near 
wake. 

Static  pressure  measurements  were  made  on  the  hydrofoil  using  30  taps  distributed  along  the  suction  and 
pressure  side  surfaces.  Each  tap  was  progressively  offset  by  in  the  spanwise  direction  by  40  mm,  and  all  taps  were 
contained  within  the  middle  50%  of  the  span.  Leading  edge  taps  at  x/C  <  0.03  were  0.8  mm  in  diameter  and  the 
remaining  taps  were  1.6  mm  in  diameter.  The  measurements  were  made  with  a  Rosemont  differential  pressure 
transducer  (3051  CD,  +/-250  kPa,  differential)  referenced  to  a  1.6  mm-diameter  tunnel  wall  tap  at  x/C  =  -2.38,  and 
routed  to  individual  hydrofoil  taps  with  a  rotary  sampling  valve  (Scanivalve  J-type).  The  pressure  transducer  was 
calibrated  with  a  pressure  standard  (Paroscientific  740  DigiQuartz,  124  kPa,  differential)  of  accuracy  +/-12  Pa. 
Significant  sources  of  static  pressure  measurement  uncertainty  include  (1)  uncertainty  in  the  calibration  method,  (2) 
transducer  non-linearity,  (3)  zero  bias,  and  (4)  hole  error  (see  Benedict  1984).  At  Rec=50M,  hole  error  governs,  and 
the  overall  uncertainty  is  +/-0.006  in  pressure  coefficient  units.  At  ReCf=4M  and  excluding  zero  bias,  calibration 
error  governs,  and  the  overall  uncertainty  is  +/-0.015.  However,  at  all  but  the  highest  Rec,  transducer  drift  between 
calibrations  produced  zero  bias  in  excess  of  the  other  uncertainties.  The  LCC’s  flow  deceleration  time  precluded 
zero-flow  measurements  at  a  frequency  sufficient  to  address  this  zero  bias.  Consequently,  corrections  for  zero  bias 
have  been  made  by  matching  the  pressure  coefficient  near  mid-chord  (average  of  the  three  mid-chord  taps)  on  the 
hydrofoil’s  pressure  side  to  that  on  the  Geometry  I  hydrofoil  at  Rec  =  50M.  This  method  produces  static  pressure 
distributions  that  are  consistent  with  the  steady  Bernoulli  equation  using  LDV-acquired  velocities  outside  the 
boundary  layer  on  the  pressure  side  of  the  trailing  edge  and  on  the  suction  side  near  mid-chord.  The  method  also 
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yields  pressure  coefficient  values  at  the  leading  edge  tap  consistent  with  the  free-stream  stagnation  pressure. 
Computation  of  pressure  gradients,  lift,  and  drag  are  not  affected  by  the  zero-bias  error  correction. 

Static  pressure  measurements  were  also  made  on  the  walls  and  ceiling  of  the  test  section.  Wall  taps 
extended  along  three  lines  from  -3.0  <  x/C  <3.0  at  y/C  =  —0.19,  0.04,  and  0.26,  respectively.  Single  taps  were 
located  at  mid-span  on  the  test  section  ceiling  (y/C=0.71)  at  x/C  =  0.0,  0.5,  and  1.0.  All  wall  and  ceiling  taps  were 
1.6  mm  in  diameter  and  were  routed  via  solenoid-actuated  ball  valves  to  a  second  differential  pressure  transducer 
(Rosemount  305 IP,  6350  mm-water  differential).  For  the  tunnel  wall  pressure  data,  zero  flow  measurements  were 
used  to  correct  the  transducer  zero  bias.  At  Rec  =  50M,  the  overall  uncertainty  is  governed  by  hole  error  and  is 
estimated  to  be  ±0.006  in  pressure  coefficient  units.  At  Rec  =  4M,  the  uncertainty  of  the  calibration  governs  and  the 
overall  uncertainty  rises  to  ±0.015. 

Flow  velocities  above  and  below  the  foil  were  measured  with  three  systems:  i)  a  two-component  Dantec 
laser-Doppler  velocimetry  (LDV)  at  0.25  span,  it)  a  two-component  Particle  Imaging  Velocimetry  (PIV)  system  at 
0.36  span  and  iii)  a  miniaturized  one-component  LDV  probe  (Fourguette,  et  al. ,  2001)  housed  within  the  hydrofoil 
body  at  x/C  =  0.43  and  0.35  span.  Table  2  provides  the  LDV  and  PIV  measurement  locations.  Results  from  all  three 
systems  were  normalized  by  the  flow  speed  upstream  of  the  foil,  U0,  measured  at  x/C  =  -2.38  with  a  dedicated  one- 
component  Dantec  LDV  system.  To  support  LDV  acquisition  in  all  three  testing  campaigns,  the  LCC  was  flood- 
seeded  with  Silicon  Carbide  particles  of  2  pm  mean  diameter.  To  support  PIV  acquisition  in  the  second  and  third 
testing  campaign,  the  channel  was  additionally  seeded  in  much  higher  density  (0.004  kg/m3)  with  silver-coated  glass 
spheres  of  16  pm  mean  diameter  (Potters  Industries,  SH400S33). 

The  fixed  and  traverse-mounted  LDV  systems  utilized  Dantec  components  including  signal  processors 
(BSA57N1 1),  fiber  optic  probes,  a  three-dimensional  probe  traverse,  and  control  software.  Laser  light  was  provided 
by  two  6- Watt  Argon-Ion  lasers  (Spectra  Physics  2016  and  2017).  The  external  2-component  LDV  used  a  Dantec 
optical  head  having  1 1 1-mm  beam  spacing  and  a  1600-mm  (in  air)  focal  length.  The  head  was  mounted  on  a  large 
traverse  and  provided  an  in-water  focal  volume  of  approximately  170  pm  in  diameter  and  6  mm  in  length.  The 
traverse  and  the  test  section  windows  allowed  LDV  data  collection  over  approximately  the  middle  1  m  of  the  test 
section  height.  This  LDV  system  was  calibrated  as  a  unit  using  a  rotating  disk  to  generate  the  velocity  standard,  and 
separate  calibrations  were  used  for  low  speeds  (Rec=17M  and  below)  and  high  speeds  (Rec=33M  and  50M).  The 
uncertainty  in  locating  the  disk  center  is  the  dominant  source  of  velocity  bias  error,  and  generates  uncertainties  of +/- 
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0.04  m/s  and  +/-0.08  m/s  (at  the  95%  confidence  level)  for  the  low  and  high-speed  calibrations,  respectively.  In 
comparison  with  bias  error,  random  error  is  negligible.  To  estimate  the  uncertainty  in  the  normalized  velocity  data, 
a  distinction  must  be  made  between  results  for  which  data  is  taken  and  normalized  by  a  single  LDV  system,  and  data 
taken  and  normalized  using  separate  systems.  In  the  former,  the  uncertainty  in  the  normalized  value  is  0.08  and 
0.004  in  normalized  units  at  Rec=1.4M  and  50M,  respectively.  In  the  later,  uncertainties  are  generally  higher. 

To  best  select  a  spanwise  location  for  the  LDV  measurements,  the  spanwise  uniformity  of  the  Geometry  I 
hydrofoil  flow  was  assessed.  Measurements  at  Rec  —  8M  and  50M  were  collected  on  a  coarse  grid  (50  mm  by  50 
mm)  in  a  planes  perpendicular  to  the  flow  direction  at  fixed  x-locations  far  upstream  of  the  foil  and  just  aft  of  the 
trailing  edge  (Bourgoyne,  2001).  Within  the  resolution  of  the  grid,  these  measurements  showed  the  time-averaged 
flow  to  be  two-dimensional  over  the  middle  90%  of  the  foil  span.  (Side-wall  contamination  extends  out  to  only  5  /a 
span  on  each  side.)  Motivated  by  the  higher  LDV  data  rates  near  the  test  section  windows,  this  finding  was  used  to 
justify  collection  of  the  remaining  LDV  at  z/S=0.25.  (Further  findings  concerning  spanwise  uniformity  are  presented 
in  Section  3.4.) 

Measurements  in  the  trailing  edge  region  were  taken  to  confirm  repeatability  between  campaigns.  These 
revealed  that  the  first  campaign  LDV  measurements  differed  from  those  of  the  second  two  campaigns  by  a 
consistent  factor  of  1.04.  This  difference  was  traced  to  the  calibration  of  the  upstream  LDV  system  that  monitored 
U0.  (Due  to  this  error,  the  normalized  velocity  data  published  in  Bourgoyne  (2001)  is  overstated  by  a  constant  factor 
of  1.04.)  The  error  was  systematically  corrected  using  comparisons  between  the  measurements  of  U0  from  the  fixed 
single-component  LDV  and  the  measurements  of  the  local  freestream  velocities  from  the  traverse-mounted  two- 
component  LDV.  The  correction  method  relies  on  the  assumption  of  negligible  Rec-dependence  in  the  local 
freestream  at  coordinates  y/C  >  0.26.  Following  this  correction  to  U0,  all  LDV  velocity  statistics  from  all  three 
campaigns  match  within  experimental  uncertainty. 

Planar  PIV  instantaneous  flow  field  measurements  in  the  vicinity  of  the  hydrofoil  trailing  edge  were  made 
using  a  LaVision  Flowmaster-3S  PIV/PTV  system,  including  two  digital  cameras  (1280  by  1024  pixels)  and  a  PC- 
based  data  acquisition  system.  Light  was  provided  by  two  flash-lamp-pumped  Nd-YAG  lasers  (Spectra  Physics  Pro 
250-10)  delivering  800  mJ  per  pulse  at  532  nm.  The  laser  sheet  was  masked  to  3  mm  thickness  and  relayed 
downward  through  the  top  of  the  LCC  test  section  to  illuminate  the  suction-side  trailing  edge  and  near  wake.  The 
pressure  side  of  the  hydrofoil  was  not  illuminated.  Due  to  test  section  geometry,  the  laser  sheet  could  not  be  installed 
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at  the  spanwise  location  of  the  LDV  and  was  instead  installed  at  the  nearest  practical  location  (z/S=0.36).  The  two 
cameras  were  operated  in  tandem  at  a  magnification  of  0.1  mm/pixel  to  capture  a  composite  field  of  view.  Camera 
depth  of  field  exceeded  the  light  sheet  thickness.  Raw  images  from  the  two  cameras  were  processed  individually 
using  LaVision’s  DaVis  v6.0.4  analysis  software  in  cross-correlated  mode  with  32  by  32  pixel  interrogation  areas. 
Accordingly,  each  PIV  vector  produced  is  the  result  of  particle-pair  averaging  over  a  cube  of  flow  measuring 
approximately  3  mm  on  a  side.  Images  were  acquired  at  a  rate  of  approximately  1  Hz. 

The  PIV  image  magnification  was  optimized  to  capture  in  a  single  frame  the  wake  physics  of  interest  and 
produce  typical  particle  image  diameters  of  3  pixels.  The  optimal  time  between  images  was  governed  by  the  low 
velocities  in  the  re-circulating  region  of  the  near  wake,  and  produced  a  field-averaged  particle  displacement  of 
approximately  7  pixels.  This  determined  the  minimum  suitable  interrogation  area  of  32  by  32  pixels,  typically 
capturing  10  particle  pairs.  The  raw  images  were  good.  For  example,  when  processed  with  a  single-pass  calculation 
for  each  interrogation  area,  they  produced  vector  fields  with  less  than  10%  erroneous  vectors  and  an  average  peak 
ratio  (ratio  of  the  magnitudes  of  highest  and  second-highest  cross-correlation  peaks)  of  1.8.  When  processed  with  an 
adaptive,  multi-pass  algorithm,  the  percentage  of  erroneous  vectors  dropped  to  less  than  3%  and  the  average  peak 
ratio  improved  to  2.6.  In  the  plots  presented  in  this  work,  the  adaptive,  multi-pass  algorithm  was  used,  and  vectors 
with  high  uncertainty  (Q<1.8)  were  dropped  from  the  sample  population. 

PIV  uncertainties  include  random  and  bias  error  in  the  velocity  measurements  and  bias  error  in  the  spatial 
location  of  these  velocities.  Assuming  a  characteristic  sub-pixel  resolution  of  ±0.25  pixels,  the  uncertainty  due  to 
random  error  in  the  velocity  magnitude  derived  from  particle  pairs  with  a  7-pixel  displacement  is  ±3%.  This 
estimate  is  supported  by  the  measured  self-consistency  of  the  PIV.  Specifically,  images  of  a  particle  field  were  taken 
simultaneously  with  different  cameras  and  processed  by  the  same  PIV  algorithm.  Comparison  of  the  resulting 
instantaneous  vector  fields  produced  an  L2  error  norm  of  3%  for  the  velocity  magnitudes.  An  identical  comparison 
of  the  time-average  of  500  instantaneous  fields  yielded  an  L2  error  norm  of  less  than  1%.  The  reduction  in  the  error 
norm  by  time-averaging  confirms  the  predominately  random  nature  of  the  error  and  supports  the  estimate  of  ±3% 
random  error. 

Bias  error  in  the  PIV  velocities  is  introduced  through  uncertainty  in  image  scaling.  Images  of  a  precision 
calibration  grid  mounted  to  the  submerged  hydrofoil  (in  the  absence  of  flow)  provided  this  scale  with  approximately 
±1%  accuracy.  Velocity  uncertainty  is  also  introduced  through  the  normalization  velocity,  U0i  acquired  with  the 
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upstream  reference  LDV,  discussed  earlier.  The  combined  effect  of  these  sources  of  error  is  a  characteristic 
uncertainty  of  the  normalized  mean  PIV  velocities  of  ±1.5%.  Note  that  this  estimate  is  overstated  in  regions  of  high 
velocity  and  understated  in  regions  of  low  velocity.  •  - '  *=  >  * 

Simultaneous  images  from  the  two  cameras  are  located  relative  to  one  another  within  ±0.2  pixels  using  a 
PIV  correlation  of  the  overlapping  region  of  the  images.  The  composite  velocity  field  is  then  located  relative  to  the 
foil  within  ±2  pixels  (±0.2  mm)  using  the  image  of  the  trailing  edge  surface  and  its  laser  sheet  shadow. 

3*  Results 

Experimental  data  were  collected  for  both  trailing-edge  geometries  for  zero  angle  of  attack  at  the  flow 
speeds  and  nominal  Reynolds  numbers  listed  in  Table  1.  This  table  also  includes  the  symbols  and  line  types  used  in 
the  remainder  of  this  paper  to  designate  Rec  values.  Geometry  I  and  II  results  are  reported  with  filled  and  open 
symbols,  respectively.  The  presentation  of  results  follows  in  subsections  devoted  to:  static  and  dynamic  foil 
deflections,  static  pressure  profiles,  outer  (potential)  flow  results,  and  inner  (viscous)  flow  results. 

3.1  Hydrofoil  Deflection  and  Vibration 

Though  the  foil  was  made  of  solid  metal  and  mounted  with  the  greatest  practical  rigidity,  the  static  lift  at 
Rec=50M  and  a=0.0  deflected  the  trailing  edge  8mm  at  the  midspan.  At  Rec=33M,  this  deflection  decreased  to  3.5 
mm.  Leading  edge  static  deflections  were  not  measured  but  should  approximate  those  of  the  trailing  edge  based  on 
the  near-symmetry  of  the  hydrofoil’s  cross  section,  lift  distribution,  and  mounting  scheme.  Also  based  on  the 
mounting  scheme  and  system  symmetries,  the  span  wise  dependence  of  the  static  deflection  should  approximate  that 
of  a  simply  supported  beam  under  a  uniform  load.  Deflection  of  this  kind  may  potentially  affect  the  hydrofoil  flow. 

However,  the  ratio  of  the  maximum  midspan  deflection  to  the  total  hydrofoil  span  is  less  than  0.002.  For  such  small 
deflections,  the  flow  field  is  presumed  to  locally  approximate  that  of  the  undeflected  shape. 

Measured  fluctuating  acceleration  levels  increased  rapidly  with  flow  speed  but  were  broadband  and  did  not 
suggest  problems  with  flow-excited  oscillations  of  foil-tunnel  structural  modes.  At  18.3  m/s  (the  worst  case),  the 
highest  RMS  surface-normal  vibratory  speed  recorded  by  any  of  the  eight  accelerometers  was  1.5  mm/s.  This  RMS 
velocity  is  less  than  0.01%  of  U0)  significantly  below  the  turbulence  level  of  the  water  tunnel.  In  addition,  the 
greatest  observed  change  in  vibration  level  occurred  between  flow  speeds  of  12  and  18.3  m/s  yet  the  differences 
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between  the  normalized  mean  flow  fields  at  these  speeds  are  well  within  experimental  uncertainty.  Together  these 
observations  suggest  that  foil  vibration  was  insignificant  in  these  experiments. 

3.2  Static  Pressures 

The  pressure  coefficient,  CP  -  ( P -Pa)  I  \pU2a  where  P  is  the  measured  static  pressure  and  p  is  the  water 

density,  was  measured  on  the  surface  of  the  hydrofoil  and  on  the  walls  of  the  LCC  test  section.  The  reference 
pressure  P0  was  measured  at  the  same  upstream  location  as  the  reference  velocity  U0.  Measured  CP  values  for  Rec 
between  4M  and  50M  are  shown  on  Figures  3  and  4  for  Geometries  I  and  II,  respectively.  Following  convention, 
the  vertical  axes  on  these  figures  display  -CP;  this  places  the  suction  side  data  in  the  upper  half  of  each  figure.  For 
clarity,  data  at  Rec  =  8M  and  33M  are  omitted  from  Fig.  3,  and  data  from  Rec  =  17M  and  33M  are  omitted  from  Fig. 
4.  In  both  cases,  the  omitted  data  fall  monotonically  between  the  plotted  points  for  neighboring  Rec  values.  The 
fitted  curves  are  cubic  splines  constrained  to  pass  within  the  uncertainty  range  of  each  data  point.  For  each  flow 
speed,  the  CP  value  at  x/C  =  1  was  extrapolated  from  the  suction  and  pressure  side  measurements  nearest  the  trailing 
edge.  Variations  in  the  static  pressure  distributions  between  measurements  at  different  Rees  are  small,  and  are  most 
discemable  near  the  foil’s  midchord  and  trailing  edge.  The  left  and  right  inset  panels  in  Figs.  3  and  4  show  the 
leading  and  trailing  edge  Cr  data  at  greater  scale,  with  the  location  of  the  suction  side  separation  shown  as  a  shaded 
region  on  the  x/C  axis.  The  sharp  changes  in  C^-curve  slope  near  the  trailing  edge  are  nearly  coincident  with  the 
initiation  of  the  suction  side  bevel  depicted  in  Fig.  2b  and  denoted  by  ‘Bevel’  on  the  right-side  insets  of  Figs.  3  and 
4.  Note  that  the  suction  side  tap  at  x/C= 0.958  with  the  Geometry  II  fell  within  the  faired  seam  of  the  trailing  edge 
modification  and  did  not  have  a  sharp-edged  hole.  As  a  result,  the  measurements  there  have  increased  uncertainty. 

A  comparison  of  Figs.  3  and  4  shows  that  the  CP  measurements  are  similar  but  there  is  less  variation  with 
increasing  Rec  with  Geometry  II.  As  might  be  expected  from  the  geometrical  difference,  Geometry  II  allows  the 
foil  to  maintain  a  negative  suction-side  Cp  to  greater  x/C  but  then  the  Cp  turns  positive  with  a  steeper  pressure 
gradient  as  x/C  approaches  unity.  On  both  trailing  edges,  the  pressure  gradient  just  prior  to  boundary  layer 
separation  becomes  less  adverse  with  decreasing  Rec,  a  trend  linked  to  the  decrease  in  the  -C^-peak  near  midchord 
on  the  suction  side.  For  both  edges,  the  CP  at  x/C  =  1  (the  base  pressure  coefficient)  decreases  with  decreasing  Rec, 
an  indication  that  the  near  wake  dynamics  are  changing  over  this  Rec  range.  Specifically,  the  dynamic  flow 
measurements  reveal  a  correlation  between  reduced  base  pressure  and  increased  proximity  to  the  trailing  edge  of 
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developed  pressure  side  vortices  (not  shown  here).  A  similar  relationship  was  demonstrated  by  Bearman  (1965)  for 
blunt  bodies  with  splitter  plates,  and  is  consistent  with  the  increase  in  drag  generally  associated  with  vortex 

shedding.  The  Re-dependence  of  the  base  pressure  is  weaker  for  Geometry  II. 

Figure  5  presents  CP  measured  on  the  tunnel  walls  above  (y/C=0.26)  and  below  (y/C=  -0.19)  the  foil  for 
Rec  =  4M  and  50M  for  Geometry  I,  and  at  Rec  =  50M  for  Geometry  II.  The  foil’s  mounting  system  prevented  the 
collection  of  wall  pressure  data  near  mid-chord,  so  an  interpolated  data  point  (*)  has  been  added  at  the  x/C  of  foil’s 
peak  CP.  This  value  is  the  average  of  the  measurements  on  the  foil  surface  and  water-tunnel  ceiling  at  the  given  x/C. 
With  this  point,  the  splined  C/>- curves  and  measured  velocities  conserve  vertical  momentum  in  a  control  volume 
enclosing  the  foil  and  described  in  Section  4.5.  Here,  variation  in  CP  with  Rec  and  trailing  edge  geometry  is  less  than 
experimental  uncertainty  except  very  near  the  foil. 

The  foil’s  lift  and  test  section  blockage  both  contribute  to  static  pressure  variations  on  the  test  section  walls. 
Given  the  length  of  the  foil’s  chord  compared  to  the  transverse  test  section  dimension,  the  lift  on  the  foil  produced 
by  bound  circulation  was  determined  by  the  upstream  flow  speed  and  the  confining  effects  of  the  test  section  walls. 
At  matching  free-stream  speeds,  foil  lift  would  be  lower  in  an  infinite  environment.  Interestingly,  the  measured  CP. 
on  the  ceiling  of  the  test  section  directly  above  the  foil  (-0.28,  see  Fig.  5a)  can  be  used  to  determine  the  foil’s  lift  to 
within  a  few  percent  using  of  a  two-dimensional  vortex  to  represent  the  foil’s  bound  circulation,  a  source  and  sink  to 
represent  the  foil’s  thickness,  and  the  method  of  images  to  account  for  the  tunnel  walls. 

3.3  Outer  mean  flow 

The  measurements  presented  in  this  section  document  the  global  flow  around  the  Geometry  I  foil, 
specifically  the  tunnel-confined  potential-flow.  Based  on  the  minimal  geometry-dependence  of  the  static  pressure 
coefficient,  this  data  is  expected  also  to  approximate  the  flow  over  the  Geometry  II  foil.  These  measurements  are 
presented  as  vertical  profiles  of  horizontal  mean  velocity,  U,  and  vertical  mean  velocity,  V,  normalized  by  U0.  Data 
are  presented  at  stations  upstream  of  the  foil,  at  the  trailing  edge,  and  downstream  of  the  foil.  The  vertical  extent  of 
the  measurements  (see  Table  2)  is  the  maximum  afforded  by  the  test  section  windows.  All  data  are  for  z/S=0.25. 

The  upstream  LDV-acquired  average  profiles  of  U  and  V  are  provided  in  Fig.  6a  and  6b,  respectively,  for 
Rec  =  8M,  and  50M  just  upstream  of  the  hydrofoil  leading  edge  at  x/C  =  -  0.014.  The  foil’s  leading  edge  lies  at 
(x/C,  y/C)=  (0,0.0092)  and  is  indicated  by  the  dashed  line  in  the  figure.  The  (U,  F)-profiles  at  Rec  =  17M  and  33M 
are  identical  within  data  scatter  to  the  results  at  Rec  =  8M  and  50M. 
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Normalized  mean  flow  velocities  (U,V)  above  and  below  the  foil  near  its  trailing  edge,  at  x/C  =  0.930, 
0.958,  and  1.000  are  shown  in  Fig.  7  for  Rec  =  8M,  and  50M.  Here,  a  small  difference  is  seen  in  the  data  from  the' 
two  plotted  Rec  values,  especially  for  U  on  the  suction  side  of  the  foil.  Measurements  from  Rec  =  17M  and  33M 
have  been  omitted  for  clarity.  If  shown,  these  data  would  fall  smoothly  and  monotonically  between  the  plotted  data 
in  Fig.  7.  Velocity  profiles  from  these  three  x/C  locations  have  been  included  because  they  show  the  strong  suction- ' 
side  outer  flow  deceleration  (compare  Fig.  7  a  and  e)  that  occurs  near  to  but  upstream  of  the  trailing-edge  bevel. 

Normalized  mean  flow  wake  profiles  downstream  of  the  foil  at  x/C  =  1.028,  1.070,  and  1.43,  are  shown  in 
Fig.  8  for  Rec  =  8M  and  50M.  As  for  the  near-trailing  edge  measurements  in  Fig.  6,  only  small  differences  are 
apparent  between  the  data  from  these  two  Rec’s,  and  the  data  from  Rec  =  17M  and  33M  (not  plotted)  fall  smoothly 
and  monotonically  between  the  plotted  results.  At  the  farthest  downstream  location  (Fig.  8e  and  f),  evidence  of  the 
lift  generated  by  the  foil  is  apparent  in  the  non-uniformity  of  U ,  and  the  negative  average  value  of  V.  Here,  the 
vertical  extent  of  the  measurements  reported  in  Figs.  6-8  hides  the  details  of  the  near- foil  boundary  layer  and  wake 
flows  that  have  thicknesses  of  order  0.01  C.  The  next  section  presents  measurements  of  the  mean  flow  close  to  the 
foil  surfaces. 

3.4  Inner  mean  flow 

This  subsection  covers  the  measured  boundary  layer  and  near-wake  mean  flow  close  to  the  foil’s  trailing 
edge.  In  this  flow  region,  the  LDV  results  are  complemented  with  PIV  results,  for  which  more  flow  conditions  were 
measured.  Unfortunately,  a  systematic  comparison  of  the  PIV  and  LDV  data  revealed  inconsistencies. 

The  foil  and  its  mounting  system  were  intended  to  produce  two-dimensional  flow  over  as  much  of  the  foil’s 
span  as  possible.  However,  unrepeatable  variation  in  the  mean  flow  was  measured  with  the  PIV  system,  particularly 
at  Rec  =  8M.  The  observed  behavior  is  attributed  primarily  to  uncontrolled  variation  with  time  and  with  spanwise 
location  in  the  downstream  location  of  suction  side  boundary  layer  transition.  Evidence  to  this  effect  is  found  in  the 
reduction  or  disappearance  of  these  measurement  variations  at  the  highest  and  lowest  RecS.  At  the  highest  Recs, 
suction  side  transition  is  believed  to  occur  uniformly  close  to  the  foil’s  leading  edge.  At  the  lower  Recs,  suction  side 
transition  occurs  uniformly  near  70%  chord  where  the  adverse  pressure  gradient  first  becomes  steep  (see  Fig.  25). 
Unfortunately,  foil  surface  imperfections  or  the  accumulation  of  PIV  seed  may  have  been  controlling  factors  for 
suction  side  boundary  layer  transition  in  the  mid-range  Rec’s  of  this  study  (4M  <  Rec  <  17M).  Though  the  effect  is 
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small  in  the  attached  boundary  layers,  the  resulting  variation  in  suction  side  boundary  layer  separation  appears  to 
amplify  the  effect.  Figure  9  shows  the  worst-case  variations  observed  in  the  near-wake  mean  velocity  profiles 
(Geomtery  I  at  Rec  =  8M).  The  thicker  wake  profile  in  Fig.  9  is  from  the  LDV  measurements  at  25%  span.  The 


thinner  wake  profile  is  from  the  PIV  measurements  at  36%  span.  The  thinner  PIV  profile  was  never  measured  by 
the  LDV  system  at  this  Rec.  In  summary,  some  variation  of  the  inner  mean-flow  profiles  at  the  trailing  edge  was 
observed  at  36%-span  in  the  mid-range  Rees  of  this  study.  To  address  this,  the  data  presented  have  been  restricted  to 
the  LDV  mean  flow  measurements  made  at  the  25%-span  location,  and  those  PIV  fields  (measured  at  36%  span) 
that  agree  with  these  LDV  measurements. 

A  measure  of  the  viscous  flow  on  the  Geometry  I  foil’s  suction  side  was  provided  by  measurements  of  Ut 
(—  the  surface-tangent  mean  velocity)  made  at  x/C  =  0.43  with  a  miniaturized  one-component  LDV  mounted  inside 
the  foil  (Fourguette  et  al.  2001).  These  measurements  are  plotted  in  Fig.  10  for  Rec  =  8M,  17  M,  33M  and  SOM 
along  with  laminar  and  turbulent  boundary  layer  profile  fits.  Following  surface-aligned  coordinates,  h  is  the  local 
surface  normal  coordinate  and  Ute  is  the  surface-tangent  mean  velocity  at  the  edge  of  the  boundary  layer  ( h  —  8) 
where  8  is  selected  for  a  best  fit  for  the  boundary  layer  profile.  Here,  8  and  other  fitting  parameters,  C/  —  skin 
friction  coefficient  and  n  =  Coles’  wake  parameter,  are  given  in  Table  3  for  the  Blasius  profile  fit  at  Rec  =  8M,  and 
the  Coles’  profile  fit  (see  White  1991), 


u: 


=  —  ln/z+  +  B  + 

K 


2U  sin2  (nhf  /25) 
k  sin2(jif/2) 


(1) 


for  the  other  Recs.  Here  the  usual  definitions  apply:  U]  =  (U J U le)^2/C^ ,  /z+  -  hUte^Cf  /2jv,  k  -  0.41,  B  - 

5.5,  v  =  0.775 xlO-6  m2/s  (at  the  average  water  temperature),  8  =  <5„  and  the /-parameter  in  (1)  was  set  to  unity. 
Neither  a  fully-laminar  nor  a  fully- turbulent  boundary  layer  profile  could  be  successfully  fit  to  the  Rec  —  17M  data 
which  suggests  a  transitional  boundary  layer  at  this  location  and  Rec-  These  measurements  indicate  that  suction  side 
boundary  layer  transition:  i)  occurred  downstream  of  x/C  ~  0.43  at  Rec  =  8M,  ii)  was  underway  but  incomplete  x/C 
=  0.43  at  Rec  -  17M,  and  iii)  was  complete  upstream  of  x/C  =  0.43  at  Rec  =  33M.  The  load-bearing  structure  of  the 
foil  and  the  LCC  precluded  velocimetry  at  other  mid-chord  locations. 

The  next  nearest  downstream  location  at  which  boundary  layer  profiles  could  be  measured  was  x/C  =  0.930. 
Figure  1 1  shows  the  boundary  layer  Ut  profiles  for  the  suction  side  (a)  and  pressure  side  (b)  at  this  location  for  the 
Geometry  I  hydrofoil.  To  present  the  data  in  surface-aligned  coordinates,  the  suction  side  data  measured  in  vertical 
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profiles  at  x/C=0.930  and  0.953  have  been  interpolated  to  approximate  the  data  along  the  surface  normal  at 
x/C=0.930.  For  Rec=1.4M  and  4M,  data  was  not  available  at  x/C= 0.953,  so  the  surface  normal  profile  is  estimated 
from  the  vertical  profile  using  y  ~  A .  A  zero  velocity  point  has  been  added  to  each  profile  at  the  foil  surface  where 

neither  PIV  nor  LDV  measurements  were  possible.  All  the  profiles  represent  turbulent  boundary  layers  and  their 
corresponding  fitting  parameters  are  provided  in  Tables  4  and  5  for  the  suction  and  pressure  sides  of  the  foil, 
respectively.  Here,  the /-parameter  in  (1)  and  the  added  degree  of  freedom  provided  by  5*5e  were  used  to 

enhance  the  smoothness  of  the  profile  fits.  Also,  a  vertical  coordinate  shift  within  the  uncertainty  in  the  LDV  s 
spatial  reference  was  allowed  to  achieve  B— 5.0  simultaneously  with  a  best  fit  to  the  linear-log  region  of  the 
boundary  layer  profiles.  Note  that  the  approximation  y  *»h  used  for  the  low  Rec  data  is  expected  to  have  a  small 

effect  on  the  boundary  layer  fits.  As  a  test  case,  the  same  approximation  was  applied  to  the  Recf=8M  data,  producing 
a  2%  overstatement  of  JT  and  negligible  effect  on  the  other  fitting  parameters. 

An  overview  of  the  mean  flow  fields  near  the  Geometry  I  and  II  foil  derived  from  the  PIV  measurements  is 
provided  in  Figs.  12-14.  The  first  two  of  these  figures  show  multiple  velocity  profiles  of  U  (Fig.  12a  &  13a)  and  V 
(Fig.  12b  &  13b)  above  and  downstream  of  the  foil,  which  is  depicted  on  the  left  side  of  each  panel.  The  gray 
vertical  lines  indicate  the  positions  at  which  the  profiles  were  measured.  One  or  two  horizontal  arrows  connects  a 
given  flow  profile  with  its  measurement  location.  A  relative  scale  for  the  flow  speeds  is  given  at  the  lower  left  of 
each  panel.  The  more  nearly  horizontal  curves  on  each  figure  indicate  the  suction-side  and  pressure-side  boundary 
layer  thicknesses  above  and  below  the  foil,  and  the  boundaries  of  the  foil  s  wake  beyond  x/C  of  unity.  Excluding 
Rec=4M,  three  general  trends  are  observed:  1)  the  wake  thins  as  Rec  increases,  2)  the  wake  is  thicker  with  Geometry 
II,  and  3)  the  reverse  flow  region  is  longer  with  Geometry  II.  Interestingly,  Geometry  II  at  Rec  =  4M  is  the  lone 
exception  to  all  three  of  these  trends.  Figure  14  provides  normalized  U  and  V  values  for  Geometry  I  (a)  and 
Geometry  II  (b)  at  the  boundary-layer  and  wake  edges  shown  on  Figs.  12  and  13.  These  velocities  are  used  to 
normalize  several  of  the  following  plots. 

The  PIV  measurements  near  the  trailing  edge  and  Tecplot®  9.0  software  were  used  to  compute  apparent 
stagnation  points  and  separating  streamlines.  These  are  shown  on  Fig.  15  where  arrows  indicate  flow  direction  and 
the  interior  of  closed  streamlines  are  shaded.  The  apparent  stagnation  point(s)  are  marked  with  black  dots,  and 
dashed  lines  represent  streamlines  that  must  exist  but  were  either  not  resolved  or  not  recorded  within  the  PIV 
measurement  suite.  These  streamline  patterns  are  subject  to  uncertainty  from  the  finite  PIV  accuracy  and  from 
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weak,  unintended  (and  unmeasured)  spanwise  flow,  particularly  in  the  recirculation  region.  In  addition,  streamlines 
under  the  foil  could  not  be  determined  because  shadowing  prevented  PIV  measurements  there. 

Both  geometry  and  Re-dependence  is  apparent  in  these  streamline  patterns,  with  Geometry  I  at  Rec  =  1.4M 
(Fig.  15a)  providing  the  most  unique  case.  Here,  there  is  no  indication  of  stagnation  or  separation  at  the  trailing 
edge  and  the  streamtube  that  emerges  just  below  the  foil  appears  to  turn  upwards  without  separating  from  the  foil’s 
trailing  edge  (a  complete  violation  of  the  Kutta  condition),  and  then  flow  backwards  along  the  foil  surface  until  it 
separates  on  the  suction  side.  This  places  an  apparent  stagnation  point  just  aft  of  the  foil’s  trailing  edge  on  the  lower 
extent  of  a  detached  recirculation  zone.  However,  the  volume  flux  in  this  contorted  streamtube  lies  at  the 
confidence  limits  of  the  PIV  measurements.  It  may  be  best  described  as  speculative  because  weak  spanwise  flow 
gradients  -  normal  to  the  PIV  plane  -  might  be  holding  this  streamtube  open  and  any  spanwise  gradient  changes 
might  close  this  streamtube  in  an  orderly  fashion.  Nevertheless,  the  Rec  =  1.4M  case  remains  distinguished  because 
the  streamline  flow-angle  farther  aft  of  the  foil,  which  should  be  much  less  susceptible  to  unintended  spanwise  flow 
gradients,  differs  from  that  shown  for  the  seven  other  panels  of  Fig.  15. 

The  remaining  flow  cases  (Fig.  15b-15h)  place  a  stagnation  point  at  the  trailing  edge  that  captures  the 
suction-side  separating  streamline  and  launches  the  dividing  streamline  into  the  wake.  The  suction-side  separating 
streamline  encloses  a  clockwise-rotating  recirculation  region.  Secondary  counterclockwise-rotating  recirculation 
regions  were  found  with  Geometry  II  (Fig.  15b  and  15d).  For  Geometry  I,  such  counter-clockwise  recirculating 
regions  and  associated  stagnation  points  are  absent  or  unresolved. 

The  location  of  suction  side  boundary  layer  separation  was  estimated  from  the  data  used  to  create  Figs.  12 
and  13  by  two  methods,  and  the  results  appear  on  Figure  16.  In  the  first  method,  the  loci  of  points  at  which  Ut  -  0 
are  plotted  in  surface-aligned  coordinates.  Excluding  near-wall  data  contaminated  by  laser  glare  and  downstream 
points  for  which  the  surface  normal  has  extended  into  the  wake,  this  loci  of  points  closely  approximates  a  straight 
line.  Taking  advantage  of  this  linearity,  the  loci  points  for  x/C< 0.998  and  over  the  ranges  0.00 15<  hfC  <  0.0050  and 
0.0015<  hIC  <  0.0080  for  Geometries  I  and  II,  respectively,  were  linearly  extrapolated  to  the  surface.  This  surface 

point  was  taken  as  the  zero-velocity  separation  point,  X^ep .  The  second  method  is  identical  to  the  first,  except  that 
the  plotted  points  are  the  loci  of  50%  probability  of  forward  tangential  velocity  and  the  extrapolation  produces  X^ep. 
These  two  methods  would  produce  identical  results  in  the  case  of  a  symmetric  probability  distribution  function  for 
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the  tangential  velocity.  However,  the  results  on  Fig.  16b  show  that  Xsep  is  consistently  larger  than  x^ep  by  a  small 

amount.  On  Fig.  16a  two  results  are  plotted  for  Geometry  I  at  Rec  =  8M  corresponding  to  the  two  trailing  edge  flow 
fields  measured  at  this  condition  (see  Fig.  9).  The  solid  dark  symbol  is  for  the  thicker  wake  that  is  consistent  with 
the  LDV  measurements. 

In  Fig.  16a,  the  two  trailing  edge  geometries  show  similar  trends  of  suction  side  separation  with  Rec  for 
Rec>3M.  However,  as  Rec  decreases  below  3M,  suction  side  separation  on  Geometry  I  moves  rapidly  downstream, 
while  separation  on  Geometry  II  is  roughly  constant.  This  behavior  is  attributed  to  differences  in  boundary  layer 
separation  on  the  pressure  side.  At  the  lowest  Rees,  the  pressure  side  boundary  layer  is  laminar  or  transitional  when 
it  reaches  the  trailing  edge,  while  the  suction  side  boundary  layer  is  turbulent.  As  these  boundary  layers  evolve 
downstream  and  encounter  an  increasingly  adverse  pressure  gradient  at  the  trailing  edge,  the  laminar  pressure  side 
boundary  layer  thickens  more  quickly  than  the  turbulent  suction  side  layer,  and  is  believed  to  detach  prior  to  the 
trailing  edge.  On  Geometry  I,  the  suction  side  boundary  layer  is  then  able  to  advance  further  downstream.  The 
Geometry  II  edge  sees  a  similar  competition  between  its  laminar  pressure  side  and  turbulent  suction  side  boundary 
layers,  but  the  steep  trailing  edge  bevel  prevents  the  suction  side  from  advancing  further  downstream.  As  a  result  the 
pressure  side  boundary  layer  remains  attached.  Early  pressure  side  separation  on  Geometry  I  provides  an 
explanation  not  only  for  Fig.  16,  but  also  for  the  unique  appearance  of  panel  (a)  among  the  other  panels  of  Fig.  15. 

The  evolution  of  the  suction  side  boundary  layer  across  separation  for  both  trailing  edges,  is  provided  in 
Fig.  17  at  x/C  =  0.930,  0.979,  0.998,  and  1.028  and  Rec  =  1.4M,  4M,  8M  and  50M.  Here,  the  trailing  edge 
coordinate  frame  (a  fixed,  rotated  Cartesian  coordinate  system  defined  by  h  and  t  at  93%  chord)  is  used  to  determine 
the  Ut  profiles  at  each  x/C  value.  This  coordinate  system,  which  is  rotated  f  =  14.12°  from  the  x-y  system,  was 
chosen  because  it  approximates  a  streamline  coordinate  system  for  all  flow  conditions  in  the  region 
0.930<x/C<1.028.  The  Ut  profiles  are  also  vertically  aligned  to  match  at  hei  the  location  of  the  boundary-layer  edge 
or  near- wake  edge  in  the  93%-chord  Cartesian  frame.  The  data  show  that  the  thinnest  boundary  layers  occur  at  ReC 
=  4M,  and  that  the  outer  part  of  the  boundary  layer  flow  evolves  more  rapidly  at  the  lower  two  Recs  than  at  the 
higher  two  Recs.  In  fact,  at  the  higher  two  Recs  the  outer  part  of  the  boundary  layers  (when  normalized  by  the  local 
external  velocity)  appears  to  be  insensitive  to  the  presence  or  absence  of  the  foil  surface.  Under  this  velocity 
normalization  and  coordinate  frame,  the  outer  portion  of  the  attached  boundary  layers  crosses  separation  with  little 
change,  and  the  Re-dependence  of  the  attached  boundary  layers  sets  the  Re-dependence  of  the  initial  wake. 
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This  geometry  dependence  is  explored  more  fully  on  Fig.  18  which  shows  trailing  edge  profiles  of  Ut  at 
Rec  =  1.4M,  4M,  17M,  and  SOM  from  the  suction  side  at  x/C=1.000  and  the  pressure  side  at  x/C  =  1.002.  Here  the 
suction  side  results  are  again  plotted  in  the  trailing  edge  coordinate  frame,  but  the  profiles  are  now  vertically  aligned 
by  the  6 -coordinate  of  the  foil  surface  at  suction  side  separation,  hsep.  The  pressure  side  profiles  simply  present  U 
plotted  in  the  x-y  coordinates  frame.  (On  the  pressure  side  U  =  U,  and  ysep  =  0,  so  Fig.  18  actually  presents  the 
pressure  and  suction  side  results  on  an  equal  footing.)  The  normalizing  velocities  can  be  derived  from  the  data  on 
Fig.  14.  The  three  higher  Recs  show  good  profile  agreement  on  the  pressure  side  and  on  the  suction  side  for  h  -  hsep 
>  0.  Thus,  the  main  geometric  effect  of  the  trailing  edge  modification  on  the  wake  is  communicated  through  the 
suction  side  separation  point.  The  Rec  =  1.4M  case  does  not  fit  this  geometrical  trend  because  the  pressure  side 
boundary  layer  separates  prior  to  the  trailing  edge  on  Geometry  I  and  at  the  trailing  edge  on  Geometry  II.  Therefore 
a  shift  accounting  only  for  suction  side  separation  does  collapse  the  geometry  dependence. 

The  flow  profiles  that  appear  in  Fig.  18  initiate  the  foil’s  wake.  Figures  19,  20  and  21  show  the  subsequent 
flow  evolution  in  x-y  coordinates  at  x/C  =  1.009,  1.028,  and  1.047,  respectively,  for  Rec  -  1.4M,  4M,  17M,  and 
50M.  Here,  the  flow  has  been  scaled  as  a  plane  wake;  the  horizontal  axes  are  AU/AU^  where  A U-U  -  Ufl,  and 
the  vertical  axes  are  (y  -  ycti)lym  where  yc,r  denotes  the  average  of  the  vertical  coordinates  corresponding  to  A  U  — 
AC/max/2,  and  2ym  is  the  full  wake  width  where  AC/=  AC/max/2.  On  each  frame  of  Figs.  19-21  a  Gaussian  profile, 
matched  at  the  half-width  points,  is  plotted  to  indicate  the  extent  of  wake  convergence  to  a  standard  symmetrical 
form.  On  Fig.  19  (x/C  =  1.009),  every  profile  shows  some  skewness  as  a  remnant  of  the  foil  boundary  layers’ 
asymmetry,  with  Rec  =  4M  being  the  most  symmetrical.  On  Fig.  20  (x/C  =  1.028),  the  profile  skewness  is  reduced 
in  every  case  and  is  nearly  absent  at  Rec  =  4M.  On  Fig.  21,  the  various  wakes  have  relaxed  to  a  self-similar 
Gaussian  profile  and  subsequent  wake  evolution  involves  little  or  no  further  profile  shape  changes. 

Given  the  similarities  shown  on  Figs.  19-21,  wake  flow  differences  can  be  readily  assessed  by  examining 
how  yclr,  At/max,  and  y\a  evolve  downstream  of  the  foil.  The  first  of  these  is  shown  of  Fig.  22  for  where  (yclr-y*)!C 
is  plotted  vs.  x/C.  Here,  the  value  of  y*  is  the  vertical  increment  obtained  by  extrapolating  the  foil  thickness  at 
suction  side  separation  downstream  to  the  trailing  edge  along  a  line  parallel  to  the  /-axis  of  the  93  /o-chord  rotated 
Cartesian  system:  y*  =  (l.O-  x^/cJtanjS-  This  vertical  shift  nearly  collapses  the  wake  center  results  for  the  two 

trailing  edge  shapes  at  each  Rec  except  for  Rec  =  1.4M.  This  further  confirms  that  the  main  impact  of  the  differing 
trailing  edges  on  the  foil’s  mean  flow  can  be  accounted  for  with  a  vertical  shift  related  to  suction  side  separation.  As 
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mentioned  earlier,  the  Rec  =  1.4M  case  involves  early  pressure  side  separation  and  a  shift  based  on  suction  side 
separation  is  insufficient  to  collapse  the  geometry  dependence. 

The  downstream  evolution  of  the  wake  velocity  deficit,  At/max,  is  shown  in  Fig.  23  for  both  trailing  edges  at 
Rec  =  1.4M,  4M,  17M,  and  50M.  As  before,  the  normalizing  velocity  U“  is  extracted  from  the  data  on  Fig.  14. 
The  curves  are  downward  trending  as  expected,  and  the  strongest  Reynolds  number  variation  reveals  a  pairing  of  the 
two  higher  Rec  and  the  two  lower  Rec-  Finally,  the  downstream  evolution  of  the  wake  half  width,  y\a,  is  shown  in 
Fig.  24  for  both  trailing  edges  and  the  same  four  Recs.  The  curves  are  upward  trending  as  expected,  and  there  is 
substantial  Reynolds  number  variation.  In  particular,  at  Rec  =  4M  with  Geometry  II  the  wake  leaves  the  trailing 
edge  as  the  thinnest  but  by  10%  of  chord  downstream  it  has  become  the  thickest,  and  this  growth  rate  is  not  matched 
nor  exceeded  by  the  flow  at  any  other  Reynolds  number. 

4.  Discussion 

The  data  presented  above  indicate  that  the  flow  near  the  trailing  edge  and  in  the  near  wake  changes  with 
both  Reynolds  number  and  trailing  edge  geometry.  These  changes  can  be  related  to  the  state  of  the  suction  and 
pressure  side  boundary  layers  upstream  of  the  trailing  edge  separation.  Here,  a  combination  of  laminar  and  turbulent 
boundary-layer  integral  calculations  were  used  to  infer  transition  location,  shape  factor  H  =  8*/d,  and  momentum 
thickness  0  for  the  two  boundary  layers  where  direct  measurements  of  these  flow  parameters  were  not  possible. 
Here  the  usual  definitions  of  5 *  and  9  apply: 


4.1  Boundary  Layer  Transition  on  the  Hydrofoil 

First,  the  averaged  pressure  distribution  (average  of  all  Recs  and  geometries)  on  the  hydrofoil  and 
Thwaites’  method  (see  White  1991)  were  used  to  compute  the  growth  of  the  laminar  boundary  layer  from  the 
leading  edge  stagnation  point  to  the  region  of  boundary  layer  transition.  The  one-step  method  of  Michel  (1952,  see 
White  1991)  was  then  used  to  predict  the  location  of  boundary  layer  transition.  The  results  of  these  calculations  are 
depicted  on  Fig.  25  against  the  speed-averaged  -Cp  curves.  The  calculated  location  of  transition  on  both  sides  of  the 
foil  moves  upstream  with  increasing  Rec-  For  Rec  below  4M,  the  flow  on  the  suction  side  is  computed  to  remain 
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laminar  until  the  steep  adverse  pressure  gradient  beginning  at  x/C  ~  0.7.  These  calculations  are  confirmed,  where 
data  is  available,  by  the  LDV-measured  velocity  profiles  from  x/C  =  0.43  (Fig.  10).  These  profiles  suggest  that 
transition  moved  upstream  of  x/C  =  0.43  between  Rec  of  8M  and  17M. 

4.2  Boundary  Layers  Approaching  TE  Region 

The  boundary  layer  parameters  obtained  from  the  Thwaites’  method  at  the  Michel-predicted  transition 
point  were  employed  as  initial  values  for  a  turbulent  boundary-layer  integral  calculation  (White  1991).  The 
momentum  thickness  was  assumed  continuous  across  the  calculated  transition  point.  Of  course,  actual  boundary 
layer  transition  occurs  over  a  finite  region  so  the  computed  results  obtained  here  are  approximate  in  this  regard.  The 
turbulent  boundary  layer  calculations  were  not  used  for  the  pressure  side  boundary  layer  below  Rec  =  3M  because 
the  PIV  results  at  x/C  =  1  revealed  a  laminar  or  transitional  boundary  layer  leaving  the  trailing  edge  on  the  pressure 
side.  At  these  conditions,  a  transition  location  was  not  calculated,  and  Thwaites  method  was  used  exclusively. 

A  sample  of  the  outcome  of  these  calculations  is  shown  on  Fig.  26  where  the  predicted  and  measured 
suction-  and  pressure-side  boundary-layer  momentum  thicknesses,  9,  at  x/C  =  0.93  are  plotted  vs.  Rec  for  Geometry 
I.  Considering  the  approximate  nature  of  the  analysis,  the  agreement  between  the  calculations  and  measurements  is 
good.  For  the  suction  side,  the  location  of  transition  is  roughly  fixed  at  Rec  below  4M  (see  Fig.  28).  As  a  result,  as 
Rec  increases  from  the  lowest  values  to  Rec=4M,  9  decreases  due  to  boundary  layer  thinning  with  increasing  Rec. 
However,  as  Rec  increases  above  4M,  the  transition  point  begins  to  move  upstream  and  increases  the  turbulent 
development  length  of  the  boundary  layer.  Thus,  9  increases  with  increasing  Rec  for  Rec>  4M.  As  a  result,  Rec=4M 
produces  a  minima  in  the  suction  side  boundary  layer  thickness.  The  Re-dependence  of  the  pressure  side  boundary 
layer  thickness  is  due  to  similar  effects,  but  without  a  fixed  transition  point.  In  this  case,  Rec=4M  corresponds 
roughly  to  a  maxima  in  the  pressure  side  boundary  layer  thickness. 

The  success  of  these  boundary  layer  calculations  has  implications  concerning  the  Re  and  geometry 
dependence  of  the  hydrofoil  flow.  Since  the  evolution  of  9  with  increasing  x/C  is  well  predicted  by  a  single  CP 
distribution,  the  Re  and  geometry  dependence  of  the  potential  flow  must  have  a  weak  influence  on  boundary  layer 
development.  This  means  that  the  upstream  feedback  from  separation  and  the  wake  does  not  appreciably  affect  the 
state  of  the  attached  boundary  layers.  Instead,  the  character  of  these  attached  boundary  layers  is  denved  primarily 
from  the  Re-dependence  of  transition  and  its  effect  on  the  laminar  and  turbulent  boundary  layer  development.  This 
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claim  is  made  more  significant  by  the  finding  of  Fig.  17,  which  shows  the  Re-dependence  of  the  outer  portions  of 
the  attached  boundary  layers  crossing  separation  with  little  change.  This  behavior,  as  well  as  other  aspects  of 
separation  and  the  near  wake  will  be  discussed  in  the  following  section. 

4.3  Separating  Boundary  Layers  and  Initial  Wake 

This  section  discusses  the  more  complex  Re  and  geometry  dependence  of  the  flow  nearest  to  the  trailing 
edge.  In  this  region,  the  boundary  layers  at  x/C  =  0.93  further  evolve  to  separate  and  form  the  near  wake.  On  the 
suction  side  of  the  foil  aft  of  x/C  =  0.93,  a  steep  adverse  pressure  gradient  leads  to  boundary  layer  separation  and  the 
boundary-layer  integral  methods  lose  validity.  However,  suction  side  boundary  layer  integral  quantities  are  still 
expected  to  exhibit  some  similarities.  As  discussed  in  Simpson  (1989),  boundary  layers  near  separation  exhibit  a 
relationship  between  a  simple  function  of  the  shape  factor,  h  =  (H-l)/H,  and  the  ratio  5*ide .  Figure  27  presents 
these  two  quantities  for  the  separating  suction  side  boundary  layers  on  both  trailing  edges  along  with  correlations 
developed  from  backward  facing  steps,  h  =  l.5(8*/8e)  and  power-law  profiles,  h  =  (2-8*ISe)  1  (see  Simpson  1989). 
The  intersection  of  these  two  correlations  indicates  the  point  of  intermittent  boundary  layer  detachment  on  a  flat- 
plate.  Here,  H  and  8*/8e  are  computed  from  Ut  in  surface-aligned  coordinates,  and  the  data  plotted  range  from 
X7  /C  -  0.004  <x/C  <  x\ep/C  + 0.011  for  Geometry  I  and  x\ep/C  -  0.003  <x/C<  x]ep/C  +  0.003  for  Geometry 

II.  The  open  symbols  on  Fig.  27a  and  the  filled  symbols  on  Fig.  27b  for  each  flow  condition  correspond  to  the 
boundary  layer  state  at  50%  probability  of  forward  flow. 

For  surfaces  of  low  curvature  and  pressure  gradient,  the  paths  of  the  boundary  layers  in  the  vicinity  of 
separation  are  roughly  linear  and  collapse  on  the  line  of  slope  1.5  indicated  in  Figure  27.  Several  of  the  data  sets 
presented  here  have  linear  regions  with  this  slope,  but  the  overall  collapse  to  a  single  line  is  imperfect,  as  might  be 
expected  for  the  curved-surface  foil  boundary-layer  flows  of  this  study.  However,  Fig.  27  does  reveal  a  relationship 
between  boundary  layer  integral  quantities  at  separation,  and  demonstrates  both  Rec  and  geometry  dependence  in 
the  flow.  The  Rec-dependence  is  apparent  in  the  variation  between  paths  at  different  values  of  Rec.  In  particular  the 
Rec  =  4M  data  generally  lie  above  the  plotted  points  for  the  other  Recs.  Trailing-edge  geometry  dependence  is 
apparent  in  the  grouping  of  the  data,  with  the  Geometry  I  results  (Fig.  27a)  lying  above  the  Geometry  II  results  (Fig. 
27b).  This  places  the  higher  surface  curvature  results  from  Geometry  II  further  from  the  flat  plate  results.  The 
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extreme  case  high  surface  curvamre  at  separation  is  achieved  by  the  foil’s  separating  pressure  side  boundary  layer 
where  H~  1.3  (as  expected,  see  Schlichting  1979). 

The  data  presented  in  Fig.  17  also  show  that  when  a  turbulent  boundary  layer  separates  from  a  rapidly 
curving  surface,  there  is  a  lag  in  the  communication  of  the  near  wall  changes  to  the  outer  flow,  a  finding  is  similar  to 
that  of  Morris  and  Foss  (2001).  In  particular,  above  Rec  =  17M,  the  momentum-containing  portions  of  the  separated 
boundary  layers  at  the  trailing  edge,  x/C  =  1.0  are  essentially  the  same  as  the  outer  portions  of  the  attached  boundary 
layers  at  x/C  ~  0.98  (when  normalized  by  the  local  freestream).  In  addition,  the  trailing-edge  geometry  dependence 
of  the  suction-side  boundary  layer  flow  as  it  leaves  the  foil  appears  to  be  determined  by  the  location  of  the  suction 
side  separation  point  for  ReC  >  4M.  In  the  trailing  edge  coordinate  frame  (93%-chord  Cartesian  coordinates),  the 
suction-side  Ut  profiles  for  Qi  -  hsep)IC  >  0  are  essentially  identical  between  the  two  geometries,  at  least  at  the  higher 
three  Reynolds  numbers,  and  the  primary  flow  field  variations  are  confined  to  (h  -  hsep)/C<  0  (Fig.  18). 

On  the  pressure  side  of  the  foil  beyond  x/C  =  0.93,  the  boundary  layer  development  story  is  much  simpler. 
Here,  the  pressure  side  boundary  layers  encounter  a  mildly-adverse  pressure  gradient  as  x/C  increases  toward  unity. 
In  the  absence  of  this  pressure  gradient,  separation  is  expected  at  the  trailing  edge.  Pressure-side  boundary  layer 
profile  measurements  made  for  0.93  <  x/C  <  1  are  limited  to  Geometry  I  at  Rec  of  8M,  17M,  33M,  and  50M. 
Results  for  Geometry  II  at  these  Rec  values  were  expected  to  be  nearly  identical  between  geometries  because  of 
their  similar  static  pressure  distributions  (see  Figs.  3  and  4).  All  pressure-side  velocity  profiles  are  very  similar  to 
flat  plate  results,  the  displacement  and  momentum  thickness  vary  according  to  established  correlations,  and  skin 
friction  coefficients  inferred  from  the  pressure  side  velocity  profiles  agree  well  with  the  Schultz-Grunow  correlation 
(see  Schlichting  1979).  Although  pressure  side  boundary  layer  profiles  were  not  recorded  for  Rec  below  8M,  the 
initial  wake  velocity  profiles  from  x/C  =  1.002  (Fig  18)  retain  the  character  of  the  pressure  side  boundary  layer. 
With  the  exception  of  Geometry  I  at  Rec  =1.4M,  these  nearest- wake  profiles  suggest  pressure-side  boundary  layer 
separation  occurs  at  the  trailing  edge. 


4.4  Development  of  the  Wake 

As  the  near-wake  of  the  foil  converts  to  a  far-wake,  its  characteristics  should  trend  toward  those  of  an  ideal 
two-dimensional  wake  with  a  constant  momentum  thickness.  Here  the  momentum  thickness,  6W,  of  the  foil  s  wake 
is  defined  by 
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where  <5/“  and  <5“  are  defined  by  the  vertical  locations  where  U  =  £/”.  The  Geometry  I  LDV  profiles  on  Fig.  8 
were  used  to  calculate  an  average  momentum  thickness  for  the  far-wake  at  Rec  =  8M  to  50M,  dw/C  -0.0039. 
Due  to  the  difficulty  in  defining  5?  and  <5“ ,  the  uncertainty  in  6W  is  high  and  a  Re-specific  calculation  is  not 
justified.  The  wake-center  velocity  deficit,  AUmix  ,  and  half-width  ym,  were  then  estimated  from  the  wake  scaling 

laws  in  Sreenivasan  and  Narasimha  (1982):  AUma/Ug  *<1.63| ~d/(x  —  and  “  0.30[(x  —  x0 >  where 

x„  is  a  hypothetical  origin  of  the  wake.  These  scaling  relationships  are  plotted  on  Figs.  23a  and  24a,  where  x0  is 
selected  for  a  best  fit  and  is  1.01  for  the  scaling  of  A  Umax  and  0.91  for  the  scaling  of  ym.  The  same  values  for  *0 

and  dw  (we  lack  far  wake  data  for  Geometry  II)  are  used  to  produce  the  curves  on  Figs.  23b  and  24b.  Although  the 
match  to  the  measured  data  is  better  for  the  wake  velocity  deficit  A  ,  the  wake  scaling  laws  appear  reasonably 


successful  the  higher  Rec  flows  beyond  x/C  -  1.03. 

Overall,  the  foil’s  wake  flow  changes  little  from  Rec  of  8M  to  50M,  but  the  findings  at  Rec  <  4M  show 
much  more  variability.  For  example,  the  wake  profile  at  Rec  =  4M  for  both  trailing  edges  relaxes  to  near-perfect 
Gaussian  form  closest  to  the  trailing  edge  when  compared  to  the  wake  profiles  at  the  other  Recs  (see  Figs.  19  and 
20).  These  observations  suggest  there  is  augmented  turbulent  transport  occurring  in  the  foil’s  near  wake  at  Rec  ” 
4M.  This  contention  is  supported  by  the  observation  of  vortex  shedding  fluctuations  in  the  foil’s  near  wake  at  this 
Rec.  Additional  discussion  of  vortex  shedding  is  deferred  to  a  future  paper  on  this  topic. 

4.5  Hydrofoil  Lift  and  Drag 

The  overall  suite  of  measurements  is  sufficient  to  allow  the  calculation  of  the  foil’s  lift  and  drag  by  two 
nearly  independent  methods  for  Geometry  I.  With  Geometry  II,  only  one  of  these  methods  was  viable.  These 
calculations  were  undertaken  to  establish  foil  performance  and  to  validate  the  internal  consistency  of  the 
measurements  made  with  Geometry  I.  Here,  lift  and  drag  are  reported  as  the  two-dimensional  lift  and  drag 
coefficients  (CL  and  CD)  with  the  foil  chord  Cused  as  the  length  scale  in  both  cases. 
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The  first  method  involved  integrating  the  foil’s  surface  stresses  to  determine  foil  lift  and  drag.  The  input 
data  here  were  the  splined  pressure  distributions  and  the  skin-friction  stress  calculated  from  the  boundary-layer 
integral  analysis.  From  this  method,  Ci  for  Geometry  I  was  determined  to  vary  approximately  linearly  with  Rec 
from  0.52  to  0.55  ±0.01  for  Rec  between  8M  and  50M,  and  CL  for  Geometry  II  was  determined  to  be  0.52  ±0.01  and 
constant  within  the  given  uncertainty  for  the  same  range  of  Rec.  Here,  Co  was  found  to  be  0.006  ±0.001  and 
constant  with  Rec  within  the  given  uncertainty  for  both  trailing  edge  geometries.  In  all  cases,  approximately  60%  of 
the  total  drag  originated  from  skin  friction,  and  40%  from  pressure.  The  uncertainty  in  lift  is  derived  from  that  of 
the  static  pressure  measurements,  discussed  earlier.  For  the  drag  calculations,  the  dominant  source  of  uncertainty  is 
the  skin  friction  values  predicted  using  the  boundary  layer  integral  method.  This  uncertainty  is  estimated  at  +/-20% 
by  comparing  the  predicted  C/ values  to  those  derived  from  the  linear-log  region  of  the  measured  velocity  profiles  at 
x/C  =  0.93. 

The  second  method  employed  a  two-dimensional  rectangular  control  volume  enclosing  the  hydrofoil  flow 
and  defined  by  the  comer  coordinates  (x/C,  y/C)  -  (-1.25,  -0.19),  (-1.25,  0.26),  (1.43,  0,26),  and  (1.43,  -0.19).  A 
uniform  inflow  velocity  and  the  LDV  far-wake  survey  data  (Fig.  8)  were  used  as  velocity  boundary  conditions  on 
the  upstream  and  downstream  control  surfaces.  Wall  static  pressure  (Fig.  5)  confirmed  that  the  upstream  pressure  is 
uniform,  but  showed  that  the  pressure  in  the  far  wake  location  was  not.  Therefore,  the  steady  incompressible 
Bernoulli  equation  was  used  to  provide  the  pressure  on  the  downstream  control  surface,  with  the  exception  of  the 
region  of  the  wake  where  measured  velocity  fluctuations  are  substantial  (--0.07  <  y/C  <  0.05).  Here,  linear 
interpolation  between  the  Bernoulli-derived  pressures  was  used.  The  splined  test  section  wall  static  pressures, 
including  the  interpolated  value  on  the  suction  side  at  x/C  =  0.7  (Fig.  5),  were  used  to  provide  the  pressure  boundary 
conditions  on  the  upper  and  lower  control  surfaces.  Again,  the  Bernoulli  equation  was  utilized  to  determine  flow 
velocity  magnitudes  along  this  surface.  The  assumed  pure  horizontal  inflow  direction,  and  the  measured  flow 
directions  at  the  leading  edge,  at  the  trailing  edge,  and  in  the  far  wake  were  then  linearly  interpolated  to  provide  a 
continuous  estimate  of  the  flow  direction  along  the  upper  and  lower  control  surfaces.  Lift  and  drag  contributions 
from  turbulent  fluxes  were  not  included.  These  control  volume  calculations  were  undertaken  for  Rec  =  8M,  17M, 
33M,  and  50M.  The  imbalance  in  mass  conservation  over  this  range  of  Rec,  for  the  control  volume  described,  was 
no  greater  than  ±0.1%  of  the  mass  flux  at  the  inflow  surface.  This  result  is  sensitive  to  the  interpolation  scheme 
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used  to  provide  the  flow  direction  on  the  upper  and  lower  surfaces,  as  well  as  the  calibration  uncertainty  in  the  LDV 
measurements. 

From  this  second  method,  CL  for  Geometry  I  was  found  to  be  0.55  ±0.05,  with  Rec  variation  falling  within 
this  larger  uncertainty.  The  greatest  source  of  this  uncertainty  is  the  interpolated  value  of  the  wall  static  pressure  at 
x/C  =  0.7  used  on  the  upper  control  surface.  Here,  the  foil’s  CD  was  found  to  be  0.009  ±0.003,  with  Rec  variation 
again  falling  within  the  uncertainty.  For  the  drag  calculations,  the  greatest  source  of  uncertainty  is  the  +/- 1% 
calibration  uncertainty  in  the  LDV,  exacerbated  by  the  use  of  separate  LDV  systems  for  the  inflow  and  far  wake 
data  acquisition.  The  upper  and  lower  control  surface  approximations  have  minimal  effect  on  the  drag  calculations. 
Most  significantly  though,  the  two  methods  for  calculating  lift  and  drag  coefficients  for  the  Geometry  I  foil, 
produced  results  that  agree  within  uncertainty. 


5.  Summary  and  Conclusions 

The  mean  flow  around  a  large,  two-dimensional  hydrofoil  has  been  experimentally  investigated  with 
special  attention  given  to  the  separated  flow  near  its  beveled  trailing  edge.  Two  trailing  edge  shapes  were 
investigated  that  are  of  interest  for  practical  applications.  The  data  reported  here  are  unique  because  of  the  high 
Reynolds  numbers  attained.  This  research  effort  leads  to  four  conclusions. 

First  of  all,  over  the  range  of  chord-based  Reynolds  number  from  1.4  to  50  million,  and  for  the  geometries 
tested,  the  Reynolds  number  and  geometric  dependencies  of  the  flow  field  are  mild  but  clearly  measurable.  The 
majority  of  the  observed  Reynolds  number  dependence  comes  for  the  development  of  the  attached  suction  and 
pressure  side  boundary  layers  where  the  location  of  transition  moves  upstream  with  increasing  Reynolds  number. 
This  dependence  is  characterized  by  a  minima  (maxima)  near  Rec=4M  in  the  variation  with  Rec  of  the  suction 
(pressure)  side  boundary  layer  momentum  thickness.  This  distinguishes  Rec=4M  as  the  test  condition  with  the 
greatest  symmetry  between  the  suction  and  pressure  side  boundary  layers.  The  characteristics  of  these  boundary 
layers  set  the  initial  conditions  for  the  wake  development.  This  is  particularly  true  in  the  outer  part  of  the  near  wake 
where  boundary  layer  properties  persist  beyond  separation. 

Second,  much  of  the  trailing-edge  geometry  dependence  of  the  mean  flow  is  a  mere  displacement  of 
suction  side  boundary  layer  flow  over  the  blunter  trailing  edge  having  the  larger  bevel  angle  (Geometry  II).  Here, 


26 


Hydrofoil  flow  at  high  Reynolds  number 


Bourgoyne,  Hamel,  Ceccio,  and  Dowling 


boundary  layer  separation  is  moved  aft  compared  to  that  on  the  thinner  trailing  edge  (Geometry  I),  but  the  thicker 
trailing  edge  leads  to  a  thicker  near  wake.  In  addition,  when  combined  with  prior  flat  plate  measurements,  the 
observed  trends  in  turbulent  boundary  layer  separation  location  show  that  the  shape-factor  at  separation  is  lower  on 
more  highly  curved  surfaces. 

Third,  when  the  separated  boundary  layers  pass  the  trailing  edge  and  merge  to  form  the  foil  s  wake,  they 
relax  quickly  with  downstream  distance  to  conform  to  classical  wake  scaling  laws.  However,  the  manner  in  which 
the  wake  develops  depends  on  its  initial  thickness  and  symmetry.  The  most  symmetric  wake  profile  was  found  at 
Rec  =  4M  for  both  trailing  edges.  These  flows  are  distinguished  by  the  thin  suction  side  boundary  layers  and  thick 
pressure  side  boundary  layers  near  the  trailing  edge  that  create  the  most  symmetrical  near  wakes  found  in  these 
studies.  Interestingly,  the  most  symmetrical  near  wake  (Rec  =  4M,  Geometry  II)  thickens  the  most  rapidly. 
Fluctuating  flow  field  results  (not  presented  here)  suggest  near  wake  vortex  shedding  occurs  to  some  degree  for  both 
trailing  edges  at  Rec  =  4M,  and  that  this  phenomena  drives  the  mean-flow  development. 

Thus,  the  final  conclusion  can  be  stated:  the  mechanism  by  which  the  initial  wake  profile  shape,  thickness, 
and  symmetry  determine  the  subsequent  wake  evolution  is  likely  to  be  found  in  the  flow’s  dynamic  components. 
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0.0412060217 
0.0462824935 
0.0517772976 
0.8577047328 
0.8648780043 
0.0709098767 
0.0782885221 
•.88S9854S07 
0.0942474131 
0.1030002385 
0.1122479579 
■. 12199276 97 
■.1322349504 
0.1429727260 
8.1542022959 
0.1651178250 
0.1711113631 
1.1907728279 
0.2838900958 
0.2174418110 
0.2314333826 
1.2458250331 
0.2606839936 
1.2757484845 
1.2912350113 
1.3078385683 
1.3231327186 
0.3394897818 
1.3560886467 
1.3728757176 
1.3898442446 
0.4069549886 
0.4241759211 
0.4414752449 
0.4588206743 
0.4761808915 
0.493S216873 
8.518*137185 
0.5288255140 
8.5451267971 
0.5620*82184 
0.S7888144S1 
0.5954792520 
0.6118556661 
0.6279868861 
0.6438473105 
0.6594176710 
0.6746771121 
0.6896071540 
1.7041910228 
1.7184136688 
8.7322616796 
1.7457234202 
■.7587889535 
1.7714499499 
1.7836198347 
8.79SS335498 
8. *064475 446 
0.8179399178 
*.82851048*2 
1.8386597277 
1.8483901650 
8.1577050638 
8.8666090179 
0.8751075763 
8.8832872964 
8.8909156404 
1.1982407983 
1.90S1917471 
8.9117780359 
1. 9180098884 
8.4238976421 
1.9294525061 


y/C 

0.0091713492 
0.00918S5203 
0.0092184614 
0.0092850871 
0.0893826812 
0.8095115988 
0.8096744745 
0.8098724687 
0.81 81 871 728 
0.8183803244 
8.01 86936883 
0.0110487415 
•.8114474330 
0.B118913S74 
8.8123821377 
0.8129213261 
8.81351*3142 
0.8141506854 
0.8148434464 
8.8155897666 
8.8163905847 
0.8172466620 
8.0181S8574S 
8. *191267018 
0.82*1512061 
0.8212320202 
8.8223688442 
8.8235611334 
8.8248080*82 
8.8261886153 
■.8274614036 
8.8288648385 
8.8383169977 
8.8318157290 
8.8333585663 
8.8349427605 
8.8365652797 
8.8382228165 
8.8399117863 
8.8416283261 
8.8433683112 
8.0451273650 
8.8469008610 
0.0486B3934O 
8.8504715020 
0.0522582799 
8.0548387845 
1.0558873652 
0.0575582264 
8.0592854434 
1.0609829920 
0.0626447781 
0.0642646656 
0.0658365138 
0.8673542058 
8.0688116892 
•.1702838183 
1.0715223925 
1.8727642800 
8.8739230544 
0.1749938442 
8.8759717*21 
0.876IS24297 
1.8776317487 
8.1783861361 
1.8788724697 
8.8793288870 
8.0796710753 
1.8798995185 
1.18801 08  823 

8. 8880012524 

0.1798746166 
*.8796373042 
1.0792885187 
0.0788276881 
0.0782552358 
0.0775721470 
0.8767800779 
0.8758813116 
■  .8748787238 
8.8737757196 
8.8725761309 
8.8712841887 
0.0699042348 
8.06*4408389 
8.0668987565 
8.06528121*5 
8.0635912686 
*.0618347489 
0.0600003341 
8.0Si066S527 
0.8560642546 
8. 0S 40345124 
0.0520017612 
0.0499831874 
0.0479915782 
0.8460372591 
0.0441287896 
8.0422729583 
•  .*404751783 
8. 03*739*312 
8.0370701211 
0.0354684980 
8.8339365905 
0.0324752719 
0.0310847357 


x/C 

0.9294525061 

0.9396887961 
0.9442334873 
0.9485715938 
0.9526350643 
0.9S64357S01 
0.9599854319 
0.9632957666 
0.9663703603 
0.9692446364 
0.97190S1613 
0.9743730401 
0.9766569294 
0.9787679756 
0.9887163281 
0.9825117916 
0.9841638361 
0.9856815693 
0.1870737361 
0.9883487196 
0.9895145175 
1. 1185 7*7718 
0. 1915487360 
0.9924313106 
0.9932338285 
8.9939600560 
8.9952130267 
0.9962327884 
0.9970567651 
0.9977176088 
0.998243S476 
0.9986588264 
0.9989840400 
0.9991184357 


y/C 

0.8310847357 
8.8296891000 
8.8215161951 
8.8273350107 
0.8262288281 
8.8251767732 
8.8241965263 
8.023235526* 
0.0222581740 
0.8212499624 
0.0202072319 
0.0191385694 
0.8180537323 
0.8169628359 
0.8151757411 
0.81 41021581 
0.013751789a 
8.8127287364 
8.8117425347 
8.8107969129 
8.8098951616 
0.0090419434 
0. 0082370483 
0.0074820817 
8.0067769713 
0.8061215404 
0.0055146117 
>.  0044404877 
•.M3S396253 
0.0827948329 
0.OO2184S494 
0.0016921275 
0.0012986822 
0.0009876963 
0.0008584542 


Geometry  II 


X/C 

0.9294525861 

8. 970OO00O0O 
0.9713500000 
0.9721490000 
0.9729480000 
0.9737480000 
1.1745470000 
0.9753460000 
0.9761450000 
0.9769440000 
0.9777430000 
0.9785420000 
0.9793410000 
0.9801480000 
0.9814140000 
0.9827900000 
0.9837490000 
0.9847080000 
0.986282O8OO 
0.987308BOOO 
0.9880210000 
0.9890090000 
0.9897930000 
0.9906340000 
0.9913670000 
0.9920040000 
O.992511OOO0 
0.9936220000 
0.9940620000 
0.9944950000 
O.9955620OOO 
O.9966290OOO 
0.9977198000 
O.9988090OOO 


y/C 


0.0310047357 
0.0296891000 
0.0212235000 
0.0208920000 
0.8206919000 
0.0204896000 
■.0202841000 
0.02 08(773 000 
0.8198667000 
0.8196527000 
0.0194349000 
0.8192132000 
0.8189872000 
0. 0187566000 
8.8185210000 
0.8181320000 
8.8176730000 
0.0173065000 
8.8168890000 
8.0160670000 
0.0154230000 
8.0149330000 
8.8142198000 
8.8136390008 
8.8129518800 
0.812211iflaa 
0.8114850000 
0.0108540000 
0.0094390000 

8.80*2530000 

0.0067331000 

0.8051300000 

0.8034623000 

0.0017979000 


Appendix  A.  Hydrofoil  geometry  in  chord-normalized  coordinates  with  the  vertex  of  the 
trailing  edge  apex  angle  defined  as  (x/C,  y/C)  =  (1,0) 


Table  1.  Upstream  flow  velocities  and  nominal  Rec  with  symbol  and  line-types  listed.  In 
the  figures  of  this  paper,  filled  and  open  symbols  correspond  to  trailing  edge  Geometry  I 
and  II,  respectively. 

Table  2.  PIV  and  LDV  velocity  survey  locations. 

Table  3.  Parameters  of  boundary  layers  near  mid-chord  (x/C=0.43)  on  the  suction  side  of 
the  hydrofoil  with  trailing  edge  Geometry  I. 

Table  4.  Parameters  of  boundary  layers  approaching  the  trailing  edge  region  (x/C=0.930) 
on  the  suction  side  of  the  hydrofoil  with  trailing  edge  Geometry  I. 

Table  5.  Parameters  of  boundary  layers  approaching  the  trailing  edge  region  (x/C=0.930) 
on  the  pressure  side  of  the  hydrofoil  with  trailing  edge  Geometry  I. 


Table  1 


U0  Rec 

m/s  millions  Symbol  /  Line-type 


0.5 

1.4 

- f| - 

1.5 

4 

- ± - 

3.0 

8 

- - 

6.0 

17 

- ► - 

12.0 

33 

. < . 

18.3 

50 

- 0. - 

^  in  to  s 


Table  2 


_ Measurement  Station _ 

1.  Inflow  Plane 

2.  Leading  Edge 

3.  Suction  Side  Boundary  Layer  near  Mid-chord 
Boundary  Layer  approaching  Trailing  Edge 

PIV  Field  at  Trailing  Edge 
Near  Wake  Plane 
Far  Wake 


Extent  of  Measurements 


Streamwise, 

x/C 

Vertical, 

y/C 

Spanwise, 

z/S 

-2.3 

-0.20<y/C<0.28 

0<z/S<0.5 

-0.014 

-0.20<y/C<0.28 

0.25 

origin  at  x/C=0.43 

0<h/C<0.01 

0.35 

0.93 

-0.20<y/C<0.28 

0.25 

0.96<x/C<1.1 1 

-0.025<y/C<0.030 

0.36 

1.028 

-0.20<y/C<0.28 

0<z/S<0.5 

1.42 

-0.20<y/C<0.28 

0.25 

Table  3 


Rec 

[millions] 

nominal 

U0  [m/s] 

+/-  0.01 

6/C  U,e/U0 

+/- 0.0001  +/-0.01 

cf 

+/- 0.0002 

n 

+/-0.03 

8 

3.0 

0.0007  1.25 

0.0004 

— 

17 

6.0 

0.0022  1.26 

0.0042 

0.00 

33 

12.0 

0.0025  1.27 

0.0030 

0.11 

50 

18.3 

0.0030  1.27 

0.0022 

0.75 

Table  4 


Rec 

[millions] 

u0  [m/s] 

S/C 

5e/6 

nominal 

+/-  0.01 

+/-  0.0002 

+/-  0.01 

1.4 

0.5 

0.0103 

1.13 

4 

1.5 

0.0072 

1.28 

8 

3.0 

0.0103 

1.18 

17 

6.0 

0.0113 

1.11 

33 

12.0 

0.0113 

1.11 

50 

18.3 

0.0113 

1.08 

ute/u0 

cf 

n 

f 

+/-0.01 

+/-  0.0002 

+/-  0.03 

1.12 

0.0037 

0.51 

1.18 

1.12 

0.0027 

0.95 

1.10 

1.10 

0.0020 

1.50 

1.06 

1.10 

0.0017 

1.70 

1.06 

1.10 

0.0015 

1.71 

1.06 

1.10 

0.0015 

1.60 

1.06 

Table  5 


Rec 

[millions] 

nominal 

U0  [m/s] 

+/-  0.01 

S/C 

+/-  0.0002 

6e/6 

+/-  0.01 

8 

3.0 

0.0140 

1.00 

17 

6.0 

0.0140 

1.00 

33 

12.0 

0.0140 

1.00 

50 

18.3 

0.0140 

1.00 

ute/u0 

+/-  0.01 

Cf 

+/-  0.0002 

n 

+/-  0.03 

f 

0.96 

0.0025 

0.62 

1.15 

0.96 

0.0023 

0.60 

1.15 

0.96 

0.0020 

0.60 

1.15 

0.96 

0.0019 

0.58 

1.16 

Figure  Captions 

Figure  1.  The  test  section  geometry  and  hydrofoil  mounting  location,  shown  in  an 
elevation  view  looking  spanwise  and  two  section  views  looking  downstream.  All 
dimensions  are  in  meters  unless  otherwise  noted,  and  the  test  section  dimensions  are 
referenced  to  the  inside  surface.  Optical  access  at  the  leading  and  trailing  edges  required 
omission  of  the  wall  faring  fillet  radii  on  the  wall  at  z/S=0  in  the  ranges  x/C<0.19  and 
x/C>0.87.  This  configuration  is  depicted  in  Section- A.  The  ends  of  the  remaining  length 
of  wall  faring  were  smoothly  tapered. 

Figure  2.  (a)The  cross  section  of  the  Geometry  I  hydrofoil  with  chord  and  max  thickness 
indicated.  The  chord  length  given  is  the  idealized  value,  measured  to  the  vertex  of  the 
trailing  edge  apex  angle.  This  point  defines  (x/C,  y/C)  =  (1,  0)  in  the  tunnel  coordinate 
system,  for  which  the  x-axis  is  taken  as  the  streamwise  axis  of  the  tunnel.  Also  depicted 
is  the  hydrofoil  coordinate  system,  for  which  the  t-axis  is  taken  as  the  local  surface 
tangent.  (b)A  detail  of  the  Geometry  I  and  II  trailing  edges.  The  0.4-mm  radius,  applied 
to  the  trailing  edge  tip,  is  not  depicted.  The  dashed  line  indicates  the  direction  of  the 
surface  tangent  at  x/C=0.930. 

Figure  3.  The  static  pressure  coefficient,  Cp,  on  the  hydrofoil  surface  and  on  the  vertical 
centerline  of  the  test  section  side  walls  (y/C=0.04)  for  Geometry  I  at  Rec  of4M,  17M, 
and  50M.  The  stagnation  pressure  (Cp  =1.0)  was  measured  at  x/C=0  but  is  not  shown. 
The  insets  present  at  greater  scale  the  data  near  the  leading  and  trailing  edges.  Within  the 
trailing  edge  inset,  the  x/C  at  which  the  suction  side  bevel  initiates  is  labeled,  and  the 
region  of  mean  suction-side  flow  separation  is  indicated  by  a  gray  bar  on  the  x/C  axis. 
Symbols  represent  the  measured  static  pressure.  Lines  represent  a  cubic  spline  of  the 
measured  values.  The  straight  line  near  zero  Cp  connects  the  wall  data  at  the  beginning 
and  end  of  the  hydrofoil  influence. 

Figure  4.  Static  pressure  coefficient,  Cp,  for  the  hydrofoil  with  the  Geometry  II  trailing 
edge,  in  a  format  similar  to  that  of  Figure  3.  Data  at  Rec  =  17M  are  not  available  on  the 
test  section  walls  and  are  replaced  with  data  at  Rec  —  8M.  In  the  trailing  edge  inset  panel, 
the  data  from  Geometry  I  at  Rec=50M  is  also  shown  for  comparison.  Also  indicated  in 
the  trailing  edge  inset  panel  is  the  streamwise  extent  of  the  smooth  seam  between  the 
hydrofoil  and  trailing  edge  modification.  The  static  pressure  tap  in  this  region  does  not 
have  a  square  edge  and  therefore  has  higher  uncertainty. 

Figure  5.  Static  pressure  coefficient,  Cp,  on  the  test  section  walls  and  ceiling,  for  the  both 
trailing  edges  at  Rec  of  4M  and  50M,  in  a  format  similar  to  Figure  3.  (a)  Data  from  the 
test  section  ceiling,  and  (b)  data  from  the  sidewalls  on  the  suction  side  (y/C— +0.26)  and 
pressure  side  (y/C=  -0.19)  of  the  foil.  Symbols  plot  measured  values  and  the  solid  line 
splines  the  average  of  the  plotted  points.  Hydrofoil  mounting  hardware  prevented  data 
collection  on  the  sidewalls  near  midchord,  and  the  suction  side  point  at  x/C=0.65 
(denoted  by  **’)  is  the  average  of  the  measured  value  on  the  suction  side  foil  surface  and 
the  test  section  ceiling  at  Rec=50M.  The  straight  line  near  zero  Cp  is  taken  from  Figure 
3,  and  connects  the  centerline  wall  data  (y/C=0.04)  at  the  beginning  and  end  of  the 
hydrofoil  influence. 


Figure  6.  The  flow  near  the  leading  edge  (x/C  =  -0.014)  of  the  Geometry  I  hydrofoil, 
presented  as  the  (a)streamwise  and  (b)vertical  components  of  the  time-averaged  velocity, 
normalized  by  the  velocity  far  upstream,  U0.  The  dashed  line  indicates  the  location  of  the 
vertical  tangent  of  the  leading  edge  (y/C=0.009).  Data  is  shown  for  Rec=8M  and  50M. 
Data  for  Rec=17M  and  33M  are  omitted  for  graphical  clarity,  but  collapse  with  the  data 

shown. 

Figure  7  The  boundary  layer  flow  on  Geometry  I,  presented  in  the  same  format  as  Figure 
6,  and  measured  at  (a,b)  xA>0.930,  (b,c)  x/C= 0.958,  and  (e,f)  x/C=1.000.  The  location  of 
the  hydrofoil  surface  is  indicated  by  the  gray  bar.  A  zero  vertical  velocity  at  the  hydrofoil 
surface  is  assumed  and  plotted  for  clarity.  Data  are  shown  for  Rec=8M  and  50M.  Data 
for  Rec=17M  and  33M  are  omitted  for  graphical  clarity,  but  trend  monotonically  in  U0 
with  the  data  shown. 

Figure  8  The  wake  flow  of  Geometry  I,  presented  in  the  same  format  as  Figure  6,  and 
measured  at  (a,b)  x/C=1.028,  (b,c)  x/C=1.070,  and  (e,f)  x/C= 1.430.  Data  are  shown  for 
Rec=8M  and  50M.  Data  for  Rec=17M  and  33Mare  omitted  for  graphical  clarity,  but 
collapse  with  the  data  shown. 

Figure  9.  The  near  wake  flow  of  Geometry  I  at  Rec=8M,  presented  in  the  same  format  as 
Figure  6,  and  measured  at  x/C=1.009.  Symbols  are  measured  values:  A,  LDV  at 

z/S=0.25;  and - A - ,  PIV  at  z/S=0.36.  Though  pressure  side  measurements  agree 

within  experimental  uncertainty,  suction  side  measurements  differ.  The  disagreement  is 
attributed  to  spanwise  variation  in  the  location  of  suction  side  boundary  layer  transition  at 
this  Rec-  The  measurements  are  consistent  with  transition  and  separation  occurring 
further  upstream  at  z/S=0.25  than  at  z/S=0.36. 

Figure  10.  Component  of  the  time-averaged  velocity  tangent  to  the  surface  at  x/C=0.43 
on  the  suction  side  of  the  Geometry  I  hydrofoil.  Symbols  are  measured  values  and  solid 
lines  are  fits  through  the  data,  using  a  laminar  boundary  layer  profile  for  Rec-8M,  and  a 
turbulent  boundary  layer  profile  (Equ.  1)  for  Rec=17M,  33M,  and  50M.  At  Rec-17M, 
the  turbulent  profile  with  Coles  parameter,  11=0,  fits  the  data  below  h/C= 0.006.  The  data 
above  M>0.006  fit  neither  a  laminar  nor  turbulent  profile,  but  are  traced  with  a  solid 
gray  line  to  improve  graphical  clarity.  Table  3  provides  the  parameters  used  in  the  curve 
fits.  All  fits  use  k=0.41,  B=5.5,  and  v=0.775  xl0'6m2/s. 

Figure  11.  Component  of  the  time-averaged  velocity  tangent  to  the  surface  at  x/C=0.930 
on  the  (a)  suction  and  (b)  pressure  side  surfaces  of  the  Geometry  I  hydrofoil.  Symbols  are 
measured  values.  Lines  are  curve  fits  using  a  turbulent  boundary  layer  profile  (as  in 
Figure  10)  with  the  fitting  parameters  provided  in  Tables  4  and  5.  For  graphical  clarity, 
symbols  are  shown  for  only  two  Rec.  The  data  scatter  at  these  Rec  is  representative  of 
that  found  at  the  other  Rec- 


Figure  12.  The  separating  boundary  layers  and  near  wake  with  Geometry  I,  presented  as 
the  (a)streamwise  and  (b)vertical  components  of  the  time-averaged  velocity,  normalized 
by  the  velocity  far  upstream,  U0.  The  trailing  edge  geometry  is  depicted  on  the  left  side  of 
the  frame  and  flow  is  from  left  to  right.  Vertical  gray  lines  are  shown  at  the  x/C 
coordinate  of  the  measurements  and  provide  the  axes  upon  which  the  velocity  values  are 
plotted  as  black  lines.  An  arrow  extends  from  each  axis  line  to  its  associated  velocity 
profile.  The  velocity  profiles  are  plotted  within  the  range  5eps  and  the  field 

values  of  5 ”/C  and  5ps  /C are  shown  as  thin  black  lines.  The  scale  used  to  extract 
velocity  values  is  provided  in  the  bottom-left  comer  of  the  frame.  Data  at  Rec=l  .4M  was 
not  available  at  the  downstream  station. 

Figure  13.  The  separating  boundary  layers  and  near  wake  with  Geometry  II,  presented  in 
the  format  of  Figure  12.  Data  at  Rec=50M  was  not  available  at  the  three  upstream 
stations. 

Figure  14.  Streamwise  and  vertical  components  of  velocity  measured  along  the  field 
values  of  6" /C  and  8ps/C  shown  in  Figures  12  and  13.  Panels  are  for  the  (a)  Geometry  I 
and  (b)  Geometry  II.  Streamwise  (vertical)  values  are  read  from  the  axes  on  the  left 
(right)  edge  of  the  panels.  Line  type  indicates  Rec. 

Figure  15.  Stagnation  and  separation  points  and  the  associated  streamlines,  derived  from 
the  time-averaged  velocity  fields  of  the  (a-d)  Geometry  I  and  (e-h)  Geometry  II. 
(a,e)Rec=1.4M.  (b,f)Rec=4M.  (c,g)Rec=17M.  (d,h)Rec=50M.  Separation  points  are  from 
Figure  16a.  Stagnation  points  are  those  coordinates  at  which  (U,V)=(0,0).  Note  that  the 
spanwise  velocity  was  not  measured,  but  is  assumed  negligible  based  on  two- 
dimensionality.  Closed  streamlines  are  filled  with  gray.  Dashed  lines  indicate  sections  of 
streamlines  which  must  exist  but  are  not  resolved.  The  dash-dot  line  in  the  upper  right 
comer  of  each  panel  runs  parallel  to  the  surface  tangent  at  x/C=0.930.  Note  that  for  all 
conditions  and  geometry,  this  line  approximates  the  direction  of  the  streamlines  for 
0.93<x/C<1.03.  No  data  are  available  for  x/C<0.993  at  Rec=50M  for  Geometry  II. 

Figure  16.  The  location  of  separation  for  both  trailing  edges  at  varying  Rec,  based  on 
PIV  and  LDV  velocity  measurements.  Separation  was  located  using  two  methods:  by 
inflection  of  the  profile  of  the  time-averaged  velocity  tangent  to  the  surface,  £/<,  and  by 
probability  of  forward  flow,  y,  equal  to  50%,  where  ‘forward’  is  defined  by  the  direction 
of  the  surface  tangent.  Provided  are  (a)  results  from  the  first  method  and  (b)  the 
difference  of  the  two  methods.  •,  Geometry  I;  O,  Geometry  II;  ©,  Separation  at 
Rec=8M  for  Geometry  I  with  delayed  suction  side  boundary  layer  transition,  discussed  in 
Section  3.4. 

Figure  17.  The  time-averaged  velocity  profiles  across  separation,  at  varying  Rec  for  both 
trailing  edges.  The  profiles  are  presented  in  a  rotated  Cartesian  coordinate  frame  aligned 
to  the  surface  tangent  at  x/C— 0.930.  To  best  collapse  the  profiles,  the  spatial  coordinate 
has  been  shifted  by  the  boundary  layer  thickness,  <5  " ,  measured  in  the  rotated  frame  and 

designated  he.  Presented  are  measurements  at  (a)Rec=1.4M,  (b)Rec=4M,  (c)Rec=T7M, 
and  (d)Rec=50M.  Closed  symbols  are  for  Geometry  I,  while  open  symbols  are  for 
Geometry  II  at  varying  x/C:  ■,  x/C=0.930;  A,  x/C=0.979;  T,x/C=0.998;  #, 
x/C=1.028.  No  data  are  available  for  x/C=0.979  at  Rec=50M  for  Geometry  II. 


Figure  18.  The  time-averaged  velocity  profiles  for  both  trailing  edges  at  varying  Rec- 
Suction  side  data  are  from  x/C=l.  Due  to  the  shadow  of  the  foil,  pressure  side  data  was 
not  available  at  this  coordinate,  and  so  data  at  x/C= 1.002  are  presented.  The  suction 
(pressure)  side  data  are  presented  in  a  rotated  Cartesian  coordinate  frame,  aligned  to  the 
suction  (pressure)  side  surface  tangent  at  x/C=0.930.  To  best  collapse  the  geometry 
dependence  of  the  flow,  the  spatial  coordinate  of  the  suction  side  data  is  shifted  by  the 
thickness  of  the  trailing  edge  at  the  separation  point,  measured  in  the  rotated  frame  and 
designated  hsep.  No  shift  is  applied  to  the  pressure  side  data.  Presented  are  measurements 
at  (a)Rec=1.4M,  (b)Rec=4M,  (c)Rec=17M,  and  (d)Rec=50M.  Closed  symbols  are  for 
Geometry  I,  while  open  symbols  are  for  Geometry  II.  Data  at  Rec=50M  for  Geometry  I 
is  incomplete  due  to  shadowing,  but  by  all  indications  collapses  with  the  data  from 
Geometry  II. 

Figure  19.  Time-averaged  streamwise  velocity  profiles  at  x/C=1.0094,  for  the  both 
trailing  edges.  Panel  show  results  at  (a)  Rec=1.4M,  (b)Rec=4M ,  (c)Rec=17M,  and 
(d)Rec=50M.  Symbols  give  measured  values  and  the  dashed  line  gives  the  Gaussian  fit 
matched  at  the  half-width  points.  Closed  symbols  are  for  Geometry  I,  while  open 
symbols  are  for  Geometry  II. 

Figure  20.  Time-averaged  streamwise  velocity  profiles  at  x/C= 1.0281,  for  both  trailing 
edges,  in  the  format  of  Figure  19. 

Figure  21.  Time-averaged  streamwise  velocity  profiles  at  x/C= 1.0469,  for  both  trailing 
edges  at  Rec=1.4M,  4M  ,  17M,  and  50M.  Symbols  give  measured  values  and  the  solid 
line  gives  a  Gaussian  fit.  Closed  symbols  are  for  Geometry  I,  while  open  symbols  are  for 
Geometry  II 

Figure  22.  The  wake  profile  center  coordinate,  yc(/C,  used  in  Figures  19-21  and  plotted 
against  x/C  for  varying  Rec  for  both  trailing  edges.  For  all  Rec  shown  except  1.4M,  a 
collapse  of  the  geometry-dependence  of  this  parameter  is  achieved  by  shifting  yctr  by  an 
amount,  y*,  derived  from  the  hydrofoil  thickness  at  the  separation  point.  The  lack  of 
collapse  at  Rec=1.4M  is  due  to  the  variation  with  geometry  of  pressure  side  separation. 
Closed  symbols  are  for  Geometry  I,  while  open  symbols  are  for  Geometry  II. 

Figure  23.  The  wake  velocity  deficit  used  in  Figures  19-21,  plotted  against  x/C  for 
varying  Rec  and  for  (a)  Geometry  I  and  (b)  Geometry  II.  The  black  lines  are  polynomial 
fits  of  the  measured  values.  The  solid  gray  line  is  the  wake  scaling  law  of  Sreenivasan 
and  Narasimha  (1982),  using  the  momentum  thickness  of  the  far  wake  and  x/C=  1.01  as 
the  hypothetical  origin  of  the  wake. 

Figure  24.  The  wake  half  width  used  in  Figures  19-21,  plotted  against  x/C  for  varying 
Rec  and  for  (a)  Geometry  I  and  (b)  Geometry  II.  The  black  lines  are  polynomial  fits  of 
the  measured  values.  The  solid  gray  line  is  the  wake  scaling  law  of  Sreenivasan  and 
Narasimha  (1982),  using  the  momentum  thickness  of  the  far  wake  and  x/C=0.97  as  the 
hypothetical  origin  of  the  wake. 

Figure  25.  Location  of  boundary  layer  transition,  predicted  from  the  mean  Cp  curve  using 
Thwaites’  method  and  the  one-step  method  of  Michel  (1952).  O,  measured  Cp  value; 

- ,  spline  of  Cp  values;  A,  predicted  transition  location  at  the  Rec  (in  millions) 

indicated. 


Figure  26.  Measured  and  predicted  boundary  layer  momentum  thickness,  0,  at  varying 
Rec.  (a)Thickness  of  the  suction  side  boundary  layers  for  Geometry  I  at  x/C=0.930, 
shown  in  Figure  10.  (b)Thickness  of  the  pressure  side  boundary  layers  for  both  trailing 
edges  at  x/C=l.  Where  available  (Geometry  I,  Rec  of  8M  and  above),  the  values  are 
computed  from  LDV  measurements  at  x/C=l.  Otherwise,  the  values  are  computed  from 
the  wake  data  at  x/0 1.002  shown  in  Figure  18;  •,  computed  from  data  from  Geometry 

I;  O,  computed  from  data  from  Geometry  II;  - ,  prediction  from  the  boundary 

layer  integral  methods. 

Figure  27.  Boundary  layer  parameters  associated  with  separation  for  (a)  Geometry  I  and 
(b)  Geometry  II  at  varying  Rec.  Symbols  give  the  boundary  layer  state  derived  from  PIV 
measurements.  The  boundary  layer  state  at  the  separation  point  (from  Figure  16a)  is 
plotted  as  a  symbol  with  inverse  fill.  Lines  provide  the  fits  (see  Simpson  1989):  —  , 

correlation  for  boundary  layers  separating  from  flat  plates  and  steps; - ,  probability 

of  forward  flow,  y  =  0.8  on  flat  plates  and  steps. 
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ABSTRACT 

High  Reynolds  number  (Re)  wall-bounded  turbulent  flows  occur  in  many  aero-  and 
hydrodynamic  applications.  However,  the  relative  dearth  of  high-Re  experimental  data  has 
hampered  the  development  and  validation  of  scaling  laws  and  modeling  techniques  applicable  to 
such  flows.  This  paper  presents  measurements  of  the  Re-dependent  turbulent  flow  features 
found  near  the  trailing  edge  of  a  two-dimensional  lifting  surface  at  chord-based  Reynolds 
numbers,  Rec,  typical  of  heavy-lift  aircraft  wings  and  full-scale  ship  propellers.  The  experiments 
were  conducted  in  the  William  B.  Morgan  Large  Cavitation  Channel  at  flow  speeds  from  0.25  to 
18.3  m/s  with  a  cambered  hydrofoil  having  a  3.05  meter  span  and  a  2.13  m  chord  that  generated 
60  metric  tons  of  lift  at  the  highest  flow  speed,  Rec  >  50  million  (50M).  The  dynamic  flow  field 
and  surface  pressure  measurements  reported  here  were  made  near  the  foil  s  trailing  edge  which 
could  be  altered  to  have  an  apex  bevel  angle  of  either  44°  or  56°.  Although  generic  turbulent 
boundary  layer  and  wake  characteristics  were  found  at  any  fixed  Rec  in  the  trailing  edge  region, 
the  flow  field  characteristics  were  found  to  be  Re-dependent  because  of  the  variable  strength  of 
near  wake  vortex  shedding.  Analysis  of  planar  particle-imaging-velocimetry  flow-field 
measurements  suggests  that  near  wake  vortex  shedding  is  strongest  when  boundary  layer 
vorticity  shed  at  the  trailing  edge  from  the  pressure  side  of  the  foil  is  most  effectively  rolled-up 
by  the  suction-side  shear  layer  which  separates  upstream  of  the  trailing  edge.  A  scaling  law  for 
this  phenomenon  based  on  Biot-Savart  induction  concepts  for  suction-side  shear-layer  vorticity 
collapses  vortex  shedding  strength  measurements  for  1 ,4M  <  Rec  <  50M  from  both  trailing  edges 
and  from  prior  measurements  on  minimally-lifting  struts  at  Rec  ~  2M. 
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1.  INTRODUCTION 


High  Reynolds  number  (Re)  flows  of  both  human  and  natural  origin  are  commonly 
turbulent  and  the  characteristics  of  the  turbulence  may  depend  on  Re,  particularly  if  the  flow  is 
wall-bounded  as  is  the  case  in  many  aero-  and  hydrodynamic  applications.  Here  the  disparate 
phenomena  associated  with  near-wall  viscous  flow  and  nearly-inviscid  outer-flow  fluctuations 
may  compete  with  each  other  to  cause  subtle  variations  as  Re  increases.  Two  examples  of  such 
mild  Re  variation  at  high  Re  are  the  gradual  decay  of  the  skin-friction  coefficient  as  the 
downstream-distance-based  Re  increases  on  a  smooth  flat  plate,  and  the  decline  of  the  friction 
factor  as  the  bulk-flow  Re  increases  inside  a  smooth-wall  pipe.  These  mild  wall-bounded-flow 
variations  with  increasing  Re  may  influence  fluctuation  levels  in  downstream  free  turbulent 
flows  that  form  when  the  wall-bounded  flow  separates  from  the  bounding  surface.  In  particular, 
shear  layers  and  wakes  formed  downstream  of  separation  points  on  bluff  or  slender  bodies  may 
develop  strong  localized  vortices  that  augment  measured  turbulence  intensity  levels  when  they 
are  present.  The  relative  presence  or  absence,  and  geometrical  arrangement  of  such  strong 
vortices  may  sometimes  be  linked  to  the  characteristics  of  the  upstream  wall  bounded  flow. 

This  paper  presents  an  experimental  study  of  this  phenomenon  at  high  Re.  It  documents 
how  mild  Re-induced  variations  in  the  boundary  layers  on  a  two-dimensional  hydrofoil  govern 
the  strength  of  near  wake  vortex  shedding  at  chord-based  Reynolds  numbers,  Rec  -  U0C/v 
(where  U0  is  the  flow  speed  far  upstream  of  the  body,  C  is  the  chord  length  of  the  surface,  and  v 
is  the  kinematic  viscosity  of  the  fluid)  from  near  1  million  (1M)  to  more  than  50  million  (50M). 
These  experiments  span  the  Re-gap  between  prior  research  studies  ( Rec  up  to  2M  or  so)  and  full- 
scale  applications,  heavy-lift  aircraft  wings  or  ship  propellers  (30M  <  Rec  <  100M).  In 
particular,  the  dynamic  flow-field  and  surface-pressure  measurements  reported  here  were 
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collected  in  the  trailing  edge  region  of  the  hydrofoil  at  low  Mach  number  (Ma).  Here,  the  term 
dynamic  is  intended  to  include  all  unsteady  features  of  the  flow  field  whether  they  occur  at  large-  v 

scale  with  some  organization  (vortex  shedding)  or  involve  a  wide  range  of  length  scales  with  less  .  /_ 

apparent  organization  (turbulence).  Our  companion  paper  (Bourgoyne  et  al.  2003,  hereafter  - 

referred  to  as  Part  I)  reports  time-averaged  flow-field  measurements  for  the  same  hydrofoil. 

Vortex  shedding  in  the  near  wake  of  hydrofoils,  airfoils,  and  turbine  blades  may  cause 
unwanted  hydroacoustic  or  aeroacoustic  noise.  Reduction,  elimination  or  control  of  such  noise 
has  motivated  much  prior  research  —  and  the  current  study  —  into  the  fluid  mechanics  leading  to 
vortex  shedding.  For  hydrofoils,  Blake  (1986)  provides  a  comprehensive  review  of  work  prior  to 
the  early  1980’s.  Included  is  a  summary  of  experiments  that  shows  trailing-edge  vortex 
shedding  to  be  pronounced  on  symmetric  and  blunt  trailing  edge  geometries.  Similar  results  are 
also  found  in  Blake  and  Gershfeld  (1989).  Vortex  shedding  downstream  of  a  blunt  vibrating 
trailing  edge  is  reported  by  Lotfy  and  Rockwell  (1993)  where  additional  contributions  to 
shedding  research  are  listed.  Prasad  &  Williamson  (1997)  review  of  the  vortex  dynamics 
pertaining  to  the  wakes  of  bluff  bodies.  Some  studies  attribute  vortex  shedding  to  absolute 
instability  of  wake  flows.  Overviews  of  global  flow  instabilities  and  blunt-body  wake  flows  are 
provided  by  Huerre  and  Monkewitz  (1990)  and  Oertel  (1990),  respectively 

The  aeroacoustics  of  a  variety  of  lifting  surfaces  has  recently  received  both  experimental 
(Swales  and  Lowson  1997,  Minniti  and  Mueller  1998,  Luire  et  al.  1998,  and  Roger  and  Moreau 
2002)  and  numerical  (Ho  and  Lakshminarayana  1997,  Knight  and  Peltier  1997,  Howe  1999, 

2000,  Manoha,  et.  al.  2000  ,  and  Wang  and  Moin  2000,  2002)  attention.  The  most  recent  of 
these  computational  studies  applies  large  eddy  simulation  to  predict  sound  generated  by  a 
minimally-lifting  strut  (Blake  1975).  However,  the  vast  majority  of  this  aeroacoustic  work  is  for 
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minimally  lifting  bodies  at  Rec  of  less  than  several  million,  below  the  Rec  range  of  full-scale 
applications. 

Similar  studies  of  vortex  shedding  and  wake  turbulence  have  also  been  conducted  to 
understand  the  performance  of  compressor  and  turbine  blades,  and  airfoil  cascades.  Part  I  lists 
prior  studies  in  this  area  that  include  mean  flow  results,  but  these  prior  studies  include  dynamic 
flow  results  as  well.  Literature  reviews  of  unsteady  wakes  from  turbomachinery  blades  are  given 
in  Cicatelli  and  Sieverding  (1995),  and  Ubaldi  et  al.  (1996).  Although  many  studies  of  this  type 
involve  significant  Ma  flow,  such  studies  still  offer  some  insight  into  the  behavior  of  equivalent 
incompressible  flows.  For  example,  Sieverding  and  Heinemann  (1990)  report  that  the  Strouhal 
number  (St)  is  not  strongly  dependent  on  Ma  in  the  subsonic  regime,  but  instead  depends  more 
on  Re  and  boundary  layer  state.  The  more  recent  experimental  work  in  this  area  (Cicatelli  and 
Sieverding  1996,  Rowe  et  al.  2001,  Ubaldi  &  Zunino  2000)  places  a  greater  emphasis  on 
determining  the  frequency  of  shedding  rather  than  on  parameterizing  its  strength,  a  primary 
objective  of  the  present  study. 

The  prior  studies  having  the  greatest  significance  for  the  findings  presented  herein 
involve  extended  (or  streamlined)  bodies  at  low  Ma.  Greenway  and  Wood  (1973)  conducted  an 
experimental  study  of  vortex  shedding  behind  two-dimensional  minimally-lifting  struts  with 
beveled  trailing  edges  at  Rec  ~  1M,  and  report  that  /')  vortex  shedding  can  be  altered  with 
variations  in  the  trailing  edge  apex  angle,  and  ii )  the  state  of  the  separating  boundary  layers  at  the 
trailing  edge  and  their  size  relative  to  their  vertical  spacing  are  important  factors  for  vortex 
shedding.  They  attribute  this  second  observation  to  a  number  of  prior  investigations  including 
that  of  Gerrard  (1966).  Trailing  edge  bevel  angle  changes  and  variable  boundary  layer 
characteristics  are  also  included  in  the  present  study.  Boldman  et  al  (1976)  showed  how  vortex 
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shedding  from  a  blunt  trailing  edge  can  be  reduced  by  unequal  freestream  velocities  on  either 
side  of  a  splitter  plate.  A  similar  ratio  of  freestream  velocities  appears  in  the  proposed  scaling 
parameter  for  vortex  shedding  strength  in  the  present  study  [see  Eq.  (4)]  but  was  not  directly 
investigated.  Blake  (1984)  examined  the  flow  over  and  downstream  of  two-dimensional 
minimally-lifting  struts  with  trailing  edges  very  similar  to  those  in  the  present  study,  and  found 
that  shedding  strength,  as  measured  by  wake  velocity  spectra,  is  affected  by  trailing  edge  bevel 
geometry.  More  recently,  Sieverding  and  Heinemann  (1990)  observed  that  trailing  edge 
geometry  and  Rec  can  have  similar  effects  on  St  in  the  flow  over  a  stmt  at  low  Ma,  and  proposed 
that  the  state  of  the  separated  boundary  layers  determines  the  shedding  process.  The  present 
study  echoes  these  results  at  higher  Rec  on  a  foil  with  significant  lift. 

All  of  this  prior  work  indicates  that  vortex  shedding  is  highly  sensitive  to  the  properties 
of  trailing-edge  boundary  layers  and  shear  layers.  Unfortunately,  this  sensitivity  limits  the 
generality  of  results  for  a  given  geometry  and  Rec .  For  example,  results  for  bluff  bodies  cannot 
be  fully  applied  to  streamlined  bodies  because  the  interaction  between  the  leading-edge 
stagnation  point  and  trailing  edge  flow  dynamics  is  more  subdued  in  for  streamline  bodies. 
Furthermore,  vortex  shedding  results  for  minimally-lifting  streamlined  bodies  cannot  be  directly 
applied  to  lifting  surfaces  because  the  sensitivity  of  vortex  shedding  to  trailing-edge  boundary 
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layer  state  confers  an  importance  to  the  pressure  gradients  associated  with  lift  that  modulate 
boundary  layer  development  upstream  of  the  trailing  edge.  Even  among  lifting  surfaces, 
important  differences  arise  from  variations  in  trailing  edge  geometry,  such  as  the  effect  of  a 
suction-side  trailing-edge  bevel,  a  geometrical  feature  commonly  implemented  on  propeller 
blades  to  ensure  structural  integrity  during  severe  off-design  conditions.  However,  even  beyond 
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issues  of  geometry,  there  is  also  potential  for  Re-dependence  in  aspects  of  boundary  layer 

\  -  -- 

growth,  separation,  and  shear  layer  development  (see  Part  I). 

To  address  these  issues,  the  present  experimental  study  employs  a  large  two-dimensional 
hydrofoil  in  a  large  high-speed  water  tunnel  and  achieves  Rec  values  which  bridge  the  gap 
between  prior  work  and  many  full-scale  applications.  The  model  incorporates  an  application¬ 
relevant  suction  side  surface,  with  a  variable-angle  trailing-edge  bevel,  that  leads  to  section  lift 
coefficients  of  0.52  to  0.55  (see  Part  I).  This  paper  documents  the  dynamic  flow  over  the 
hydrofoil  -  including  turbulence  statistics  and  vortex  shedding  characteristics  -  from  just 
upstream  of  suction  side  boundary  layer  separation  through  the  development  of  the  near  wake  for 
Rec  -  1.4M  to  50M.  In  addition,  a  simple  scaling  relationship,  which  may  be  generally 
applicable  to  trailing-edge  near-wake  flows,  is  deduced  from  the  experimental  measurements. 
This  relationship  successfully  collapses  the  available  shedding  strength  measurements  and  may 
provide  a  means  for  assessing  the  likely  strength  of  vortex  shedding  based  on  foil  geometry  and 
time-averaged  flow  parameters  alone. 

The  remainder  of  this  paper  is  divided  into  five  sections.  The  next  section  (§2)  covers  the 
experimental  setup  and  techniques.  The  third  section  presents  the  flow-field  measurements  that 
document  how  variation  in  Rec  or  trailing-edge  geometry  may  increase  or  decrease  the 
prominence  of  trailing  edge  vortex  shedding.  The  fourth  section' describes  our  attempt  to  extract 
and  quantify  the  mechanism(s)  underlying  vortex  shedding.  The  fifth  section  explains  the 
construction  of  a  scaling  relationship  that  may  have  predictive  capabilities.  The  final  section 
(§6)  summarizes  this  research  effort  and  presents  the  conclusions  drawn  from  it. 


2.  EXPERIMENTAL  SETUP  AND  TECHNIQUES 

The  test  facility,  test  model,  instrumentation,  calibration,  and  estimates  of  uncertainty  in 
the  mean  flow  measurements  are  folly  described  in  Part  I.  This  section  provides  a  brief  summary 
of  that  material  along  with  additional  experimental  details  relevant  to  the  measurements  of  flow 
dynamics  near  the  trailing  edge  of  the  hydrofoil. 

The  hydrofoil  was  tested  in  the  William  B.  Morgan  Large  Cavitation  Channel  (LCC),  a 
low-turbulence  recirculating  water  tunnel  with  a  3.05  m  by  3.05  m  by  13  m  test  section  capable 
of  steady  flow  speeds  (U0)  from  0.25  m/s  to  18.3  m/s.  The  hydrofoil  (C  =  2.134  m  chord,  0.171 
m  maximum  thickness,  depicted  in  Fig.  1)  was  centered  in  and  spanned  the  LCC  test  section  for 
a  6%  blockage  ratio.  Its  cross  section  was  a  modified  NACA-16  suction  side  that  ended  in  a 
trailing  edge  bevel  with  a  44°  (Geometry  I)  or  56°  (Geometry  II)  apex  angle.  The  pressure  side 
of  the  hydrofoil  was  flat  aft  of  28%  chord.  The  foil’s  surface  was  polished  to  a  nominal 
roughness  of  0.25  pm,  and  neither  the  suction  nor  pressure  side  boundary  layers  were  tripped. 
The  foil’s  angle  of  attack  (measured  with  respect  to  the  flat  pressure  side)  for  these  tests  was  0°. 
The  foil  generated  590  kN  (60  metric  tons)  of  lift  with  U0  =  18.3  m/s.  The  LCC  water 
temperature  varied  from  24  °C  to  40  °C  during  these  tests  so  the  Rec  values  quoted  herein  are  for 
the  average  water  temperature  of  (32  °C)  and  have  a  nominal  ±10%  variation. 

Single-point  two-component  Laser-Doppler  velocimetry  (LDV)  and  planar  two- 
component  particle-imaging  velocimetry  (PIV)  were  used  to  measure  fluid  velocity  in  the 
vicinity  of  the  foil’s  trailing  edge.  Seed  particles  for  the  LDV  system  were  silicon  carbide  with  a 
nominal  diameter  of  2  pm.  Seed  particles  for  the  PIV  system  were  16  pm  silver-coated  glass 
spheres.  The  LCC  was  flood  seeded  with  both  types  of  particles.  The  LDV  system  (built  by 
Dantec)  produced  a  long  narrow  focal  volume  (170  pm  diameter,  6  mm  length)  oriented  parallel 
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to  the  foil  span.  The  PIV  system  used  two  side-by-side  LaVision  Flowmaster  1024  x  1280 
digital  cameras,  LaVision  software,  and  two  Spectra-Physics  flash-lamp-pumped  Nd-YAG  lasers 
delivering  800  mJ  at  532  nm.  The  laser  sheet  thickness  was  masked  to  3  mm.  The  imaging  ratio 
of  0.1  mm  to  1  pixel  and  the  32  by  32  interrogation  window  -  typically  capturing  ten  particle 
pairs  -  set  the  PIV  resolution  at  one  fluid  velocity  vector  for  a  cube  3  mm  on  a  side.  The  image 
capture  rate  was  ~1  Hz  and  500  to  2000  images  were  acquired  at  each  flow  condition. 

Uncertainties  in  the  LDV-acquired  and  PIV-acquired  velocity  fluctuations  include  both 
bias  and  random  error.  While  bias  error  in  the  LDV  velocity  fluctuations  is  negligible,  the 
fluctuations  are  presented  as  normalized  quantities,  and  bias  enters  through  the  normalization 
value,  U0.  This  bias  is  estimated  to  be  ±0.2%  of  the  freestream  velocity.  Random  error  is 
negligible  by  comparison.  Uncertainty  in  the  instantaneous  velocity  vectors  from  the  PIV,  and 
the  spatial  location  of  these  velocities,  is  discussed  in  Part  I.  Convergence  of  the  turbulence 
statistics  is  addressed  in  Section  3.3  of  this  paper. 

In  the  present  study,  the  two-component  LDV  was  used  to  acquire  temporal  velocity 
spectra  at  fixed  points  in  the  foil’s  near  wake.  Dantec  software  (BSA  Flow)  was  used  to  estimate 
the  power  spectrum  from  uncorrelated  particle  bursts  having  random  arrival  times.  This  software 
employs  sample-hold  re-sampling,  Hanning  windowing,  and  Fast  Fourier  Transform  (FFT) 
spectrum  calculations.  The  data  re-sampling  has  the  effect  of  a  low-pass  filter.  No  high-pass 
filtering  was  employed,  but  the  maximum  re-sampling  frequency  was  maintained  at  2-3  times  the 
mean  data  rate  as  a  compromise  between  frequency  resolution  and  high-frequency  aliasing. 

In  addition  to  the  velocity  measurements,  dynamic  surface  pressure  measurements  were 
made  with  an  array  of  fifteen  flush-mounted  pressure  transducers  (PCB  138M101)  located  at  the 
PIV  measuring  plane  and  in  the  vicinity  of  the  hydrofoil  trailing  edge.  These  sensors  were 


8 


arrayed  in  an  L-shape  with  lines  of  transducers  set  parallel  and  perpendicular  to  the  flow 
direction.  The  results  presented  in  this  paper  are  from  a  single  transducer  located  on  the  foil  s 
pressure  side  nearest  the  trailing  edge  (x/C,  y/C,  z/S)  =  (0.990,  0,  0.36),  where  the  x  andy  j  •  -gA 
coordinates  are  defined  in  Fig.  1  and  z  completes  a  right-handed  set  of  Cartesian  axes.  The 
/dynamic  pressure  signals  were  analog  band-pass  filtered  from  2  Hz  to  5  kHz  and  sampled  at  10 


kHz. 

Although  the  model  and  its  mounting  scheme  were  design  to  be  as  rigid  (and  durable)  as 
possible,  some  model  vibration  did  occur  and  was  monitored  with  an  array  of  eight 
accelerometers.  A  careful  comparison  of  fluctuation  spectra  recorded  by  the  trailing  edge 
pressure  transducers  and  the  nearest  accelerometer(s)  did  not  reveal  any  significant  correlations 
in  the  frequency  range  of  interest.  Furthermore,  when  the  measured  acceleration  fluctuations 
were  converted  to  root-mean-square  (rms)  surface  velocities,  their  normalized  levels  (0.01%, 


worst  case)  were  found  to  be  well  below  that  of  the  nominal  peak  freestream  turbulence  level  of 
the  tunnel  (0.5%),  and  two  or  more  orders  of  magnitude  below  the  turbulent  velocity  fluctuations 
measured  near  the  foil  surface  (~  4  to  7%).  Thus,  the  measurements  reported  here  are  believed  to 
be  free  of  model  vibration  contamination. 


3.  RESULTS 

Experimental  data  were  collected  for  both  trailing  edge  geometries  at  the  flow  speeds  and 
Reynolds  numbers  listed  in  Table  1.  This  table  also  includes  the  symbols  and  line  types  used  to 
designate  results  at  the  various  Rec  values.  Geometry  I  (44°  trailing  edge  apex  angle)  and 
Geometry  II  (56°  trailing  edge  apex  angle)  results  are  reported  filled  and  open  symbols, 
respectively. 
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Fluid  velocity  measurements  could  be  made  with  either  the  LDV  or  PIV  systems  (or 
both),  and  the  relative  advantages  of  each  were  exploited  to  optimize  test  time  utilization.  The 
PIV  system  provided  instantaneous  two-dimensional  velocity  fields  from  which  spatial 
correlations,  velocity  gradients  (vorticity),  spatial  spectra,  and  time-averaged  results  could  be 
obtained.  However,  the  PIV  frame  rate  was  insufficient  to  measure  the  temporal  evolution  of  the 
flow,  and  its  modest  spatial  resolution  suppressed  small-scale  velocity  fluctuations.  The  LDV 
system  complemented  the  PIV  system,  and  provided  single-point  measurements  with  higher 
spatial  and  temporal  resolution  from  which  turbulence  statistics,  temporal  spectra,  and  time- 
averaged  results  could  be  obtained.  The  agreement  between  the  time-averaged  velocity 
measurements  from  each  system  was  good.  In  addition,  dynamic  surface  pressure  was  measured 
made  near  the  foil’s  trailing  edge. 

The  main  goal  of  this  paper  is  to  present  the  Re  and  geometry  dependence  of  the  dynamic 
flow  features  in  the  vicinity  of  the  foil’s  trailing  edge.  The  following  subsections  cover  the 
attached  boundary  layers  approaching  the  trailing  edge  region  (§3.1),  separating  boundary  layers 
and  the  near  wake  (§3.2),  comparison  of  LDV  and  PIV  results  (§3.3),  near  wake  velocity 
fluctuations  (§3.4),  near-wake  velocity  spectra  (§3.5),  trailing  edge  surface  pressure  spectra 
(§3.6),  and  vortical  flow  features  related  to  wake  velocity  fluctuations  (§3.7). 

3.1  Attached  boundary  layers  approaching  the  trailing  edge  region 

LDV-measured  characteristics  of  the  attached  boundary  layers  approaching  the  trailing 
edge  region  (x/C=0.93)  are  presented  in  Fig.  2,  for  the  (a)  suction  side  and  (b)  pressure  side  of 
the  Geometry  I  hydrofoil  at  Rec  =  1.4M,  4M,  17M,  and  50M.  Here,  U  and  V  are  the  x-  and  y- 
direction  mean  velocities,  u  and  v'  are  the  x-  and  y-direction  velocity  fluctuations,  y0  represents 
the  foil  surface  in  Fig.  2a,  and  Ue  is  the  local  horizontal  flow  speed  at  the  outer  edge  of  each 
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boundary  layer.  All  profiles  represent  turbulent  boundary  layers.  The  suction  side  boundary 
layer  has  higher  fluctuation  levels  at  the  lower  two  Reynolds  numbers.  The  differences  in  the 
fluctuation  levels  are  believed  to  originate  in  the  differing  lengths  of  laminar,  transitional,  and 
turbulent  boundary  layer  flow  on  the  suction  side  of  the  foil  at  each  Re.  The  pressure  side 
boundary  layer  measurements  at  x/C  =  0.93  are  consistent  with  flat  plate  results.  Pressure  side 
boundary  layer  LDV  measurements  were  not  made  at  the  lower  two  Rees 


3.2  Separating  boundary  layers  and  the  near  wake 

An  overview  of  LDV-measured  velocity  fluctuation  levels  with  Geometry  I  is  provided  in 
Fig.  3.  The  three  panels  show  multiple  profiles  of  the  normalized  turbulent  stresses,  u'2/ul , 

v*/ul ,  and  urv'/ul ,  upstream  and  downstream  of  the  trailing  edge,  shown  in  the  left  half  of  each 

panel.  The  gray  vertical  lines  indicate  the  positions  at  which  the  profiles  were  measured.  A 
relative  scale  for  the  levels  is  given  at  the  lower  left  of  each  panel.  The  smoothly-bending  more- 
nearly-horizontal  curves  on  each  panel  indicate  the  outer  edges  of  the  suction-side  and  pressure- 
side  boundary  layers  and  of  the  foil’s  wake,  defined  as  the  locations  at  which  the  streamwise 
velocity  fluctuations  fall  below  the  noise  level  of  the  LDV  measurements.  In  all  cases,  this  noise 
level  (approx.  1%  of  the  freestream  speed)  was  incoherently  subtracted  from  the  reported  u  and 
v'  variances.  These  fluctuation  profiles  have  been  minimally  smoothed  to  improve  graphical 
clarity  using  a  three-data-point  sliding  window.  Results  from  Rec  —  8M  and  50M  are  shown  and 
display  a  mild  Re-dependence  consistent  with  later  downstream  separation  of  the  suction  side 
boundary  layer  at  Rec  -  50M.  Results  for  Rec  -  17M  and  33M  (not  shown)  follow  a  monotonic 
trend  between  the  plotted  results  at  Rec  -  8M  and  50M. 

Similar  overviews  of  the  PIV-measured  velocity  fluctuation  levels  at  Rec=  1-4M,  4M, 
and  50M  for  Geometry  I  and  Geometry  II  are  provided  in  Figs.  4  and  5,  respectively.  The 
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presentation  format  is  the  same  as  Fig.  3,  including  the  xlC  locations  at  which  the  data  are 
presented,  with  one  exception.  Shadowing  by  the  hydrofoil  prevented  PIV  data  acquisition  on 
the  foil’s  pressure  side  and  compromises  the  very-near-wake  data  at  x/C<  1.002  and  y/C  <  0,  so 
the  profiles  plotted  are  from  x/C=  1.002  instead  ofx/C=1.000.  For  Geometry  I,  PIV  data  are  not 
available  at  Rec=  1 .4M  at  the  downstream  station.  For  Geometry  n,  data  are  not  available  at  Rec 
=  50M  at  three  upstream  stations  and  at  x/C  -  1.002  andy/C  <  0. 

As  with  the  LDV  results  shown  on  Fig.  2,  differences  exist  in  the  data  and  are  most 
evident  when  comparing  the  lower  and  higher  Rec  measurements.  The  results  most  different 
from  the  others  are  those  from  Rec  =  1.4M  with  Geometry  I,  and  Rec  =  1.4M  and  4M  with 
Geometry  II.  As  discussed  in  Section  3.4  of  Part  I,  the  pressure  side  boundary  layer  on 
Geometry  I  appears  to  separate  upstream  of  the  trailing  edge  at  Rec  -  1.4M;  at  all  other 
conditions  the  pressure  boundary  layer  separates  at  the  trailing  edge.  Such  early  separation  could 
lead  to  the  unique  downward  shift  of  the  Rec  =  1.4M  Geometry-I  pressure-side  fluctuation 
profiles.  For  Geometry  II,  the  elevated  turbulence  levels  at  the  lower  two  Rec  arise  from  strong 
vortex  shedding  at  these  conditions.  In  particular,  Fig.  5b  shows  vaju\  to  be  much  higher  at 

these  two  conditions,  particularly  in  the  region  downstream  of  the  trailing  edge  near  y  =  0.  Such 
elevated  vertical  velocity  fluctuations  figure  prominently  in  a  simple  scaling  analysis  of  vortex 
shedding  presented  in  §5. 1 .  At  Rec  >  17M,  the  PIV-measured  fluctuation  profiles  show  the  same 
mild  Re-dependence  found  in  the  LDV  data  (see  Fig.  2).  However,  the  PIV-data  at  Rec  =  17M 
and  33M  are  so  similar  to  those  at  50M  that  they  were  omitted  from  Figs.  4  and  5  for  graphical 
clarity 
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3.3  Comparison  of  LDV and PIV results 

A  comparison  of  the  Geometry  I  LDV  and  PIV  results  at  Rec  =  50M  (see  Figs.  2  and  4) 
reveals  similar  trends  and  profile  shapes  but  disagreements  in  magnitude.  The  differences  are 
attributed  to  differing  levels  of  statistical  convergence,  processing  artifacts  in  the  PIV 
measurements,  and  the  superior  spatial  resolution  of  the  LDV  measurements.  The  statistical 
convergence  of  each  measurement  was  set  by  a  compromise  between  test  time  at  each  condition 
and  the  total  number  of  conditions  tested.  Here,  the  longest  fluctuation  time  scale  in  the  foil’s 
near  wake  flow  is  the  vortex  shedding  period,  T  =  l/fp  /Ua  (Blake  1986)  where  A™m  is 

the  nominal  near  wake  thickness.  Single-point  LDV  measurements  typically  involved  500  to 
12,000  particle  bursts  in  a  sampling  period  of  0.6  min.,  approximately  600 T  to  3,600T  at  Rec  = 
8M  and  50M,  respectively.  For  the  PIV  measurements,  1,000  to  2,000  images  were  recorded  at 
each  test  condition  in  multiple  runs  of  80  images  each  over  a  test  time  of  30  to  60  minutes,  a 
sampling  period  of  ~1057\  Figure  6  shows  sample  convergence  histories  for  the  PIV 
measurements  at  Rec  =  50M,  and  these  results  indicate  that  2,000  PIV  images  yield  well 
converged  mean  values  (Fig.  6a)  and  adequately  converged  turbulence  statistics  (Fig.  6b). 

The  poorest  agreement  between  the  LDV  and  PIV  flow  fluctuation  measurements 
occurred  at  Rec  =  50M  in  the  profiles  at  x/C  —  1.0094  (see  Fig.  7).  In  general,  the  PIV 
measurements  fall  below  the  LDV  results  because  of  the  additional  spatial  filtering  associated 
with  the  larger  PIV  interrogation  volumes.  However,  the  peak  pressure-side  PIV  results  for 
fiju1  and  u'v'jul  do  occasionally  exceed  the  LDV  values.  This  occurs  at  this  location  because 
the  PIV  interrogation  volume  spans  a  region  of  high  streamwise-velocity  gradient.  A  random 
distribution  of  particles  in  such  an  interrogation  volume  and  the  recursive  PIV  image-processing 
software  tend  to  produce  extreme  —  as  opposed  to  spatially  averaged  -  instantaneous  horizontal 


velocities,  thus  raising  the  fluctuation  levels  at  this  location.  Fortunately,  this  problem  appears  to 
be  confined  to  1.00  <  x/C  <  1.01  and  y  ~  0.  Figure  8  provides  a  further  quantitative  comparison 
between  the  LDV  and  PIV  results.  Here,  the  ratios  of  LDV-measured  fluctuations  divided  by 
PIV-measured  fluctuations  are  provided  vs.  Rec  for  the  peak  suction-side  values  of  u,2/u20 , 

v^/ul ,  and  tfv'/ul  at  x/C  =  1.0094.  For  these  suction  side  ratios,  the  LDV  measurements 
produce  consistently  higher  fluctuation  levels  and  PIV-measurement  spatial  averaging  in  both 
the  streamwise  and  vertical  directions  must  be  considered  in  order  to  quantitatively  account  for 
the  differences  shown  on  Fig.  8.  Nevertheless,  the  agreement  between  LDV  and  PIV 
measurements  improves  with  decreasing  Rec,  as  would  be  expected  if  their  differences  are 
related  to  varying  resolution  of  the  energy-containing  flow  scales. 

The  remainder  of  this  paper  relies  on  the  PIV  measurements  for  which  more  conditions 
were  measured  and  from  which  the  spatial  characteristics  of  the  instantaneous  flow  fields  can  be 
deduced.  However,  these  results  must  be  interpreted  in  light  of  the  spatial  averaging  inherent  in 
the  PIV  measurements. 

3.4  Near-wake  velocity  fluctuations 

The  Re  and  geometry-dependence  of  the  near  wake  shown  in  Figs.  4b  and  5b  may  be 
illustrated  concisely  by  plotting  the  curves  of  the  y-location  at  which  the  fluctuation  levels  are 
maximum  for  a  given  x-coordinate,  and  the  fluctuation  levels  recorded  along  such  curves.  The 
results  for  the  normalized  vertical  velocity  variance,  v'2/t/2,  are  given  in  Figs.  9  and  10  for 
Geometry  I  and  II,  respectively.  Values  for  1.000  <  x/C  <  1.005  are  not  plotted  because  the 
suction  and  pressure  sides  peaks  of  v'2/u2  are  comparable  and  consistent  selection  of  a  single 
peak  is  problematic.  In  Fig.  9a,  the  Rec  =  1.4M  case  is  clearly  different,  and  this  is  again 
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attributed  to  pressure  side  boundary  layer  separation  prior  to  the  trailing  edge.  Results  at  the 
other  Rees  are  more  similar,  though  close  examination  of  Fig.  9b  reveals  some  Re-variation  in 
the  downstream  location  of  the  maximum  fluctuations.  By  comparison,  dramatic  Re-variation  in 
the  location  of  maximum  fluctuations  is  shown  for  Geometry  II  in  Fig.  10.  Here,  the  point  of 
maximum  fluctuations  resides  near  x/C=  1.03  at  Rec=  1.4M,  progresses  upstream  towards  the 
trailing  edge  as  Rec  is  increased  to  4M,  and  then  progresses  downstream  again  when  Re 
increases  to  17M  or  higher.  A  comparison  of  Figs.  9b  and  10b  also  shows  enhanced  vertical 
velocity  fluctuations  occur  with  Geometry  II. 

Figures  11  and  12  follow  the  same  format  as  Figs.  9  and  10,  but  give  the  downstream 
evolution  of  the  normalized  Reynolds  shear  stress,  -«v/l/2 ,  taken  from  Figs.  4c  and  5c.  The 

stress  profiles  have  two  extremes  of  opposite  signs,  and  the  coordinates  and  values  of  both  are 
plotted.  Since  the  locations  of  peak  Reynolds  stress  mark  the  shear  layers  centers,  Figs.  1  la  and 
12a  give  an  indication  of  the  evolution  of  the  wake  thickness.  As  found  by  other  researchers 
(Blake  1986),  the  vertical  separation  of  the  shear  layers  initially  decreases  but  then  increases,  and 
the  location  of  closest  approach  coincides  with  the  downstream  location  of  the  Reynolds  shear 
stress  extremes.  As  with  Fig.  9a,  Fig.  11a  shows  Geometry  I  at  Rec-  1.4M  as  a  distinguished 
case.  The  initiation  of  the  peak  pressure-side  Reynolds  shear  stress  curve  aty/C  =  -0.002  is 
consistent  with  early  pressure  side  separation.  In  Fig.  12  for  Geometry  II,  as  in  Fig.  10,  the 
downstream  distance  at  which  the  peak  stresses  are  encountered  vary  significantly  with  Rec- 
And  again,  comparison  of  Figs.  1  lb  and  12b  shows  enhanced  shear  stress  occurs  in  the  pressure 
side  shear  layer  with  Geometry  II. 

Figures  13,  14  and  15  show  the  near  wake  evolution  using  velocity  profiles  at  xlC  = 
1.002,  1.009,  1.028,  and  1.047  for  Rec=  1.4M,  4M,  17M,  and  50M.  The  data  are  replotted  from 
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Figs.  4  and  5,  but  here  have  been  scaled  as  a  plane  wake:  the  horizontal  axes  are  normalized  by 
AC/max  where  AC/  =  U(y)  -  Uess,  (Uess  is  the  suction  side  wake-exterior  horizontal  velocity)  and 
the  vertical  axes  are  (y  -  yctr)/ym  ,  where  yctr  denotes  the  average  of  the  vertical  coordinates 
where  A U  =  At/max/2.  The  full  wake  width,  2ym,  is  that  measured  at  AU=  AC/max/2.  These  same 
coordinates  were  used  for  the  mean  flow  profiles  presented  in  Part  I,  wherein  the  values  for  yctr, 
AC/max,  and  ym  are  provided.  On  the  right-most  frames  of  Figs.  13  and  15,  fluctuation  profiles 
for  the  self-similar  wake  of  a  symmetric  airfoil  at  Rec  ~  104  (Wygnanski  et  al.  1986)  are  plotted 
for  comparison.  The  plane  wake  coordinates  are  effective  in  collapsing  the  fluctuation  results 

between  the  two  trailing  edge  geometries  at  x/C  =  1.002  except  for  u'2  at  Rec-  4M  (Fig.  13b). 
However,  as  the  wakes  evolve  downstream,  fluctuation  profile  variations  between  the  two 
trailing  edge  geometries  begin  to  increase,  particularly  at  the  lower  two  Rec.  Most  striking  is 
Rec  =  4M,  for  which  the  fluctuations  associated  with  the  pressure  side  shear  layer  grow  more 
rapidly  on  Geometry  II  than  on  Geometry  I.  This  is  noteworthy  because  the  normalized  pressure 
side  boundary  layers  are  nearly  independent  of  Rec ■  and  trailing  edge  geometry  above  Rec  =  4M. 
This  unusual  near  wake  behavior  at  Rec  =  4M  with  Geometry  II  is  caused  by  increased  vortex 
shedding  in  the  near  wake,  a  phenomenon  linked  to  trailing-edge  geometry  and  suction  side 
shear  rate  in  §5.1. 

Use  of  wake  coordinates  in  Figs.  13-15  also  serves  to  emphasize  Re-dependence  in  the 
foil’s  near  wake  evolution.  Most  notably  at  x!C=  1.047,  the  normalized  strength  of  the 
fluctuations  decreases  by  more  than  a  factor  of  two  with  increasing  Rec ■  While  some  of  this 
variation  is  undoubtedly  due  to  Re-dependence  in  the  PIV  spatial  averaging  (see  §3.3,  and  Fig. 
8),  the  remainder  appears  to  be  a  genuine  Re-effect.  Furthermore,  the  present  results  at  the 
lowest  Rec  of  this  study  (1.4M)  are  the  best  match  with  the  prior  measurements  in  a  foil  wake  at 
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Rec  ~  104.  However,  in  all  cases,  the  velocity-fluctuation  and  Reynolds-stress  profiles  have  not 
achieved  self-similar  form  by  x/C  =  1.047  (the  furthest  downstream  location  presented  on  Figs. 
13-15)  nor  by  x/C  =1.11  (the  furthest  downstream  location  at  which  such  measurements  were 
made,  not  presented  here). 


3.5  Near -wake  velocity  spectra 

Based  on  the  turbulence  fluctuation  results,  vertical  velocity  spectra  were  measured  in  the 
foil’s  near  wake  to  aid  in  identifying  the  relative  contributions  of  nearly-periodic  flow 
oscillations  (vortex  shedding)  and  more-nearly-random  fluctuations  (turbulence).  Temporal 
spectra  of  v'  were  acquired  at  two  downstream  positions  in  the  near  wake  using  the  LDV  system. 
Spatial  spectra  of  v'  were  computed  from  the  PIV  measurements  along  a  line  segment  that  passes 
through  or  near  the  temporal-spectra  measurement  points. 

Figures  16  and  17  present  the  LDV-measured  vertical  velocity  spectra  in  the  foil’s  near 
wake  at  jc/C=1.009  and  1.070,  respectively,  for  Geometry  I  (Figs.  16a,  17a)  and  Geometry  II 
(Figs.  16b,  17b)  at  Rec  =  1.4M,  4M,  17M,  and  50M.  The  y-coordinate  for  each  spectral 
measurement  was  chosen  by  searching  for  the  maximum  spectra  peak  height,  at  the  given  x/C, 
over  the  range  -0.025  <  y/C  <  +0.025.  In  all  cases,  the  maximum  spectral  peak  height,  and 
maximum  spectral  peak  area,  were  found  in  the  immediate  proximity  of  the  pressure  side  shear 
layer  with  lesser  spectral  peaks  typically  associated  with  the  suction  side  shear  layer.  The 
coordinates  of  the  peak  spectral  activity  are  provided  in  Table  2  and  are  nearly  coincident  with 
the  loci  of  maximum  v,2/ul,  shown  in  Figs.  9  and  10.  The  temporal  spectra,  <£„(/)>  were 

normalized  so  that  v*/u20  =  where  N°m  is  a  constant  nominal  wake  thickness 

of  0.0 1C,  <!>„(/)  is  the  normalized  spectra,  and  /  is  frequency  in  Hz.  The  dimensionless 

frequency  on  the  horizontal  axes  on  Figs.  16  and  17  may  also  be  interpreted  as  a  normalized 
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wavenumber,  kAn°m  ■  fA™a/U0 .  Only  spectral  values  below  the  effective  Nyquist  frequency  set 

by  the  LDV  data  rate  are  plotted.  The  data  from  Rec=  50M  was  not  as  fully  converged  and  was 
smoothed  for  graphical  clarity.  This  smoothing  had  negligible  effect  on  peak  characteristics. 

Figure  18  presents  the  PIV-measured  vertical  velocity  spectra  in  the  foil’s  near  wake 
measured  on  the  line  segment  from  ( x/C,y/C)  =  (1.000,  0.009)  to  (1.100,  -0.005)  for  Geometry  I 
(Fig.  18a)  and  Geometry  II  (Fig.  18b)  at  Rec=  1.4M,  4M,  17M,  and  50M.  These  spatial  spectra, 
®v(k),  were  computed  from  between  500  and  1000  PIV  fields,  and  were  normalized  so  that 

v*/u2  =  A'"”"  f*°°<bdk,  where  <tv(k)  is  the  normalized  spatial  spectrum,  and  k  is  the  wavenumber 

in  cycles/m-1.  Unfortunately,  these  spatial  spectra  are  potentially  contaminated  at  low 
wavenumbers  by  two  factors  that  do  not  influence  the  low  frequency  portion  of  the  temporal 
spectra.  First,  the  foil’s  near  wake  lacks  a  statistically  homogeneous  direction  in  which  to 
compute  spatial  spectra.  And  second,  the  limited  spatial  extent  of  the  PIV  data  prevents 
adequate  spectral  estimation  at  the  wavelengths  typical  of  vortex  shedding.  Although  these 
factors  limit  the  utility  of  these  spatial  spectra  in  the  study  of  vortex  shedding,  the  spectra  should 
be  unbiased  above  kA™m  =  10~’  or  so.  The  downward  sloping  lines  on  Figs.  16-18  have  a  slope 

of -5/3  and  are  provided  for  reference. 

The  vertical  velocity  spectra  show  several  common  features.  At  low  frequencies 
(/A™m /U0  less  than  ~0.1)  both  trailing  edges  produce  flat  spectral  levels  as/decreases.  And,  at 

high  frequencies  and  wavenumbers  {fAnoymjlJ0  and  kAn°m  greater  than  ~0.5)  the  expected  power- 

law  form  for  turbulent  fluctuations  is  recovered  within  the  limitations  of  the  measurements.  In 

between,  every  spectrum  displays  a  peak  -  to  greater  or  lesser  degree  -  that  is  associated  with  the 

strength  of  vortex  shedding.  In  the  temporal  spectral,  these  peaks  occur  at  different  normalized 

frequencies,  with  differing  heights  and  widths,  but  always  near  the  anticipated  Strouhal  number, 
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St  =  fpA™m /Ua  — 1/2 Jt,  for  near-wake  vortex  shedding  from  foils  and  struts  (Blake  1986).  In 
particular,  Geometry  H  produces  higher  spectral  peaks  than  Geometry  I  at  every  Rec,  with  the 
highest  spectral  peaks  occurring  at  Rec=  4M  for  both  trailing  edge  geometries. 

For  graphical  clarity,  Figs.  16-18  show  results  from  only  four  of  the  seven  Rees  of  this 
study.  Using  the  parameters  defined  in  Fig.  19,  vortex-shedding-peak  results  from  the  temporal 
spectra  at  ail  seven  Rees  are  presented  in  Figs.  20  and  21  for  xIC=  1.009  and  1.070,  respectively. 
The  last  panel  on  these  figures  gives  the  frequency  of  the  spectral  peak,^,  in  the  dimensionless 
form  laSt.  Some  of  the  weaker  spectral  peaks  lack  a  distinctive  maximum  so  the  frequency  of 
the  centroid  of  the  peak  was  substituted  for  fp  in  these  cases.  To  ease  comparisons  between 
conditions,  the  normalized  frequency  f  -  f /  fp  was  used  in  all  calculations  leading  to  Figs.  20 

and  21,  so  that  V2ju20  ■ 

The  Rec  and  geometry  dependencies  noted  on  Figs.  16-18  are  further  evident  in  Figs.  20 
and  21.  The  spectral  peak  for  Rec  =  4M  with  Geometry  II  has  the  greatest  magnitude  and  the 
minimum  width.  Comparison  of  Figs.  20  and  21  also  reveals  aspects  of  the  downstream 
evolution  of  the  wake.  The  peak  height  expressed  in  terms  of  the  energy  in  the  upstream  velocity 
(the  normalization  depicted)  is  reduced  at  x/C  =  1.070  in  comparison  to  x/C  =  1.009.  However, 
peak  height  compared  to  the  energy  at  all  frequencies,  obtained  by  dividing  the  peak  values  of 
panels  (b)  by  the  variance  given  in  panel  (a),  is  actually  increased  at  x/C  =  1 .070  in  comparison 
to  x/C=1.009.  This  indicates  that  a  greater  fraction  of  the  total  fluctuations  occur  at  the  shedding 
frequency  as  the  flow  evolves  downstream.  This  trend  and  the  diminished  fluctuation  level 
suggests  that  smaller  vortical  structures  are  dissipated  or  progressively  grouped  into  larger  -  and 
more  slowly  moving  -  wake-scale  structures  in  the  manner  expected  for  a  developing  two- 
dimensional  vortex  street. 
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3.6  Trailing  edge  surface  pressure  spectra 

The  shedding  trends  in  the  wake  velocity  spectra  agree  with  the  independently  measured 
temporal  spectra,  <hp( /),  of  surface  pressure  fluctuations  on  the  trailing  edge.  Figure  22 

presents  results  from  a  flush  mounted  dynamic  pressure  transducer  located  on  the  foil’s  pressure 
side  at  x!C  =  0.99.  These  spectra  are  normalized  in  a  manner  similar  to  the  velocity  spectra: 
p^/<202  b  A™m /U0  f^&pdf ,  where  <tp  is  the  normalized  spectrum,  Qo  »  pU2Jl,  and  p  is  the 

density  of  water  at  32  °C.  The  normalized  variances  computed  from  the  integral  of  the  power 
spectra  are  less  than  0.0001  for  Rec  >  3M  for  both  trailing  edges.  Each  plotted  spectrum  on  Fig. 
22  was  calculated  from  approximately  3  million  data  points,  partitioned  with  50%  overlap  to 
provide  376  spectral  windows  of  214  (16,384)  data  points,  with  an  estimated  spectral  amplitude 
uncertainty  of +/-12%  (Vetterling,  et.  al.,  1992).  A  linear  least-squares  fit  was  subtracted  from 
each  partition  to  remove  any  drift  and  to  zero  the  partition  mean.  Neither  acceleration 
contamination  nor  noise  contributions  have  been  removed  from  these  spectra,  and  a  1/1 00th 
decade  filter  has  been  applied  to  clarify  the  high-frequency  portion  of  these  spectra  (with 
negligible  effect  on  peak  magnitudes).  The  upward  translation  of  the  curves  on  Fig.  22  as  Rec 
falls  below  8M  reflects  the  increasing  fraction  of  the  transducer  output  resulting  from  noise  that 
does  not  scale  with  Q0.  The  diagonal  lines  below  the  data  on  Fig.  22  have  slope  of-1. 

Two  fluid  dynamic  phenomena  are  evident  in  the  plotted  pressure  spectra.  The  surface 
pressure  fluctuations  produced  by  an  attached  turbulent  boundary  layer  (TBL)  lead  to  broadband 
unpeaked  spectra  (see  Willmarth  1975,  or  Farabee  and  Casserella  1991)  having  a  shallow 
spectral  slope  slightly  above  -1  (Goody  2002),  and  such  a  spectral  signature  is  clearly  evident  in 
all  the  results  on  Fig.  22.  In  addition,  pressure  fluctuations  induced  by  organized  near  wake 
vortex  shedding  lead  to  narrow  spectral  peaks  at  f=fp,  and  such  peaks  are  evident  in  the  spectra 
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measured  for  Geometry  II.  To  be  apparent,  the  vortex-shedding  component  must  be  stronger 
than  the  TBL  fluctuations  at  the  same  frequency.  The  pressure  spectra  for  Geometry  I  in  Fig. 
22a  lack  significant  peaks,  indicating  that  vortex  shedding  was  not  detectable  above  the  TBL 
fluctuations.  However,  peaks  in  the  pressure  spectra  are  found  with  Geometry  II  for  Rec  <  8M 
(see  Fig.  22b),  with  a  maximum  detected  shedding  condition  at  Rec  —  4M.  Even  though  the 
curve  for  Rec  =  3M  on  Fig.  22b  is  compromised  by  poor  signal-to-noise  and  suffers  from  a 


strong  noise  peak  near  /An;m/t/0  -  0.06,  it  does  indicate  a  vortex  shedding  peak  near  the  correct 

frequency.  Pressure  fluctuation  data  from  speeds  below  Rec  =  3M  are  dominated  by  noise  and 
are  not  presented.  Overall,  given  the  threshold  effect  of  the  TBL  pressure  fluctuations,  these 
pressure  spectra  are  fully  consistent  with  the  vertical  velocity  spectral  results  reported  on  Figs. 
16-21. 

3. 7  Vortical  flow  features  related  to  the  velocity  fluctuations 

The  spatial  PIV  velocity  measurements  allow  the  size,  strength,  and  spacing  of  vortices  in 
the  foil’s  near  wake  to  be  investigated.  Visualization  of  vortices  and  regions  of  shear  flow  may 
be  achieved  directly  from  computations  of  vorticity.  However,  for  this  presentation  of  results, 
the  swirling  strength  (Zhou  et  al.  1996,  1999,  Adrian  et  al.  2000)  is  used  as  a  more  selective 
method  for  visualization  of  vortices  alone.  Swirling  strength,  x  (with  units  of  t  ’like  vorticity),  is 
the  absolute  value  of  the  imaginary  part  of  the  eigenvalues  of  the  two-dimensional  velocity 
gradient  tensor.  It  specifies  the  orbital  rotation  rate  of  a  fluid  element  centered  on  the  point  of 
interest,  but  is  zero  in  pure  shear  flow. 

Figure  23  presents  grayscale  images  of  normalized  swirling  strength  computed  from 
selected  instantaneous  velocity  fields  for  Geometry  I  at  Rec  -  4M  and  for  Geometry  13  at  Rec  = 
1.4M,  4M,  17M,  and  50M.  Here  clockwise  rotation  has  been  treated  as  positive  (white).  The 
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requisite  velocity  derivatives  are  computed  on  a  1.6  mm  grid  from  a  central  difference  stencil  of 
width  3.2  mm  (the  length  scale  of  the  PIV  interrogation  area).  The  swirling  strength  is  then 
computed,  spatially  smoothed  with  a  box-car  filter  of  3.2  mm  x  3.2  mm,  and  thresholded  to 
remove  low-level  values.  To  better  visualize  the  weaker  suction  side  swirling  strength,  the  gray 
levels  used  in  Fig.  23  are  not  symmetric  about  zero.  In  addition,  the  instantaneous  fields  chosen 
best  illustrate  the  prevailing  differences  between  the  five  test  conditions  shown.  Only  Rec-  4M 
is  presented  for  Geometry  I  (Fig.  23a)  because  results  at  other  Rec  for  this  geometry  are  similar 
in  appearance.  The  instantaneous  velocity  field  used  to  compute  the  swirling  strength  in  Fig.  23c 
was  also  used  to  compute  the  vorticity  of  the  field  shown  in  Fig.  24a;  note  the  similarity  of  these 
two  fields  for  x/C  >  1.02. 

In  Fig.  23  the  frequency  of  the  most  energetic  velocity  fluctuations  is  apparent  in  the 
spacing  of  the  largest  pressure  side  vortices.  For  example,  observe  the  three  largest  pressure  side 
vortices  for  Rec  =  4M  with  Geometry  II  in  Fig.  23c.  The  streamwise  spacing  of  the  vortices  at 
x/C=  1.045  and  x/C  =  1.10  is  approximately  6A"°m ,  which  gives  the  expected  Strouhal  number 

of  approximately  l/2nr.  The  greater  degree  of  organization  at  Rec=  4M  for  Geometry  II  (c)~ 
compared  to  Geometry  I  (a)  at  the  same  Rec  reflects  the  geometry  dependence  seen  in  the  spectra 
ofFigs.  16-18. 

The  Re-dependence  of  the  Geometry  II  flow  seen  in  the  velocity  spectra  (Figs.  17b,  18b) 
is  also  evident  in  the  swirling  strength  fields,  Fig.  23b-e.  Most  obvious  is  the  variation  with  Re 
in  the  streamwise  spacing  of  the  strongest  pressure  side  vortices,  a  feature  clearly  linked  to  Re 
variation  in  St.  The  pressure  side  vortex  spacing  at  Rec  —  4M  matches  2nSt  ~  1,  while  the  other 
Re  have  more  closely  spaced  vortices  implying  a  2 itSt  >1.  In  particular,  the  Rec  -  50M 
condition  in  panel  (e)  shows  two  predominant  values  of  vortex  spacing.  Here,  three  groups  of 


vortices  have  a  spacing  of  several  wake  thicknesses,  while  vortices  within  these  groups  are 
spaced  more  closely.  This  closer  spacing  may  reflect  vortex  roll-up  on  a  length  scale  determined 
by  the  pressure  side  boundary  layer  instead  of  the  larger  wake  length  scale. 

In  addition,  Fig.  23  shows  near  wake  conditions  for  which  the  pressure-side  time- 
averaged  shear  layers  (Part  I,  Fig.  18)  are  very  similar  at  the  trailing  edge  but  then  go  on  to 
develop  into  different  wake  vortex  distributions.  For  example  at  Rec  =  8M  and  50M  for 
Geometry  II,  the  flow  characteristic  that  differs  most  between  these  two  conditions  is  the  state  of 
the  suction-side  shear  layer  at  the  trailing  edge,  which  varies  with  Rec  because  of  Re-dependent 
suction-side  boundary  layer  development.  As  a  second  example,  the  flows  with  Geometry  I  (Fig. 
23a)  and  Geometry  II  (Fig.  23c)  at  Rec=  4M  have  nearly  identical  pressure-side  shear  layers  at 
the  trailing  edge  and  very  similar  suction  side  shear  layers,  yet  these  conditions  also  differ  in 
pressure-side  vortex  distribution.  In  this  case,  the  flow  characteristic  that  differs  most  between 
the  conditions  is  the  wake  thickness,  as  set  by  the  different  trailing  edge  geometries.  These 
observations  suggest  that  the  suction  side  shear  layer  and  the  wake  thickness  both  affect  the  roll¬ 
up  of  the  pressure  side  shear  layer  into  wake  vortices.  - 

4.  Discussion:  Analysis  of  Vortex  Fields 

The  turbulence  fluctuation  profiles  of  the  evolving  wake  display  both  Re  and  geometry 

dependencies,  particularly  in  the  vertical  velocity  fluctuations  at  the  lower  two  Rec ■  The  spectra 

of  the  vertical  velocity  and  of  the  trailing  edge  surface  pressure  fluctuations  link  this  Re  and 

geometry  dependence  to  vortex  shedding.  PIV-acquired  instantaneous  swirling-strength  fields 

confirm  this  connection.  These  facts  motivate  an  analysis  of  the  Re  and  geometry-dependence  of 

the  vorticity  distributions  in  the  wake  and  the  relationship  between  identifiable  vortices  and 

velocity  fluctuations.  Specifically,  the  analysis  presented  here  attempts  to  separate  the  near- 
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wake  dynamic  contribution  of  strong  large-scale  vortices  from  the  background  of  weaker 
smaller-scale  turbulence.  The  goal  of  this  effort  is  develop  a  physical  model  of  vortex  shedding 
that  can  indicate,  or  even  predict,  when  near-wake  vortex  shedding  will  be  strong.  This  analysis 
is  based  on:  assessing  of  the  mean  flux  of  vorticity  into  the  wake  (§4.1),  identifying  vortices  in 
the  wake  (§4.2),  determining  vortex  characteristics  (§4.3),  establishing  vortex  locations  (§4.4), 
and  reconstructing  velocity  fluctuations  from  the  identified  vortices  alone  (§4.5).  The  physical 
model  and  the  scaling  laws  derived  from  it  are  presented  in  §5. 

4.1  Mean  flux  of  vorticity  into  the  wake 

The  vorticity  in  the  wake  comes  from  the  foil’s  suction  and  pressure  side  boundary 
layers.  Assuming  a  negligible  correlation  between  instantaneous  vorticity  fluctuations  and 
instantaneous  foil-surface-tangent  velocity  fluctuations,  the  time-averaged  circulation  flux  into 
the  near  wake  may  be  computed  from  the  time-averaged  boundary  layer  profiles  (see  Fig.  2,  and 

Part  I  Fig.  18).  Here,  the  time-averaged  flux  of  circulation  is  defined  as  fa  J^U^dh,  where 
Q  =  dUjdh\  Fig.  1  defines  the  surface-normal,  h,  and  surface-tangent,  t,  coordinates;  he  specifies 
is  the  edge  of  the  boundary  layer;  and  Ut  is  the  mean  tangent  velocity.  The  computed  circulation 
fluxes  are  presented  in  Table  3  as  ratios  of  the  non-dimensional  circulation  flux,  f1  jU 0 ,  at  a 
given  condition,  normalized  by  the  circulation  flux  value,  t  =  {^ps/^l)Rec_^M>  fr°m  the  pressure 

side  boundary  layer  for  Geometry  II  at  Rec  =  4M.  The  value  of  f  *  is  given  at  the  bottom  of 
Table  3.  Unfortunately,  the  measured  data  were  insufficient  to  estimate  suction  side  circulations 
fluxes  with  Geometry  II.  Mismatch  in  the  pressure  and  suction  side  circulation  fluxes  with 
Geometry  I  are  presumed  to  be  caused  by  the  neglected  velocity-vorticity  fluctuation  correlations 


and  the  fact  that  the  suction  side  flux  was  calculated  at  x/C  =0.93,  well  upstream  of  the  foil’s 
trailing  edge. 

Surprisingly,  the  strongest  vortex  shedding  (Rec=  4M)  occurs  with  the  lowest  normalized 
pressure-side  circulation  flux  for  both  trailing  edges.  For  Geometry  I,  where  suction-side  data 
are  available,  the  Rec=  4M  case  does  not  have  the  lowest  or  the  highest  circulation  flux.  From 
these  observations,  it  is  concluded  that  circulation  flux  from  the  attached  boundary  layers  is  not 
by  itself  a  predictor  of  the  strength  of  vortex  shedding.  This  conclusion  holds  even  when  the 
local  external  velocity,  Ue,  is  used  in  lieu  of  U0  for  normalization.  This  conclusion  is  reasonable 
because  the  magnitude  of  f  accounts  only  for  the  amount  of  circulation  entering  the  wake 
without  addressing  the  distribution  of  that  circulation.  In  the  near  wake,  the  circulation 
distribution  will  be  determined  by  the  characteristics  of  individual  vortices  and  their  spatial 
organization. 

4.2  Vortex  identification 

Extracting  vortex  characteristics  from  the  PIV  field  measurements  involved  several 
processing  steps.  First  the  PIV  fields  were  sifted  for  strong  vortices  by  computing  the  vorticity, 
co,  via  centered  spatial  first  differences  (see  §3.7),  thresholding  the  computed  vorticity 
magnitudes,  identifying  the  remaining  vorticity  peaks  as  candidate  vortices,  and  ranking  these 
candidate  vortices  from  highest  to  lowest  vorticity  magnitude.  Then,  starting  with  the  strongest 
vorticity  peak,  each  candidate  vortex  was  assigned  a  radius  by  starting  from  its  peak-vorticity 
location  and  repetitively  computing  the  spatially-averaged  vorticity  within  circles  of  increasing 
radius  until  the  peak-divided-by-average  vorticity  ratio  reached  1.77.  When  applied  to  an 

idealized  vortex  with  a  Gaussian  vorticity  distribution,  co(r )  a  oomaxe~^a)  and 

T(r)  =  omMJto2 (\  -  e~^a)1 )  with  r  as  a  radial  coordinate,  this  method  yields  a  core  radius  of 
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1.121  a,  which  is  the  radius  at  which  the  vorticity-induced  velocity  is  maximum.  After  this  core 

sizing  was  completed,  the  candidate  vorticity  peak  was  idealized  as  a  Gaussian  vortex  having  the 
\  11 
same  core  diameter  and  far-field  circulation.  Before  continuing  on  to  the  next  candidate  vortex, 

the  vorticity  put  into  the  new  idealized  Gaussian  vortex  was  subtracted  from  the  PIV  vorticity 

field  and  the  threshold  was  reapplied  to  the  remaining  candidate  vortices.  In  this  way,  field 

vorticity  was  not  double-counted  and  multiple  vorticity  peaks  in  close  proximity  were  absorbed 

into  a  single  idealized  vortex.  For  this  procedure,  the  vorticity  threshold  was  chosen  to  extract 

approximately  33%  of  the  spatially-averaged  vorticity  magnitude  at  each  flow  condition.  This 

threshold  typically  identified  15  to  20  vortices  in  a  given  PIV  field. 

Figure  24  illustrates  this  process.  It  shows  an  example  of  a  raw  vorticity  field  (a)  and  the 
results  of  the  vortex  identification  scheme,  presented  as  an  idealized  vorticity  field  (b). 
Comparison  of  Fig.  24a  and  b  confirms  that  the  strongest  vortices  have  been  appropriately 
located  and  sized.  In  the  following  section  the  statistics  of  the  identified  vortices  will  be  used  to 
explore  the  variations  in  the  field  vorticity  distribution  and  its  connection  to  the  vertical  velocity 
fluctuations. 

4.3  Characteristics  of  individual  vortices 

The  convection  of  vortices  across  the  flow  field  will  induce  velocity  fluctuations  at  any 
fixed  point  in  the  field.  In  two  spatial  dimensions,  the  magnitude  of  the  induced  velocity,  |u|, 
from  a  point  vortex  is  given  by  the  Biot-Savart  law, 

\n\=r/(2xR),  (i) 

where  T  is  the  circulation  in  the  vortex,  and  R  is  the  distance  between  the  vortex  center  and  the 
point  of  interest.  Thus,  the  variables  that  control  the  induced  velocity  from  a  field  of  vortices 
include  the  size  and  strength  of  the  each  vortex  and  its  location  relative  to  the  point  of  interest. 


For  the  strong  vortices  identified  by  the  procedure  outlined  in  §4.2,  the  far-field 
circulation,  I\  and  core  diameter,  d  =  2  a,  are  shown  as  functions  of  downstream  distance  in  Fig. 
25  for  Geometry  I  at  Rec  =  4M  (a  panel)  and  Geometry  II  at  Rec  =  1.4M,  4M,  17M,  and  50M 
(panels  b-e).  On  both  Figs.  25  and  26,  a  vertical  dashed  line  indicates  the  location  of  the 
maximum  vertical  velocity  fluctuations  taken  from  Figs.  9a  and  10a. 

For  all  conditions,  the  vortices  are  largest  in  terms  of  both  circulation  and  diameter 
nearest  the  trailing  edge  and  reduce  in  strength  and  size  downstream.  Recall,  however,  that  the 
vorticity  within  1.00  <  x/C  <  1.02  may  be  due  both  the  swirling  and  shear  flow,  so  vortex 
statistics  there  may  be  biased  upward  by  the  shear  component.  If  the  calculations  leading  to  Fig. 
25  are  redone  after  excluding  identified  vortices  with  low  swirling  strength  within  their  cores,  the 
curves  are  shifted  downward  slightly  at  the  trailing  edge,  but  the  trends  remain  the  same. 
Beyond  these  common  trends,  two  conditions  in  Figure  25  are  distinctive.  The  flow  at  Rec  - 
1.4M  with  Geometry  II  (Fig.  25b)  shows  the  greatest  asymmetry  in  both  vortex  strength  and  size. 
This  asymmetry  is  due  to  the  transitional  state  of  the  pressure  side  boundary  layer  and  fully 
turbulent  state  of  the  suction  side  boundary  layer  at  this  Rec ■  All  other  conditions,  which 
represent  fully  turbulent  boundary  layers  from  both  sides  of  the  foil,  have  greater  symmetry. 
Interestingly  at  Rec  -  4M,  vortex  circulation  for  Geometry  II  exceeds  that  of  Geometry  I  by  a 
greater  factor  than  that  suggested  by  the  ratios  of  mean  circulation  fluxes  of  Table  3. 

Figure  26  presents  the  average  location,  and  vertical  separation,  of  the  centers  of  the 
suction  and  pressure  side  vortices  as  a  function  of  x/C  along  with  the  location  of  the  maximum 
vertical  velocity  fluctuations  taken  from  Figs.  9a  and  10a  (the  panel  layout  on  Fig.  26  is  identical 
to  that  of  Fig.  25).  The  vortex-center  curves  are  not  a  path  lines  for  individual  vortices  but  they 
do  indicate  likely  vortex  trajectories.  For  all  conditions,  these  vortex-center  curves  correspond 


well  to  the  location  of  peak  Reynolds  shear  stress  shown  in  Figs,  lib  and  12b.  Here,  both 
trailing  edge  geometries  at  Rec  =  4M  display  a  greater  downward  slope  in  the  vortex-center 
curves  (Fig.  26a  and  c)  than  is  seen  at  the  other  conditions. 

In  addition,  the  maximum  vertical  velocity  fluctuations  track  the  pressure-side  vortex- 
center  curves  for  all  conditions  except  those  at  Rec  ~  4M.  The  distance  from  the  vortex-center 
curves  to  the  location  of  maximum  velocity  fluctuations  gives  an  indication  of  the  relative 
contribution  of  the  suction  and  pressure  side  vortices  to  the  peak  velocity  fluctuations.  Thus,  it 
can  be  concluded  from  Fig.  26  that  the  pressure  side  vortices  dominate  the  velocity  field 
fluctuations  away  from  the  shedding  condition,  but  that  suction  and  pressure  side  vortices 
contribute  more  equally  to  the  velocity  fluctuations  when  shedding  occurs.  However,  closest  to 
the  foil,  1.00  <  x/C  <  1.01,  the  pressure  side  vortices  are  found  to  dominate  for  all  conditions 
tested.  Figure  26  also  shows  that  the  location  of  maximum  vertical  velocity  fluctuations  is 
typically  found  at  the  downstream  location  where  the  vertical  vortex  separation  (or  wake 
thickness)  is  smallest.  This  finding  is  consistent  with  the  interpretation  of  this  x-location  as  the 
vortex  street  formation  length,  where  the  vortices  become  well  staggered  and  move  into  closer 
vertical  proximity.  Thus,  in  the  remainder  of  this  paper,  the  downstream  distance  from  the 
trailing  edge  to  the  peak  in  vertical  velocity  fluctuations  is  referred  to  as  the  vortex  street 
formation  length. 

4.4  Vortex  street  organization 

From  simple  geometric  considerations,  the  vertical  velocity  fluctuations  in  the  middle  of 
the  near  wake  caused  by  suction  and  pressure  side  vortices  are  enhanced  when  these  vortices  are 
staggered  -  as  in  the  classical  Von  Karman  vortex  street  -  and  suppressed  when  they  are 
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vertically  aligned.  Thus,  fields  with  the  identical  individual  vortices  may  develop  very  different 
net  vertical  velocity  fluctuations  based  on  the  arrangement  or  organization  of  the  wake  vortices. 

Here  vortex  street  organization  was  determined  by  conditionally  averaging  instantaneous 
vorticity  fields  based  on  the  presence  or  absence  of  a  pressure  side  vortex  near  x/C  =  1.05.  First, 
instantaneous  vorticity  fields  with  no  pressure  side  vortex  centered  within  a  search  window  of 
1.045  <  x/C  <  1.055  were  discarded,  leaving  typically  half  of  the  original  field  measurements. 
Then,  the  remaining  fields  were  then  shifted  along  the  streamwise  axis  so  that  the  pressure  side 
vortex  of  greatest  peak  vorticity  within  the  search  window  was  centered  on  x/C  =  1.050.  Finally, 
these  shifted  fields  were  ensemble  averaged.  This  conditional  averaging  technique,  which  is 
meant  to  mimic  conventional  phase  averaging,  relies  solely  on  pressure  side  vortex 
identification,  so  any  suction  side  structure  found  is  indicative  of  vortex  street  organization. 

Figure  27  displays  the  conditionally-averaged  PIV-field  results  as  normalized  swirling 
strength  using  the  same  panel  layout  as  Fig.  25.  The  black  central  feature  in  each  panel  of  this 
figure  represents  the  conditionally-selected  (or  reference)  pressure-side  vortex.  The  result  in  Fig. 
27a  for  Geometry  I  at  Rec  =  4M  are  similar  to  those  at  the  other  Rec  with  this  trailing  edge 
geometry.  Based  on  a  comparison  Fig.  27b-e,  Geometry  II  at  Rec~  4M  (Fig.  27c)  produces  the 
most  clearly  structured  vortex  street  with  a  staggered  arrangement  of  vortices.  The  pressure-side 
vorticity  at  Rec=  1.4M  (Fig.  27b)  is  almost  as  structured  as  that  at  Rec=  4M,  but  the  suction  side 
vorticity  is  less  coherent  at  the  lower  Rec .  Results  at  both  of  the  higher  Rees  (Fig.  27d  and  e) 
indicate  little  geometric  coherence  between  pressure-  and  suction-side  vorticity.  The  results 
shown  in  Fig.  27  objectively  confirm  the  subjective  assessments  of  relative  organization  made 
from  the  individual  fields  of  Fig.  23.  In  addition,  the  vortex  street  at  Rec  =  4M  becomes  well 


organized  at  the  shedding  frequency  immediately  downstream  of  the  trailing  edge  so  that  it 

provides  the  most  constant  St  with  increasing  downstream  distance  (see  Figs.  20e  and  21e). 

( 

4.5  Reconstruction  of  vortex-induced  velocity  fluctuations 

To  determine  the  extent  to  which  the  large  vortices  determine  the  vertical  velocity 
fluctuations,  an  idealized  velocity  field  was  constructed  from  the  Biot-Savart  law  (1)  and  the 
idealized  vortices  identified  by  the  scheme  described  in  §4.2.  A  sample  vertical  velocity  field  is 
shown  on  Fig.  24c,  presented  as  the  instantaneous  vertical  velocity  minus  the  average.  Just 
below  it,  Fig.  24d,  is  the  reconstructed  idealized  vertical  velocity  for  the  same  field.  Working 
from  the  identified  vortices,  vertical  velocity  results  similar  to  Fig.  24d  may  be  generated 
separately  for  pressure  (vp)  and  suction  (v,)  side  vortices. 

Using  a  combined  Reynolds  and  vortex  decomposition,  v  =  v  +  v'=  vp  +  v„  the  mean- 
square  vertical  velocity  fluctuation  can  be  determined  in  terms  of  suction  and  pressure  side 
vortical  velocity  components. 

V7  -7)^+2^ -v?,\  (2) 

The  last  term  in  (2)  gives  the  contribution  from  correlation  of  the  suction  and  pressure  side 
vortices  and  quantifies  the  impact  of  vortex  street  stagger  or  alignment.  The  results  of  the 
decomposition  suggested  in  (2)  are  shown  in  Fig.  28  with  the  same  panel  layout  as  Fig.  25.  The 
separate  contributions  are  from  the  suction  side  vortices  (solid  line),  pressure  side  vortices 
(dashed  line),  and  the  con-elation  term  (dotted  line).  Total  vortex-induced  fluctuations  are  also 
shown  (dash-dot-dot  line).  The  symbols  denote  the  PIV-measured  vertical  velocity  fluctuations 
(from  Figs.  9b  and  10b).  A  comparison  of  the  total  idealized  fluctuations  to  the  actual 
fluctuations  for  Geometry  II  shows  the  relative  contribution  of  the  identified  vortices  to  the 
overall  fluctuation  levels.  The  peak  fluctuation  levels  are  well-matched,  and  the  trends  with 
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downstream  distance  are  similar.  This  confirms  that  the  induced  velocities  from  the  largest  - 
vortices  are  the  primary  factor  in  setting  near-wake  vertical-velocity  fluctuation  levels.  This  type  '/ 

r  J  s  -~r'"r 

■  ■  ■■  -o'-  •’  “*■ 

of  field  reconstruction  is  less  effective  with  Geometry  I  where  the  measured  vorticity  is  less  _  .  J_' 

coherent.  -  , . 

The  qualitative  trend  with  Re  and  trailing-edge  geometry  in  the  location  of  maximum 
fluctuations  (formation  length)  is  captured  by  the  reconstructed  fields.  The  formation  length  is 
shortest  for  Geometry  II  at  Rec=  4M  (c)  and  increases  for  the  other  conditions.  In  all  cases,  the 
agreement  between  the  actual  and  reconstructed  fluctuations  is  improved  by  lowering  the  vortex 
peak  threshold,  thereby  increasing  the  fraction  of  field  vorticity  included  in  the  identified 
vortices.  A  lower  vorticity  threshold  is  of  particular  benefit  for  the  Rec=  1.4M  case  which  has 
the  greatest  asymmetry  between  the  suction  and  pressure  side  vorticity.  Here,  as  the  threshold  is 
reduced,  more  suction  side  vorticity  is  captured  with  little  change  to  the  identified  pressure  side 
vorticity.  The  net  effect  of  the  added  vorticity  is  an  increase  in  the  formation  length  in  the 
reconstructed  fields. 

The  vortex  decomposition  shown  on  Fig.  28  also  confirms  that  the  pressure  side  vortices 
are  the  main  contributors  at  the  vertical  velocity  fluctuation  maximum.  The  x/C  location  of  peak 
vertical  velocity  fluctuations  coincides  in  all  cases  with  the  peak  in  the  contribution  of  the 
pressure  side  vortices,  and  in  no  cases  corresponds  to  the  peak  in  suction  side  contributions. 

However,  for  Rec  -  4M,  the  suction  side  contributions  and  their  correlation  with  the  pressure 
vortices  side  do  play-  a  secondary  but  significant  role.  The  highest  total  fluctuations  occur  at  Rec 
=  4M  with  Geometry  II,  even  though  the  highest  pressure  side  fluctuations  occur  at  Rec-  1 .4M. 
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5.  Discussion:  Prediction  and  Scaling  of  Vortex  shedding 

The  prior  analysis  of  PIV  fields  shows  that  vertical  velocity  fluctuations  near  the  trailing 
edge  are  primarily  due  to  the  pressure  side  vortices.  The  distribution  of  the  pressure  side  vortices 
appears  connected  to  the  suction  side  shear  layer  properties  and  the  wake  thickness.  In  this 
section  this  connection  is  put  on  a  quantitative  footing.  The  vertical  velocity  induced  at  the 
trailing  edge  by  the  suction  side  vortices  is  presumed  to  drive  perturbations  of  the  unstable 
pressure  side  shear  layer  and  influence  its  roll-up  on  the  wake-thickness  scale.  The  pertinent 
properties  of  the  suction  side  vortices  are  tied  to  parameters  from  the  suction-side  mean  flow. 
The  physical  model  and  the  proposed  scaling  law  are  presented  in  the  next  subsection  (§5.1), 
which  is  followed  by  short  discussions  of  the  roles  of  suction  (§5.2)  and  pressure  (§5.3)  side 
mean  shear  rate,  predictive  application  of  the  proposed  scaling  law  (§5.4),  and  the  potential  for 
wake  instability  concepts  to  explain  the  findings  of  this  study  (§5.5). 

5.1  Shedding  strength  and  suction  side  mean  shear 

The  measurements  reported  here  are  consistent  with  a  physical  model  of  vortex  shedding 
that  correlates  the  Biot-Savart  induced  vertical  velocity  at  the  trailing  edge  and  the  roll-up  into 
vortices  of  the  pressure  side  shear  layer.  Consider  the  idealized  vortex  street  depicted  in  Fig.  29. 
The  suction  side  flow  is  best  expressed  in  streamline  coordinates,  which  are  well  approximated 
by  the  h-t  Cartesian  frame  (Fig.  1)  aligned  to  the  suction  side  surface  tangent  at  x/C  =  0.93.  The 
corresponding  streamline  coordinates  for  the  pressure  side  are  simply  identical  to  the  tunnel  x-y 
coordinate  frame.  On  Fig.  29,  the  hydrofoil  is  depicted  with  suction  and  pressure  side  external 

mean  velocities  U“  and  Ups .  The  center  lines  of  the  suction  and  pressure  side  shear  layers  are 
indicated  with  dashed  lines  and  the  characteristic  convection  velocities  of  these  shear  layers  are 
denoted  Usss and  Ups ,  respectively.  The  vortices  of  alternating  sign,  which  form  the  vortex 
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street  are  shown  as  gray  circles,  and  the  vertical  distance  from  the  suction  side  vortices  to  the  ; . .  ’ 


pressure  side  shear  layer  at  the  trailing  edge  is  denoted  Ar  The  instant  in  the  shedding  cycle 
shown  is  the  time  of  pressure  side  vortex  roll-up  at  the  trailing  edge.  The  suction  side  vortex^ 
with  contribution,  v",  to  the  induced  vertical  velocity  at  the  trailing  edge  is  acting  with  clock- 


wise  circulation,  T,  at  distance  R.  This  vortex  is  presumed  to  dominate  the  flow  fluctuations  at 


the  trailing  edge  at  the  instant  shown.  The  perturbation  from  this  vortex  on  the  pressure  side 
shear  layer  at  the  trailing  edge  is  represented  by  the  ratio  of  the  induced  vertical  velocity  to  the 
streamwise  velocity  of  the  shear  layer.  From  (1),  this  ratio  can  be  approximated  as: 


-i - - — cose  V) 

Usps  2 nRUps 

The  circulation  (strength)  of  the  suction  side  vortex  can  be  approximated  by  the  mean  ^ 
suction-side  circulation  flux  passing  the  trailing  edge  in  one  shedding  cycle,  integrated  over  a  «•- 

vertical  distance  Ay  and  centered  on  the  suction  side  shear  layer:  T  « [{dU,/ dh)s  U5*Ay  ]//p  •  The 


appropriate  mean  velocity  gradient,  ( dUjdh )“,  and  convection  velocity,  are  those  values 
characteristic  of  the  center  of  the  mean  suction  side  shear  layer.  In  this  approximation,  the 
correlation  between  instantaneous  vorticity  fluctuations  and  velocity  fluctuations  in  the  shear 
layer  are  assumed  negligible  in  comparison  to  the  product  of  their  means.  This  circulation  flux 
does  not  reflect  the  total  flux  available  from  the  attached  suction  side  boundary  layer,  but  instead 
the  lesser  amount  residing  near  the  shear  layer  centerline  at  the  trailing  edge.  For  the  flows  of 
this  study,  R  ~  Ay  and  0-^/4  for  the  moderately  and  weakly  shedding  conditions,  though 
occasional  instantaneous  fields  with  zero  6  were  observed  for  the  strongest  shedding  conditions. 
Overall,  variation  in  6  between  conditions  was  weak,  so  for  simplicity,  6  can  be  treated  as  a 
constant  and  (3)  becomes: 
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(4) 


y"  o:  wacf  u: 

Ur  2 nfp  UT 

When  the  shear  layer  velocity  profiles  are  of  similar  shape,  the  external  and  shear  layer 
velocities  should  be  proportional,  U?JU?  ~U?  jlJ?  .  This,  and  the  further  simplification, 

U™  «  Ups  «  U0 ,  allows  (4)  to  be  adjusted  and  cast  in  terms  of  St: 

va  1  (dUjdhf  (5) 

Ua  2nSt  U0/Ay  ’ 

where  2nSt  -1  for  the  flows  of  this  study.  The  scaling  suggested  by  (5)  relates  of  v“/U0  to  a 
dimensionless  suction-side  shear  rate,  (dUt/dh)a /(U0/Ay).  This  parameter  can  be  calculated 

from  the  mean  velocity  profiles  at  the  foil’s  trailing  edge  (see  Part  I,  Fig.  18).  Here  ( dUjdh )" 
was  calculated  as  the  average  velocity  gradient  from  3.2  mm  below  to  3.2  mm  above  the  vertical 
coordinate  of  peak  streamwise  velocity  fluctuations  along  the  line  A-B  shown  in  Fig.  29.  Figure 
30  shows  the  resulting  values  of  (dUjdh)ss /(uj A™m)  vs.  Rec  for  Geometry  I  (a)  and  Geometry 
II  (b).  For  both  geometries,  the  dimensionless  shear  rate  peaks  at  Rec  =  4M,  the  maximum 
shedding  strength  condition. 

The  form  of  right  side  of  (5)  is  similar  to  the  shedding-strength  parameter  b/d'  suggested 
in  Blake  (1986),  where  b  is  the  vertical  thickness  of  the  body  and  <5*  is  the  average  of  the  suction 
and  pressure  side  boundary-layer  displacement  thicknesses.  Based  on  a  review  of  experimental 
data  up  through  the  early  1980’s,  strong  vortex  shedding  was  found  to  occur  when  b/ <5*  s  0.3  for 
turbulent  flow  airfoils  (Rec  near  or  above  approximately  2M).  The  present  results  provide  a 
refinement  of  these  ideas  for  higher  Rec- 

As  the  flow  evolves  beyond  the  instant  shown  on  Fig.  28,  the  velocity  ratio  v*" /U 0  sets 

the  input  perturbation  level  to  the  Kelvin-Helmholtz  instability  of  the  pressure  side  shear  layer. 
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Such  instabilities  have  been  found  active  in  other  separated  turbulent  flows  (Song  and  Eaton, 
2002).  Therefore,  the  strength  of  large  near-wake  vortices  formed  at  the  shedding  frequency 
should  increase  when  this  ratio  increases.  Here,  the  strength  of  vortex  shedding  is  quantified 
from  the  measured  peaks  in  the  LDV-acquired  near- wake  velocity  spectra  at  x/O 1.009  (see 
Figs.  16,  17  and  19-21).  The  advantage  of  this  choice  for  shedding  strength  is  that  it  allows 
comparisons  with  the  results  of  Blake  (1984)  for  minimally  lifting  struts. 

The  parameters  of  the  vertical-velocity  spectral  shedding  peaks  are  plotted  against  the 
dimensionless  shear-rate  parameter  of  (5)  with  2xSt  =1  in  Fig.  31.  Here,  the  wake  thickness, 
(indicated  on  Fig.  29),  is  the  measured  vertical  distance  between  the  suction  side  Reynolds  stress 
peak  at  x/C=  1.009  and  the  pressure  side  Reynolds  stress  peak  at  x/C  =  1.002.  For  this  definition, 
the  wake  thicknesses  for  Geometry  I  are  0.0085C  at  Rec  =  3M,  0.0080C  at  Rec  =  4M,  and 
0.0090C  at  the  remaining  Rec.  The  wake  thicknesses  for  Geometry  II  are  0.0095C  at  all  Rec. 
The  separate  panels  of  Fig.  31  show  normalized  peak  height  (a),  normalized  peak  width  (b),  and 
the  fraction  of  total  spectral  area  found  under  the  peak  (c),  and  the  total  spectral  area  (d).  Results 
from  Geometry  I  and  II  are  shown  as  black  and  open  circles,  respectively,  with  one  exception. 
The  datum  on  each  panel  from  Geometry  I  at  Rec  —  1  .4M  is  depicted  as  a  black  square  to 
represent  its  disparity  of  pressure-side  boundary  separation  upstream  of  the  trailing  edge.  The 
gray  filled  triangles  in  Fig.  31a  are  data  from  three  minimally  lifting  struts  (Blake  1984)  at  Rec~ 
2M,  two  having  rounded  suction  side  trailing  edge  bevels  of  25°  and  45°,  and  a  third  having  a 
sharp  trailing  edge  bevel  of  25°.  In  these  strut  studies,  the  boundary  layers  were  tripped  near  the 
leading  edge  and  spectral  measurements  were  made  with  hot  wire  anemometry  approximately 
2.5A y  downstream  of  the  various  trailing  edges  (compared  to  ~1.0  ^  f°r  this  study).  A  zero 
shedding  strength  is  presumed  for  the  two  25°  edges  based  on  the  absence  of  peaks  in  their 
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reported  surface  pressure  spectra.  The  error  bars  for  the  25°-edge  data  points  (plotted  at  the 
bottom  of  Fig.  31)  are  omitted  for  graphical  clarity,  but  are  similar  to  that  shown  for  the  45° 
edge. 

Figure  31a  and  c  show  that  vortex  shedding  strength  measurements  collapse  well  when 
plotted  against  the  dimensionless  suction-side  shear  rate.  Data  from  Geometry  I  and  II,  as  well 
as  minimally-lifting  struts,  are  consistent  with  a  single  trend.  The  rapid  increase  in  the  shedding 
strength  at  a  dimensionless  shear  rate  of  unity  is  due  both  to  an  increase  in  vortex  strength  and  to 
a  reduction  in  the  formation  length  itself.  These  combined  effects  are  nonlinear,  and  no  physical 
basis  is  proposed  for  the  functional  forms  of  the  curve  fits  given  in  the  caption  of  Fig.  3 1 . 

The  proposed  role  of  the  induced  vertical  velocity  at  the  trailing  edge  is  consistent  with 
the  established  characteristics  of  vortex  shedding.  First,  trailing  edge  motion,  which  increases 
the  trailing  edge  vertical  velocity  ratio,  has  been  shown  to  enhance  vortex  shedding  and  decrease 
the  formation  length  in  wind  tunnel  tests  of  a  strut  with  a  beveled  trailing  edge  at  Rec  ~  1M 
(Greenway  and  Wood  1973).  Second,  the  use  of  splitter  plates  (Roshko  1955,  Morkovin  1964, 
Bearman  1965)  and  base  bleed  (Bearman  1967,  Wood  1967)  are  known  to  reduce  vortex 
shedding,  and  modification  of  the  vertical  velocity  ratio  at  the  trailing  edge  may  be  the 
mechanism  behind  these  effects.  Finally,  Kuethe  (1972)  showed  for  both  a  flat  plate  and  a  lifting 
surface  at  Rec  approaching  1M,  that  the  diversion  of  a  portion  of  the  boundary  layer  vorticity 
into  streamwise  vortices  (by  a  wave-like  static  surface  feature)  suppressed  vortex  street 
development  in  the  near  wake.  Such  streamwise  vortices  may  disrupt  vortex  shedding  by 
destroying  spanwise  coherence  of  the  vertical  velocity  ratio  at  the  trailing  edge. 


5.2  Pressure  side  mean  shear  and  vortex  shedding 

It  is  surprising  that  the  scaling  embodied  in  (5)  and  illustrated  in  Fig.  31a  and  c  lacks 
variables  connected  to  the  pressure  side  shear  layer  since  trailing-edge  vortex  shedding  might 
reasonably  be  expected  to  depend,  in  general,  on  the  properties  of  both  shear  layers.  One 
possible  explanation  involves  the  more  downstream  location  of  pressure  side  separation  in  the 
flows  studied  here.  The  roll-up  of  the  suction  side  shear  layer  initiates  well  upstream  of  trailing 
edge.  Thus,  The  more  fully  formed  suction  side  vortices  may  dominate  the  interaction  between 
the  two  shear  layers  at  the  trailing  edge.  A  second  explanation  follows  from  the  lack  of  variation 
in  the  foil’s  pressure-side  shear  layers  in  this  study;  the  pressure  side  circulation  flux  only  varies 
by  14%  across  all  conditions,  and  only  3%  across  the  conditions  with  Rec>  2M.  If  the  effect  of 
the  pressure  side  variables  were  weak  in  comparison  to  those  of  the  suction  side  for  the 
geometries  studied,  small  variations  caused  by  pressure  side  characteristic  changes  might  go 
undetected.  Only  early  pressure  side  separation  at  Rec  =  1.4M  appears  to  have  influence  on  the 
shedding  characteristics  shown  on  Fig.  31.  However,  pressure  side  variables  are  expected  to 
increase  in  importance  for  more  symmetric  foils,  or  when  suction  and  pressure  side  separation 
points  and  initial  shear  layer  conditions  are  more  similar. 

5.3  Suction-side  mean  shear  and  wake  frequency  selection 

The  scaling  presented  by  (5)  suggests  that  near  the  trailing  edge,  there  exists  a  negotiation 
between  the  suction  and  pressure  side  shear  layers  to  determine  whether  or  not  the  pressure  side 
shear  layer  will  respond  at  a  forcing  frequency  set  by  the  suction  side.  When  the  both  sides 
agree  on  the  frequency,  there  is  strong  shedding,  such  as  Rec  =  4M  with  Geometry  II.  When  the 
two  shear  layers  initially  disagree,  vortex  street  formation  may  be  delayed  or  even  suppressed. 
When  strong  near-wake  fluctuations  are  absent,  the  near-wake  shear  layers  tend  to  roll  up  on 
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length  scales  set  by  the  separating  boundary  layers  (Morris  and  Foss  2001).  Such  small-scale 
roll-up  is  evident  in  a  flow  visualization  of  Prasad  and  Williamson  (1997)  for  a  cylinder  with 
vortex  shedding  at  Red=  104,  on  which  the  shear  layers  roll  up  first  on  the  smaller  boundary  layer 
scale,  and  then  further  downstream  location  on  the  larger  wake  scale.  This  phenomenon  is  also 
evident  for  Geometry  II  flow  at  Rec  =  50M  where  the  preferential  frequency  in  the  initial 
pressure  side  shear  layer  roll  up  is  higher  than  the  expected  vortex-shedding  Strouhal  frequency. 
However,  the  suction  and  pressure  side  shear  layers  remain  so  ill-matched  that  significant  wake- 
scale  vortex  organization  does  not  occur  upstream  of  x/C  =  1.10.  These  observations  suggest 
that  near-wake  symmetry  is  important  in  setting  the  vortex  shedding  strength,  and  the  strongest 
shedding  case  found  in  these  studies  does  occur  when  the  mean  velocity  profiles  at  x/C  —  1 .002 
are  most  symmetric.  However,  both  shear  layers  are  in  general  turbulent  so  the  suction  side  layer 
must  have  enough  time  to  form  a  sufficiently  strong  vortex  to  convince  the  pressure  side  shear 
layer  to  cooperate.  Strong  vortex  shedding  cannot  occur  at  faster  rate  than  this  because  the 
suction  side  vortices  cannot  accumulate  enough  vorticity  to  be  strong  enough.  Shedding  does 
not  occur  at  a  much  slower  rate  either  because  the  near  wake  vortices  must  be  near  enough  to 
each  other  to  interact,  but  they  accelerate  rapidly  away  from  the  trailing  edge  as  they  leave  the 
separated  flow  region  and  cannot  interact  if  formed  too  infrequently.  Thus,  a  fluid-dynamic 
compromise  between  suction-side  shear-layer  influence  and  pressure-side  shear-layer  response 
appears  to  determine  the  dominant  roll-up  frequency  at  the  trailing  edge,  St  ~  1/2k. 

5.4  Predictive  application  of  the  shedding-shear  correlation 

Of  the  forms  investigated,  the  dimensionless  shear  rate  in  (5)  provides  the  best 
compromise  between  collapse  of  the  vortex  shedding  data  and  potential  for  predicting  vortex 
shedding  at  untested  flow  conditions. 
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From  Fig.  31a  or  c,  vortex  shedding  is  predicted  to  occur  when: 

{dUjdhf  __  Q  . 

U.IK 


(6) 


This  equation  can  be  used  to  form  predictions  if  the  mean  suction  side  shear  rate  and  wake 
thickness  can  be  estimated  from  the  flow  speed  and  the  foil’s  size  and  trailing  edge  geometry. 
Although  potential  flow  and  boundary  layer  integral  calculations,  or  even  Reynolds-averaged 
Navier-Stokes  (RANS)  computational  fluid  dynamics  (CFD)  calculations  may  be  nearly  up  to 
this  task,  the  main  problem  with  using  (6)  in  a  predictive  fashion  is  that  ( dUjdh )  is  influenced 
by  the  presence  or  absence  of  vortex  shedding.  However,  vortex  shedding  tends  to  decrease 


(dUt/dh)ss  so  it  will  tend  to  be  overpredicted  by  any  calculation  scheme  that  cannot  itself  predict 
vortex  shedding.  Thus,  any  flow  solution  -  perhaps  obtained  from  a  RANS  model  -  that 
potentially  captures  the  relevant  physics  of  the  weakly  shedding  cases  might  be  used  with  (6) 
and/or  Fig.  31  to  estimate  an  upper  bound  on  the  strength  of  vortex  shedding  in  the  actual  flow. 


5.5  Wake  instability 

Stability  theory  potentially  offers  an  alternative  means  to  explain  the  importance  of  the 
suction  side  mean  shear  and  the  observed  differences  in  vortex  street  formation  length.  The 
strongest  shedding  case  is  also  the  condition  with  wake-scale  vortex  organization  nearest  to  the 
trailing  edge.  This  may  reflect  a  state  of  absolute  instability  in  the  underlying  flow,  where  wake 
scale  disturbances  are  able  to  travel  upstream  and  attach  themselves  to  the  foil’s  trailing  edge. 
However,  a  stability  analysis  starting  from  the  trailing-edge  mean-velocity  profiles  was  beyond 
the  scope  of  this  work,  but  it  is  a  possible  next  step  toward  elucidating  the  role  of  the  suction- 
and  pressure-side  mean  shear  in  vortex  shedding. 
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6.  Summary  and  Conclusions 

An  extensive  experimental  investigation  of  the  dynamic  flow  in  the  near  wake  of  a  twcH 
dimensional  hydrofoil  at  low  Mach  number  and  high  Reynolds  number  is  presented  in  this  paper. 

The  LDV,  PIV  and  surface  dynamic  pressure  measurements  reported  here  were  focused  on 
understanding  the  unsteady  flow  in  the  foil’s  trailing  edge  region  where  near  wake  vortex 
shedding  occurs.  Results  for  chord  based  Reynolds  numbers  from  1 .4  million  to  50  million  are 
reported  for  two  trailing  edge  shapes  have  44°  and  56°  apex  angles.  This  effort  has  lead  to  four 
main  conclusions. 

/ 

First  of  all,  the  profiles  of  turbulent  fluctuations  through  the  near  wake  are  Reynolds 
number  dependent  because  of  the  varying  strength  of  structured  near-wake  vortex  shedding,  even 

at  the  Reynolds  numbers  of  this  study  that  reach  those  of  full-scale  aero-  and  hydrodynamic  . .  v 

lifting  surface  applications.  The  relative  strength  of  the  vortex  shedding  was  found  to  depend  on  * 
trailing  edge  geometry  as  well.  This  conclusion  is  drawn  from  the  statistical,  spectral,  and  '  ; 

instantaneous  field  results  shown  on  Figs.  2-5,  9-18, 22-24  and  27.  '  * 

Second,  strong  vortices  from  both  the  suction  and  pressure  side  shear  layers  are  organized 
into  a  staggered  vortex  street  that  dominates  the  near  wake  fluctuations  when  vortex  shedding 
occurs.  This  conclusion  is  based  on  the  conditionally-averaged  vorticity  field  results  and  the 
relative  success  of  using  a  finite  and  countable  number  of  idealized  vortices  to  reconstruct  the 
measured  vertical  velocity  fluctuations  in  the  foil’s  near  wake  (see  Figs.  24,  27  and  28). 

* 

However,  when  shedding  is  weak  or  absent,  pressure  side  vortices  dominate  the  fluctuations  in 
the  foil’s  near  wake  and  there  is  little  or  no  prevalent  geometrical  structure  of  wake  vortices. 

Third,  the  amplitude  and  frequency  content  of  vertical  velocity  fluctuations  at  the  foil’s 
trailing  edge  appear  to  control  the  strength  of  near-wake  vortex  shedding.  This  finding  is 
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consistent  with  the  near-wake  dynamics  in  the  current  measurements  and  with  past  shedding 
enhancement  and  suppression  techniques. 

And  finally,  the  proposed  scaling  for  vortex  shedding  strength  based  on  the  mean  suction- 
side  shear  rate,  the  upstream  velocity,  and  the  wake  thickness  is  successful  throughout  the  Rec 
parameter  range  of  this  study  for  both  trailing  edge  geometries,  and  is  consistent  with  prior 
studies  at  lower  Reynolds  number  on  minimally  lifting  struts  with  asymmetric  trailing  edges. 
Although  this  shedding  strength  scaling  relies  on  measured  parameters,  it  may  be  possible  to  use 
it  with  mean-flow  CFD  calculations  to  predict  an  upper  bound  on  vortex  shedding  strength. 

In  summary,  the  experimental  results  presented  in  this  paper,  together  with  those 
presented  in  Part  I,  provide  new  insights  into  the  relationships  between  the  various  steady  and 
unsteady  flow  phenomena  on  a  two-dimensional  but  application-relevant  lifting  surface  at  high 
Reynolds  number.  These  experimental  results  also  provide  a  unique  dataset  for  the  development 
and  validation  of  scaling  laws  and  computational  models  for  wall-bounded  high-Reynolds 
number  flows 
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Tables 
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Table  1.  Line  and  Symbol  Types  for  Flow  Speeds  and  Rec  Values 


Uo 

m/s 

Rec 

(millions) 

Symbol/Line-type 

0.5 

1.4 

- ■ - 

1.5 

4 

▲ - 

3.0 

8 

- y - 

6.0 

17 

. ►- - 

12.0 

33 

•  •••••••  ••••••••• 

18.3 

50 

_  •  •  —  •  •  -  •  •  -  • 

Table  2 .y/C  Spectral  Measurement  Locations  of  the  LDV  Focal  Volume 


Rec 

x/C=  1.0094 

x/C  =  1.0703 

(millions) 

Geometry  I 

Geometry  II 

Geometry  I 

Geometry  II 

1.4 

0 

0.0023 

-0.0046 

-0.0094 

2 

0 

0.0023 

-0.0023 

-0.0046 

3 

0.0023 

0.0023 

-0.0023 

-0.0023 

4 

0.0023 

0.0023 

0 

-0.0023 

8 

0.0023 

0.0023 

0 

0 

17 

0.0023 

0.0023 

0 

0 

50 

0.0023 

0.0023 

0 

0 

Table  3.  Boundary  Layer  Circulation  Fluxes  at  Separation 

Rec 

Suction  Side: 

f”/r* 

Pressure  side: 

ryf* 

(millions) 

Geometry  I 

Geometry  II 

Geometry  I 

Geometry  II 

1.4 

1.46 

— 

1.01 

1.08 

4 

1.45 

— 

0.97 

1.00 

17 

1.35 

— 

1.03 

1.02 

50 

1.36 

- 

1.05 

1.05 

t'  =  [tpslul) 

\  '  0 1  Re  -4M 

0.400 

» 
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Figure  Captions 

Figure  1.  (a)The  cross  section  of  the  Geometry  I  hydrofoil  with  chord  and  max  thickness 

indicated.  The  chord  length  given  is  the  idealized  value,  measured  to  the  vertex  of  the 
trailing  edge  apex  angle.  This  point  is  defined  as  (x/C,  y/C)  =  (1, 0)  in  the  tunnel 
coordinate  system:  x-axis  parallel  to  streamwise  axis  of  the  tunnel,  and  y -axis  vertical. 
Also  depicted  is  the  hydrofoil  coordinate  system,  for  which  the  t-axis  is  taken  as  the 
local  surface  tangent  and  the  /i-axis  is  the  local  surface  normal,  (b)  Trailing  edge 
details  showing  Geometry  I  and  Geometry  II.  The  0.4-mm  radius,  applied  to  the 
trailing  edge  tip,  is  not  depicted.  The  dashed  line  indicates  the  direction  of  the  surface 
tangent  at  x/C=0.930. 

Figure  2.  LDV-measured  boundary  layer  characteristics  at  x/C  =  0.93  and  Rec  =  1 -4M,  4M, 
17M,  and  50M  on  the  (a)  suction  and  (b)  pressure  side  surfaces  of  the  Geometry  I 
hydrofoil.  Panels  show  the  normalized  mean  velocity  components,  velocity 
fluctuation  variances,  and  Reynolds  shear  stress,  all  in  the  tunnel  coordinate  frame. 
Symbols  are  defined  in  Table  1.  No  data  are  available  at  the  lower  two  Rec  on  the 
pressure  side. 

Figure  3.  LDV-measured  separating  boundary  layer  and  near- wake  fluctuations  for  Geometry  I 
at  Rec  =  8M  and  50M,  presented  as  normalized  (a)  streamwise  fluctuations  u’2/ul ,  (b) 

vertical  fluctuations  v'2jul ,  and  (c)  Reynolds  shear  stress,  u'v'/ul .  The  trailing  edge 
geometry  is  depicted  in  the  left  half  of  each  the  frame.  Vertical  gray  lines  are  shown 
at  the  x/C  coordinate  of  the  measurements  and  provide  the  vertical  axis  for  plotted 
profile.  The  scale  used  to  set  the  horizontal  extent  of  each  profile  is  provided  in  the 
bottom-left  comer  of  the  frame.  Symbols  are  defined  in  Table  1.  The  velocity 
profiles  are  terminated  at  the  locations  of  99%  boundary  layer  or  wake  thickness(es). 
Results  from  Rec  =  17M  and  33M  are  omitted  for  clarity  but  fall  monotonically 
between  the  data  presented. 

Figure  4.  PIV-measured  separating  boundary  layer  and  near-wake  fluctuations  for  Geometry  I 
at  Rec  =  1 -4M,  4M  and  50M,  presented  in  format  similar  to  Fig.  3.  Data  are  presented 
at  x/C= 1.002  in  lieu  of  x/C=1.000.  Line  types  are  defined  in  Table  1.  Data  at 
Rec=1.4M  were  not  available  at  the  downstream  station. 

Figure  5.  PIV-measured  separating  boundary  layer  and  near-wake  fluctuations  for  Geometry  II 
Rec  =  1 .4M,  4M  and  50M,  presented  in  the  format  of  Fig.  3 .  Data  at  Rec=  50M  was 
not  available  at  the  three  upstream  stations. 

Figure  6  Convergence  history  for  the  PIV  measurements  at  (x/C,  y/C)  =  (1.0094,  0.0080)  in  the 
near  wake  with  Geometry  I  at  Rec  =  50M.  Here,  N  is  the  number  images  used  to 
compute  the  normalized  mean  velocity  components  (a),  and  velocity  variances  and 
and  Reynolds  shear  stress  (b). 
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Figure  7.  Comparison  of  velocity  variances  and  Reynolds  shear  stress  at  x/C— 1.0094  with 
Geometry  I  at  Rec  =  50M  between  LD V-measured  results  (•)  and  PIV-measured 
results  ( - ). 


Figure  8. 


Ratio  of  LD  V-measured  to  PIV-measured  velocity  variances  and  Reynolds  shear 
stresses  measured  by  at  x/C=  1.0094  with  Geometry  I  at  Rec= 17M,  33M,  and  50M. 
Data  at  Rec=  8M  are  unavailable.  Streamwise  velocity  variance  ratio  ■  ' 


U.'1  LDV  / ua  piv  ;  vertical  velocity  variance  ratio  *  v'2  ldv / v?2  piv  ;  Reynolds 

shear  stress  ratio  ®  u'v'ldv  j u'v'piv  . 


*.  : '  % 


Figure  9.  Loci  of  maximum  vertical  velocity  variance  at  Rec  -  1.4M,  4M,  17M,  and  50M, 
taken  from  Fig.  4  for  Geometry  I  (a),  and  vertical  velocity  variance  values  at  these 
loci  (b). 


Figure  10.  Same  as  Fig.  9  except  these  data  are  for  Geometry  II. 

Figure  11.  Same  as  Fig.  9  except  these  data  are  for  the  normalized  Reynolds  shear  stress.  Data  at 
Rec=  17M  are  omitted  for  clarity,  but  are  similar  to  Rec-  50M. 


Figure  12.  Same  as  Fig  1 1  except  these  data  are  for  Geometry  II. 

Figure  13.  PIV-measured  streamwise  velocity  variance  in  the  near  wake  for  Geometry  I  (filled 
symbols)  and  Geometry  II  (open  symbols)  at  varying  x! C  for  Rec~  1 -4M  (a),  4M  (b), 
17M  (c),  and  50M  (d).  The  data  are  plotted  in  wake  coordinates  and  scaled  by  the 
time-averaged  velocity  deficit,  ALrmax.  The  x/C  coordinate  is  indicated  in  each  panel. 
The  solid  line  in  the  far  right  panel  gives  the  self-preserving  wake  profile  of  a 
symmetric  airfoil  at  Rec  ~  10  (Wygnanski  et  al.,  1986). 

Figure  14.  Same  as  Fig.  13  except  these  data  are  for  the  normalized  vertical  velocity  variance. 

Figure  15.  Same  as  Fig.  13  except  these  data  are  for  the  normalized  Reynolds  shear  stress. 

Figure  16.  Near- Wake  LD  V-measured  temporal  power  spectra  of  vertical  velocity  fluctuations  at 
x/C  =  1 .0 1  at  Rec=  1 .4M,  4M,  1 7M  and  50M  for  Geometry  I  (a)  and  Geometry  I  (b). 
These  data  were  taken  at  the  location  of  peak  vertical  velocity  fluctuations;  Table  2 
lists  their  y/C  locations.  Numerals  indicate  Rec  in  millions.  The  line  to  the  right  of 
the  data  has  a  slope  of-5/3. 

Figure  17.  Same  as  Fig.  16  except  these  data  are  from  x/C=  1.07. 

Figure  18.  PIV-measured  spatial  power  spectra  of  vertical  velocity  fluctuations  along  the  line 
segment  defined  by  (x/C,  y/C)  =  (1.000,  0.009)  and  (1.100,  -0.005)  at  Rec=  1.4M, 

4M,  17M  and  50M  for  Geometry  I  (a)  and  Geometry  I  (b).  This  line  falls  near  the 
locus  of  peak  vertical  velocity  fluctuations  for  1.03<x/C  <1.10.  Each  spectrum  is 
computed  from  between  500  and  1000  vector  fields.  Numerals  indicate  Rec  in 
millions.  The  line  to  the  right  of  the  data  has  a  slope  of-5/3. 
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Figure  19.  Illustration  of  the  parameters  used  to  quantify  the  LDV-measured  vertical  velocity 
spectral  peaks  in  the  near  wake.  Letters  correspond  to  the  panels  in  Figs.  20  and  21. 

Figure  20.  Characteristics  of  LDV-measured  vertical- velocity  near- wake  spectral  peaks  at  x/C  = 
1 .0 1  vs.  Rec  in  millions:  Geometry  I  ,  and  Geometry  II  -O- .  Parameters  are 
illustrated  in  Fig.  19.  Panels  show:  the  normalized  vertical  velocity  variance,  the 
height  of  the  spectral  peaks  above  the  spectral  background  (b),  the  spectral 
background  value  (c),  the  spectral  peak  width  (d),  and  the  frequency  of  the  spectral 
peak  (e)  presented  as  a  St  multiplied  by  2k. 

Figure  21.  Same  as  Fig.  20  except  these  data  are  from  x/C=  1.07. 

Figure  22.  The  power  spectra  of  the  pressure-side  surface  pressure  fluctuations  at  x/C  =  0.99  at 
various  Rec  for  Geometry  I  (a)  and  Geometry  II  (b).  Numerals  indicate  Rec  in 
millions.  The  solid  line  below  the  data  has  a  slope  of-1. 

Figure  23.  Instantaneous  PIV  fields  of  normalized  swirl,  n°m  /U0  for  Geometry  I  at  Rec=  4M 

(a),  and  Geometry  II  at  Rec  =  1 .4M  (b),  4M  (c),  17M  (d),  and  50M  (e). 

Figure  24.  Actual  and  idealized  instantaneous  velocity  and  vorticity  fields  for  Geometry  I  at  Rec 
=  4M.  Panels  present  grayscale  images  of  measured  normalized  vorticity  a>An°m  /Ua 

(a),  vorticity  of  the  idealized  vortices  (b),  measured  normalized  vertical  velocity 
fluctuation,  (v  -  V)/U0  (c),  and  normalized  vertical  velocity,  v/U0 ,  computed  from  the 

idealized  vortex  field  of  (b). 

Figure  25.  Vortex  circulation  and  core  size  of  suction  side  vortices  ( '  )  and  pressure  side 

vortices  ( - )  vs.  downstream  distance  in  the  near  wake.  The  panels  follow  the 

same  arrangement  as  Fig.  23.  The  thicker  lines  represent  normalized  core  circulation, 
read  from  axes  on  the  left.  Thinner  lines  represent  core  diameter,  read  from  the  axes 
on  the  right.  The  vertical  dashed  line  in  the  left  part  of  each  frame  indicates  the  x/C 
location  of  the  peak  vertical  velocity  fluctuations  from  Figs.  9  and  10. 

Figure  26.  Vortex  center  locations  vs.  downstream  distance  in  the  near  wake;  suction  side 

vortices  - ;  pressure  side  vortices - ;  vertical  distance  between  suction 

side  and  pressure  side  vortices  .  (illustrated  in  panel  a);  and  vertical  location 

of  maximum  vertical  velocity  fluctuations  at  the  given  x/C  .  The  panels 

follow  the  same  arrangement  as  Fig.  23 

Figure  27.  Conditionally-averaged  PIV  fields  of  normalized  mean  swirl,  rendered  as 

grayscale  images.  The  panels  follow  the  same  arrangement  as  Fig.  23. 

Figure  28.  Measured  and  reconstructed  vertical  velocity  fluctuations  vs.  downstream  distance  in 
the  near  wake;  fluctuations  induced  by  suction  side  vortices  v'l/u]  ; 

fluctuations  induced  by  pressure  side  vortices  v'2  ju \  ;  contribution  of 

correlation  between  suction  and  pressure  side  induced  fluctuations 


48 


2(v,v  ~v,vp)/U* - ;  total  fluctuations  from  identified  vortices,  v'2/u20 - ; 

and  PlV-measured  variance  of  the  vertical  velocity  fluctuations  v'2/t/02  • .  The 
panels  follow  the  same  arrangement  as  Fig.  23. 

Figure  29.  Illustration  of  the  vertical  velocity  induced  at  the  trailing  edge  by  a  suction  side 
vortex.  Here  the  induced  velocity  is  assumed  to  influence  the  formation  of  a 
pressure-side  vortex.  The  suction  side  vortex  strength  may  be  related  to  the  suction 
side  circulation  flux  across  surface  A-B. 

Figure  30.  The  normalized  mean  suction  side  shear  rate  at  the  trailing  edge  vs.  Rec  from 
Geometry  I  (O)  and  Geometry  I  (•). 

Figure  31.  Vortex  shedding  characteristics  for  the  trailing  edge  flows  of  this  study  and  for  a 
minimally  lifting  strut  (Blake,  1 984)  vs.  the  dimensionless  shear  rate  parameter, 
{dUjdh)’ l\j0/ Ay  )•  Shedding  parameters  (y-axes)  are  extracted  from  the  LDV 

temporal  spectra  at  x/C=  1.01  (see  Figs.  16,  19,  and  20).  The  mean  suction  side  shear 
rate  is  from  Fig.  30  but  is  scaled  by  the  measured  wake  thickness,  Ay ,  given  in  the 

text  to  form  the  x-axes  in  this  figure.  The  grey  triangles  indicate  strut  data:  strut  with 
45°  rounded  trailing  edge  A;  25°  rounded  trailing  edge  V ;  and  strut  with  25°  sharp 
trailing  edge  ^ .  The  solid  lines  indicate  curve  fits  in  each  panel:  y=0.001+l(T  e 
{z),y= 1.02-0. 77x  (b),y=0.023+0.29x6  (c),  and y=-0.039+0.076x  (d). 
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